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Abstract: The asymptotically- logarithmically-AdS black-hole solutions of 5D dila- 
ton gravity with a monotonic dilaton potential are analyzed in detail. Such theories 
are holographically very close to pure Yang-Mills theory in four dimensions. The 
existence and uniqueness of black-hole solutions is shown. It is also shown that 
a Hawking-Page transition exists at finite temperature if and only if the potential 
corresponds to a confining theory. The physics of the transition matches in detail 
with that of deconfinement of the Yang-Mills theory. The high-temperature phase 
asymptotes to a free gluon gas at high temperature matching the expected behavior 
from asymptotic freedom. The thermal gluon condensate is calculated and shown to 
be crucial for the existence of a non-trivial deconfining transition. The condensate 
of the topological charge is shown to vanish in the deconfined phase. 
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1. Introduction and Summary 

In the past decade we have witnessed a rebirth of theoretical efforts to address the 
strong coupling problem in gauge theories. The tool has been the large- N c expansion 
but a new twist has emerged, whereby the relevant dual string theories live in an 
appropriately curved higher dimensional space-time. 

The prototype example has been the AdS/CFT correspondence as exemplified 
by the (well studied by now) duality of H = 4 super Yang-Mills theory and IIB 
string theory on AdS 5 x S 5 [|I], |2], ||] . Further studies focused on providing examples 
that are closer to real world QCD, []|, It is fair to say that we now have a good 
holographic understanding of phenomena like confinement, chiral symmetry and its 
breaking as well as several related issues. The finite temperature dynamics of gauge 
theories, has a natural holographic counterpart in the thermodynamics of black-holes 
on the gravity side, and the thermal properties of various holographic constructions 
have been widely studied, [§], || || [7], |], |10[, exhibiting the holographic version of 
deconfinement and chiral restoration transitions. 

The simplest top-down string theory model of QCD involves D4 branes with 
supersymmetry breaking boundary conditions for fermions |4[], as well as a flavor 
sector that involves pairs of D 8 — D 8 probe branes inserted in the bulk, ||11|| . The 



qualitative thermal properties of this model closely mimic what we expect in QCD, 
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[^]. Although such theories reproduced the qualitative features of IR QCD dynamics, 
they contain Kaluza-Klein modes, not expected in QCD, with KK masses of the same 
order as the dynamical scale of the gauge theory. Above this scale the theories deviate 
from QCD. Therefore, although the qualitative features of the relevant phenomena 
are correct, a quantitative matching to real QCD is difficult. 

Despite the hostile environment of non-critical theory, several attempts have 
been made to understand holographic physics in lower dimensions in order to avoid 
the KK contamination, based on two- derivative gravitational actions, [12 1 . Indeed, 
large N QCD is expected to be described by a 5-dimensional theory. The alternative 
problem in non-critical theories is that curvatures are of string scale size and the 
truncation of the theory to the zero mode sector is subtle and may be misleading. 

A different and more phenomenological bottom-up approach was developed and 
is now known as AdS/QCD. The original idea described in [Tj| was successfully 
applied to the meson sector in [TJ|], and its thermodynamics was analyzed in || |T5[ |. 
The bulk gravitational background consists of a slice of AdSs, and a constant dilaton. 
There is a UV and an IR cutoff. The confining IR physics is imposed by boundary 
conditions at the IR boundary. This approach, although crude, has been partly 
successful in studying meson physics, despite the fact that the dynamics driving 
chiral symmetry breaking must be imposed by hand via IR boundary conditions. 
Its shortcomings however include a glueball spectrum that does not fit very well the 
lattice data, the fact that magnetic quarks are confined instead of screened, and 
asymptotic Regge trajectories for glueballs and mesons that are quadratic instead of 
linear. 

A phenomenological fix of the last problem was suggested by introducing a soft 
IR wall, Although this fixes the asymptotic spectrum of mesons and meson 

dynamics is in principle self-consistent, it does not allow a consistent treatment of 
the glue sector both at zero and finite temperature. In particular, neither dilaton nor 
metric equations of motion are solved. Therefore the "on-shell" action is not really 
on-shell. The entropy computed from the BH horizon does not match the entropy 
calculated using standard thermodynamics from the free energy computed from the 
action, etc. Phenomenological metrics for the deconfined phase were also suggested, 
[ |T7| , [T3 | capturing some aspects of the expected thermodynamics. 

The theoretical advances were paralleled by a very successful experimental effort 
at RHIC, JH)]. The consensus on the existing data is that shortly after the collision a 
ball of quark-gluon plasma (QGP) forms that is at thermal equilibrium, and subse- 
quently expands until its temperature falls below the QCD transition (or crossover) 
where it finally hadronizes. Relativistic hydrodynamics describes very well the QGP 
[pO[l , with a shear- viscosity to entropy density ratio close to that of Af = 4 SYM, [pi 



The QGP is at strong coupling, and it necessitates a treatment beyond perturbative 
QCD approaches, j2"2]| . Moreover, although the shear viscosity from M = 4 seems to 
be close to that "measured" by experiment, lattice data indicate that in the relevant 



- 4 - 



RHIC range 1< £■ < 3 the QGP seems not to be a fully conformal fluid. Therefore 

J- c 

the bulk viscosity is expected to play a role near the phase transition |23|, [24[. The 
lattice techniques have been successfully used to study the thermal behavior of QCD, 
however they are not easily extended to the computation of hydrodynamic quanti- 
ties. They can be used however, together with parametrizations of correlators in 
order to pin down parameters [23]. To date it seems that the holographic approach 
is a promising one in this direction. 

In the bottom-up holographic model of AdS/QCD, the bulk viscosity is zero as 
conformal invariance is essentially not broken (the stress tensor is traceless). In the 
soft-wall model, no reliable calculation can be done for glue correlators and therefore 
transport coefficients are ill-defined. Similar remarks hold for other phenomeno- 
logically interesting observables as the drag force and the jet quenching parameter 

In order to go beyond the inadequacies of existing bottom-up holographic models, 
input has been put together both from string theory and QCD in order to craft an 



improved holographic QCD model, ||28|| . It is a five-dimensional Einstein dilaton 
system, with an appropriately chosen dilaton potential. 

The vacuum solution involves an asymptotically logarithmically AdS solution 
near the boundary. The bulk field A, dual to the 't Hooft coupling, is vanishing 
logarithmically near the boundary in order to match the expected QCD behavior. 
This implies that the potential must have a regular Taylor expansion as A —>■ 0, and 
that A = is not an extremum. This is unlike almost all asymptotically AdS solutions 
discussed so far in the literature. In particular the canonically normalized scalar (the 
dilaton) is diverging at the boundary r — > as <p ~ — log(— logr). The coefficients 
of the UV Taylor expansion of the potential are in one-to-one correspondence with 
the holographic /3-function. 

In the IR, the potential must have an appropriate behavior so that the theory 
is confined, has a mass gap and a discrete spectrum. This selects a narrow range of 
asymptotics that roughly obey 

V(X) ~ X 2Q , A -> oo with -<Q <-. (1.1) 

3 3 

The vacuum solution always ends in a naked singularity in the bulk. Demanding 
that this is a "good" singularity in the classification of Gubser [f29| implies Q < 4/3. 
However, here we use a narrower notion of what we mean by "good" singularity: we 
accept only singularities that are repulsive to physical fluctuations, i.e. such that no 
extra boundary conditions are needed there. This requirement further restricts Q < 
in Ql.lj) . Simple interpolations between the UV and IR asymptotics reproduce 
very well the low-lying glueball spectrum as well as the perturbative running of the 
't Hooft coupling p8] . 



Five-dimensional Einstein dilaton systems with a monotonic dilaton potential (no 
extrema) provide an interesting class of gravitational theories that display diverse 
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behaviors as a function of the IR asymptotics of the potential. Potentials with 
asymptotics growing slower than ( |1 . 1|) with Q = 2/3 do not exhibit confinement. 
Potentials with asymptotics growing faster than (|0|) with Q = 2^/3 do exhibit 
confinement but the naked singularity is too strong and extra boundary conditions 
are needed at the singularity in order to study the spectrum of fluctuations. 

In this paper we will analyze the existence and structure of black-hole solutions, 
and their thermodynamics, in the class of gravitational theories described above. We 
will take the horizon to be a flat three-dimensional torus, but our techniques extend to 
the case where the horizon is a three-sphere. Our preliminary results in this direction 
have been published in []30|]. The thermodynamics of similar systems has also been 
analyzed in |31|. Our aim is to eventually describe the finite-temperature structure 
of a holographic model closely resembling pure large-iV c Yang Mills theory. The 
latter is widely analyzed using lattice techniques (see e.g. []32], |33|), which indicate 
that the theory exhibits a first order confinement-deconfinement phase transition at 
a non-zero critical temperature of the order of the strong coupling scale A. While 
one of the main motivations for having a realistic holographic description of finite- 
temperature QCD is the computation of transport coefficients and other quantities 
relevant for heavy ion collision experiments, in this paper we will only be concerned 
with equilibrium thermodynamics, as this is already a daunting task. We leave the 
hydrodynamics for the very near future |34j . 



As we show in this paper, the correspondence between the 5D Einstein-dilaton 
setup advocated in and 4D large N c pure Yang-Mills extends to the finite tem- 
perature regime in a remarkable way. One of the main results of this work is the 
proof of the the existence, in the case of confining theories, of a first order Hawking- 
Page phase transition between the thermal gas and black-hole solutions. Moreover, 
the 5D black-holes in confining theories provide a realistic holographic description of 
the thermodynamics of the deconfined phase of 4D, large N c pure Yang-Mills, that 
emerges from lattice studies. 

The black-holes that we discuss obey the same UV asymptotics as the zero- 
temperature solution, namely they are asymptotically-logarithmically AdS with a 
logarithmically vanishing dilaton. Close to the AdS boundary the metric is that of 
a 5D AdS-Schwarzschild black-hole in Poincare coordinates (i.e. with flat horizon), 
up to logarithmic corrections. Although asymptotically AdS black holes in Einstein- 
dilaton theories have received considerable attention (see e.g. P5|), these examples 
were always associated with the existence of an exactly AdS solution with constant 
dilaton, corresponding to an extremum of the dilaton potential. In contrast, the 
5D black- holes we discuss here are of a new type, since the dilaton potential in our 
model is always monotonic. The AdS point is at infinity in field space, therefore the 
models we discuss do not have a pure AdS solution. 

In this paper we derive a series of general results about the thermodynamics of 
the 5D system, that do not depend on the specific form of the potential but only 
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on its small A and large A asymptotics. All these results point to the remarkable 
similarity between the thermodynamics of 5D models in the confining class, and the 
thermodynamics of 4D large N c Yang-Mills. While the detailed quantitative com- 
parison between a specific model and the lattice results for thermal Yang-Mills will 
appear elsewhere [|36|| , these general results are the main focus of this paper. Below 
we briefly summarize them. 



Parameters of the solutions 

The parameters of the action are the 5D Planck mass M p , and the various parameters 
that determine the shape of the dilaton potential V(A). In particular, the value V(0) 
sets the AdS length scale I. In addition, every black hole solution is characterized by 
the five integration constants of the 5th order dilaton-gravity system of equations. 

We show how to identify these parameters in the dual gauge theory and eventu- 
ally how to fix them. 

• We keep the form of the potential generic, except that it should be a monotonic 
function of A and it should obey the small and large A asymptotics we mentioned 
above. A specific potential can be fixed by requiring that the zero-T spectrum 
agrees with the lattice data, as was done in |28fl . 

• The specific value of the AdS length i is irrelevant for any physical observable, 
and sets the overall units of the 5D solution. 

• The 5D Planck scale (in AdS units) is fixed by matching the free field asymp- 
totics of the QCD free energy in the high-temperature limit. This universally 
fixes M p £ = l/(457r 2 ). The physical Planck mass, that determines the strength 
of gravitational interactions (and of the interactions between glueballs in the 
dual theory) includes an extra factor of N^ 3 , that guarantees the suppression 
of quantum corrections in the large N c limit in our setup. 

• Among the 5 integration constants of the equations of motion, only two are 
physical and independent: the value of the dilaton at the horizon, and an overall 
scale A, related to the dilaton asymptotics near the UV boundary. The former 
determines the black hole temperature and entropy; the latter is also present 
in the vacuum solution, and it is dual to Aq C d- Fixing the UV asymptotics is 
equivalent to selecting a specific value for A. 

Since our model should be thought of as coming from a non-critical string theory, an 
additional parameter is the string scale £ s . This is invisible in the 5D Einstein-dilaton 
setup we assume throughout this paper, and we will not discuss it any further. It can 
be fixed by comparing the effective string tension (calculated from the linear part of 
the static quark potential) to the lattice data, as was done in [E3|. 
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Existence and uniqueness of the black-hole solutions 



The T = solution defines a vacuum background. Once we specify the UV asymp- 
totics of a black-hole solution to be the same as for the vacuum background, there 
is only one additional parameter that characterizes the solution, that we take to be 
the value of the dilaton at the horizon, A^. For any monotonic dilaton potential 
V(X) we show that, for each X h , ranging between zero and infinity, there exists a 
unique black-hole solution, with a temperature T, entropy S (horizon area) and free 
energy T functions of Xh only. Thus, the thermodynamics is encoded in the func- 
tions T(\ h ), S(Xh) and JF(A/i), namely the temperature, entropy and free energy. 
Moreover, we show that the limit A^ — » oo of the black-hole solution coincides with 
the unique zero-temperature solution that, for a given potential, displays a "good" 
(i.e. repulsive) singularity. 

Although our models allow an infinite number of black-hole solutions (each with 
a different value of A) with the same mass, this does not imply that these black- 
holes admit scalar hair. The reason is that each different value of A corresponds to 
a different asymptotic for the metric and dilaton. In other words, A plays the role of 
an extra "charge" that can be measured at infinity. Moreover, due to the absence of 
extrema in the dilaton potential, there is no pure ArfS'-Schwarzschild solution with 
a constant dilaton. Thus, our black hole solutions satisfy the no-hair theorems for 
asymptotically AdS gravity coupled to scalars (see |37j for a recent discussion). 



Deconfinement phase transition 

We show that any 5D theory, whose zero-temperature solution is dual to a UV-free 
and IR-confining 4D theory, also exhibits a Hawking-Page type of phase transition, 
dual to a deconfining phase transition in 4D. The transition is always first order, 
except in the special case when the IR behavior of the vacuum solution is at the 
borderline between confining and non-confining geometries: in this case the (string 
frame) solution is asymptotically a flat metric with a linear dilaton, and the phase 
transition is second order. Conversely, non-confining theories do not exhibit a phase 
transition: black-holes dominate the thermal ensemble for any non-zero temperature 
T. 

The phase structure of black-holes in confining theories is similar to what is 
found in the original Hawking-Page situation |R3|, namely: 

• Black hole solutions exist only above a certain temperature T m i n ; 

• Generically, for any T > T m in there exist (at least) two different black-hole 
solutions with the same T and the same UV asymptotics (large and small 
BHs); 
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• Above a certain critical temperature T c > T min , it is (one of the) large black- 
holes that dominate the thermal ensemble, while for < T < T c it is the "ther- 
mal gas" solution (with the same metric and dilaton as the zero-temperature 
solution) that dominates. On the other hand, the small black-holes never dom- 
inate the ensemble. 

• Typically, the big black- holes have positive specific heat, and are thermody- 
namically stable, whereas the small black-holes have negative specific heat and 
are unstable. There may be however exceptions to this rule, where in a limited 
range of the small black-holes are also stable. 

• In the borderline case (asymptotically linear dilaton), there is only a single 
black- hole for T > T m , n , and the second order transition occurs exactly at 
T • 

-*■ mm- 

There is one important difference with the Hawking-Page case, however. There, the 
black- holes are global AdS-Schwarzschild with horizon topology S 3 , and the theory 
is dual to a conformal field theory on the S 1 x S 3 [Q; Here instead we are dealing with 
black- holes whose horizon has topology T 3 , and the phase transition is associated 
with dynamical confinement instead of the non-trivial topology of space. 



Similarities of the black-hole phase and the deconfined phase in Yang- 
Mills 

The thermodynamics of black-holes in 5D duals of confining theories shares many 
features with the deconfined phase of 4D Yang-Mills at large N c . 

• The appropriately regularized free energy TjN 2 c acts as an order parameter for 
the phase transition. This is similar to the case of = 4 SYM on the sphere, 
and as expected in pure YM in flat space; 

• Another order parameter is the Polyakov loop which vanishes in the confined 
phase and becomes non-trivial above the deconfinement transition. This is 
paralleled by the dual gravity computation. A string worldsheet that encircles 
the Euclidean time direction and extends in the radial direction has infinite 
action in the confining geometry, hence the vev of the loop vanishes. On the 
other hand it becomes finite in the black-hole geometry yielding a finite value 
for the associated vev[|J. 

• The latent heat per unit volume is of order N^Tf] 

• At very high temperature the thermodynamic quantities behave like in a con- 
formal theory, although the approach to conformality is logarithmically slow. 
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With a suitable choice of the relation between the 5D Planck scale and the 
AdS length, in the limit T — > oo we find a free gas, as appropriate for a gravity 
dual of pure Yang-Mills and unlike strongly coupled theories like M = 4 sYM . 

• The speed of sound is small near the phase transition, and it approaches the 
conformal value c 2 s — > 1/3 at high temperature 1 . 

The topological vacuum density 

The vacuum expectation value of the topological density, (TtF A F) can be com- 
puted by including, on the gravity side, a 5D axion, dual to the YM theta parameter 
[ffijfl. We show that in the black-hole phase (deconfinement) the profile of the axion 
is necessarily trivial, unlike the low temperature phase (confinement). This causes 
the vev of the topological density to vanish. It is in agreement with the large N c 



expectations and with the lattice calculations in finite temperature Yang-Mills [41 



Explicit calculation of the free energy and the role of the gluon conden- 
sate 

We compute explicitly the free energy of the black-hole solutions, relative to the vac- 
uum, as the difference between the on-shell actions. A crucial role in this computation 
and in the dynamics of the phase transition is played by (the holographic dual of) the 
thermal vev of TrF 2 . This quantity appears in the near-boundary expansion of the 
difference between the black-hole metric scale factor b(r) and its zero-temperature 
analogue, b Q (r). If r denotes the conformal coordinate of both metrics, the AdS 
boundary is at r = 0. We show that, once the UV asymptotics are fixed to be the 
same for all the solutions, then as r — > 0: 

b(r) - b (r) = b a (r) + subleading^ . (1.2) 

The quantity Q is proportional to (TrF 2 )r— (TrF 2 ) D . It is a function of temperature, 
and it provides a measure of the deviation from conformality. It plays a crucial role 
for the existence of the phase transition. 

Indeed, an explicit calculation shows that the free energy difference between a 
black-hole and the vacuum solution is given schematically (omitting 3-space volume 
and other numerical factors) by: 

r~g--?-. (i.3) 



1 Reproducing this behavior was the main motivation of pi"] , and it emerges quite naturally in 
our setup 
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The second term is negative definite, and it is the only one present in a conformal 
field theory. The gluon condensate is therefore crucial for the existence of a phase 
transition. 

One important point is that the free energy written above receives contributions 
from the UV boundary alone: to compute T it is sufficient to know the metric close 
to the AdS boundary. This is important for at least two reasons: 

1. In the context of AdS /CFT all the information about observable quantities 
must be encoded in near-boundary data. No explicit contributions are obtained 
from the IR regime. The IR influences only indirectly by fixing normalizable 
modes near the boundary via regularity conditions. 

2. The specific dynamics of the IR singularity is irrelevant for the computation 
of the thermodynamics. By contrast, in other studies claiming the existence of 
deconfining phase transitions in simple 5D models || [|, it is the IR boundary 
or singularity that gives the required positive contribution to the free energy. 

Matching the (3 -function to the trace anomaly 

We provide a non-trivial check of the gauge/gravity correspondence applied to the 



5D theories of : the matching of the (flat space) conformal anomaly, encoded in 
the YM /3-function, to lowest order in a small A expansion. 2 

From the free energy, we can compute the trace of the thermal stress tensor in 
the deconfined phase, which turns out to be proportional to the gluon condensate: 

(Tg) ~ 0. (1.4) 

On the other hand, in 4D Yang-Mills theory the same quantity (Tjf) obeys the 
dilatation Ward identity: 

TO = ^<TrF 2 >. (1.5) 



By a holographic computation we can find the relation between (TrF 2 ) and Q (the 
latter being defined by eq. (|1.2)), to lowest order in the A — > limit, and can show 



that the two expressions for coincide precisely. 



Part of the analysis in this paper is performed with the help of some new technical 
tools that we believe are interesting by themselves: 

• The thermal generalization of the superpotential: This is widely used in the 
zero-temperature counterpart (see e.g. ]39], ||(J). A superpotential W allows 



to recast Einstein's equations for the dilaton and scale factor in the first order 
form, and to decouple them from the equations governing the evolution of the 
thermal function appearing in the metric; 



2 



This is the initial and simplest step in a rigorous program for the renormalization of asymptot- 
ically logarithmically AdS theories, that will appear shortly, |49[. 
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• The scalar variables: this is a pair of functions of A, X and Y, that are invariant 
under radial diffeomorphisms. They satisfy a coupled system of first-order 
differential equations. These functions encode all the information about the 
UV and IR properties of the full solution. From them one can easily derive all 
the thermodynamic observables and relations. 

The paper is organized as follows. In Section we review the setup, the vacuum 



solutions, and the results about confinement found in ||28|| . We analyze the possible 
types of singularity and we give a more exhaustive analysis of this issue as compared 
to |28 . 



In Section 3 we describe the black-hole solutions their existence and uniqueness 
properties and define the relevant thermodynamic quantities. We then compute the 
free energy difference between the black-hole and vacuum solutions, as a function of 
entropy, temperature, and the value of the gluon condensate. 

In Section f| we show that, to lowest order in A as A — > 0, the trace anomaly 
computed from the equation of state matches the holographic computation of the 
the vev of TrF 2 . 

In Section |5| we prove that confinement at zero temperature is in one-to-one cor- 
respondence with the presence of a phase transition at a finite temperature T c . We 
find the explicit form of the black-hole solutions in the two opposite regimes when 
the black-hole horizon is very close to or very far from the UV boundary. We then 
show that in confining theories there is always a finite, minimum black-hole tem- 
perature, whereas in non-confining theories black-holes exist with arbitrarily small 
temperatures. This fact, together with the first law of black-hole thermodynamics, is 
used to prove the main statement of this section in the particular case when there are 
only two black-hole solutions for each temperature. The proof in the most general 
case is left to Appendix G. 

In Section 6 we study the dynamics of the 5D axion, dual to the Yang-Mills 
vacuum angle, showing that above the critical temperature the axion profile is nec- 
essarily constant, and the topological density has zero vev. 

In Section 7 we define the scalar variables, and show that their use helps in 
computing all the thermodynamic quantities, as well as the UV and IR asymptotic 
properties of the solutions. In particular, in Subsection 7.2 we compute the near- 
boundary expansion of the black-hole metric components and dilaton profile, with 
respect to the vacuum solution. Section 8 contains a brief outlook. 

Most technical details are left to the Appendices. In Appendix A we give details 
about Einstein's equations in various frames, and the relation between fluctuations 
in different frames. In Appendix B we revisit the case of a constant potential and 
derive AdS and dilaton flow solutions and the corresponding black-holes. In Ap- 
pendix C we give the details of the computation of the black-hole on-shell action 
and ADM mass. In Appendix D we provide the high-T asymptotics of the gluon 
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condensate. Appendix E is devoted to the discussion of the general solution to the 
zero-temperature superpotential equation, and the classification of zero-temperature 
singularities. In Appendix F we introduce the finite temperature generalization of 
the superpotential, and use it to solve explicitly the black-hole equations for large A^. 
In Appendix G we give the general proof of the statement in Section 5, in the case 
where more than 2 black-holes exist for certain temperature ranges. In Appendix H 
we provide some details of the computations with scalar variables. In particular in 
Subsection H.4 we show that a black-hole solution with a regular horizon always con- 
nects with the UV boundary, thus providing the proof of the existence of black-hole 
solutions for arbitrary A^. In Appendix I we give the high-T asymptotics of various 
quantities. Finally in Appendix J we show some interesting analytical solutions of 
the system. 



2. Review of vacuum solutions 

The holographic duals of large N c Yang Mills theory proposed in PS[ are based on 



five-dimensional Einstein-dilaton gravity with a dilaton potential. The basic fields 
for the pure gauge sector are the 5D metric (dual to the 4D stress tensor) and a 
scalar field <3> (dual to TrF 2 ). The action for these fields is taken to be 3 : 

S 5 = -M*N 2 [ d 5 x^/g R - \{d§) 2 + V($) + 2M*N 2 / d 4 xVh K. (2.1) 
J L 3 J J aM 

Here, M p is the five-dimensional Planck scale and N c is the number of colors. The 
last term is the Gibbons- Hawking term, with K being the extrinsic curvature of the 
boundary. The effective five-dimensional Newton constant is G5 = l/(167rMp N 2 ), 
and it is small in the large- N c limit. 

The vacuum solutions are of the form 

ds 2 = b{r) 2 (dr 2 + r^cfeW) , $ = $(r), (2.2) 

where the metric is written in conformally-flat coordinates. 

The radial coordinate r corresponds to the 4D RG scale. In the holographic 
dictionary, we identify the 4D energy scale E with the metric scale factor, E = E Q b(r), 
up to an arbitrary energy unit Eq. Also, we identify A = e* with the running 't Hooft 
coupling At = N c gy M , up to an a priori unknown multiplicative factor, A = n\ t . All 
physical observables are independent of the parameter k, as explained in Appendix 

& 

3 See Appendix A for our sign conventions. 

4 More precisely, this statement applies to quantities that can be computed within the Einstein 
frame. Quantities that involve the string frame metric may depend on re in a nontrivial way. In this 
paper we will not be concerned with any such quantity, and we leave a more detailed discussion of 
this issue for an upcoming work [p6[. 
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With the above identifications, one can give a holographic definition of the (3- 
function of the system in terms of the background solution: 



Above and throughout this paper a dot stands for a derivative with respect to the 
radial (conformal) coordinate r. 

With the ansatz (|2.2[), Einstein's equations are 



The dilaton field equation is not an independent equation, but it follows from (|2.4j ). 

It is sometimes useful to work with the domain wall coordinates, in which the 
metric reads: 

ds 2 = du 2 + e 2A(u) 7 ]ij dx l dx j , dr = e~ A{u) du, b{u) = e A{u \ (2.5) 
In this frame Einstein's equations take the form: 

3A" + 12A' 2 = V($), A" = (2.6) 

9 

where a prime denotes a derivative w.r.t. u. In particular, the second equation 
implies that the scale factor of an asymptotically AdS§ spacetime (for which A ~ 
—u/£ as u — > — oo) is monotonically decreasing, and it therefore provides a consistent 
definition of the holographic energy scale. 

Einstein's equations can be put in first order form by defining a superpotential 
W($), i.e. one solution to the equation: 



With this definition, Einstein's equations ( |2.6| ) become: 

$'(«) = A'(u) = ~W(9). (2.8) 

The system of eqs. ( |2.7| - |2.8| ) has three integration constants, one of which is an 
artifact due to reparametrization invariance 5 . The solution is completely specified 
by a choice of W($>), up to an integration constant that consists in a simultaneous 
rescaling of b(r) and r, and only affects the overall scale of the system. In other 
words, all nontrivial physics is encoded in VT(<I ) ). 



5 This can be seen by choosing $ as a coordinate: then the first equation in ( |2.8| ) becomes 
vacuous, and only two first order equations remain. 
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The general solution of eq. (|2.7|) is discussed in detail in Appendix [E]. As we will 
discuss at the end of this section, for any V($), there is a single choice of PU($) that 
satisfies some reasonable physical conditions. We are thus left with a one-parameter 
family of solutions, distinguished only by a choice of scale. This parallels the situation 
in the gauge theory. 

Another useful reformulation of the Einstein's equations is in terms of the loga- 
rithmic derivative of W($), which is directly related to the /3-function: 

xw = -Y^l = m ( 2. 9 ) 

The complete solution of the system is encoded in this function. It is determined 
from the potential by solving a first-order equation: 

d®-~3 {1 - X) { 1+ 8X^T ) ■ (2 ' 10) 

Once X is known the scale factor and the dilaton are obtained from it by solving the 
first order equations ( |H.1| ) and ( |H.2| ). Thus, this formulation reduces the Einstein 
equations to three first order equations. In Section |7|, we shall present a natural 
generalization of this formulation to the black-hole solutions. 
The precise relation between A($) and VT(<I?) is given by, 

W($) = ^e-*J-°° xm , (2.11) 

where I is the asymptotic AdS length (see next subsection). This relation holds both 
for the zero T theory and at finite T. In the latter case, the functions X and W are 
replaced by their finite T counterparts, as explained later. 

The small-A and large-A asymptotic of W(X) (or X(X)) determine the solution in 
the UV and the IR of the geometry, corresponding to the large- and small-6 regions, 
respectively. 

2.1 UV asymptotics 

In the UV, asymptotic freedom with logarithmically running coupling requires the 
background to be asymptotically Anti-de Sitter. The perturbative /^-function, (3 ~ 
—b$X 2 — biX 3 + . . ., requires an expansion of X in the form: 

V(A) = -|a-|a 2 + G(A 3 ) (2.12) 

where bo and b\ are the /3-function coefficients. Using ( |2.9|) one finds the expansion 
of the superpotential as, 

W{X) = ^ (1 + w X + Wl X 2 + ...), (2.13) 
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which implies a potential of the form 



1 2 

V(X) = -(l + v X + v 1 X 2 + ...). 



(2.14) 



Here i is the AdS length, and the dimensionless parameters Wi, u$ are fixed in terms 
of the /^-function coefficients. In particular the small-A expansion parameters Wi 
of the superpotential are universal, and do not depend on the particular choice of 
solution of eq. ( |2.7| ): as shown in Appendix [E], different solutions of eq. ( |2.7| ) differ 
by subleading non-analytic terms. For a general potential (|2.14|) , the /3-functions 
coefficients and the parameters of the potential are related as follows: 



9 



-v 



9 



bi 



9 



-Vi - 



207 

256 ? 



9 



9 

-Wi — -w, 



o- 



Let us here also define the ratio, 



b = I2> 



hi 



(2.15) 



(2.16) 



which will prove useful in what follows. Note that b is invariant under the rescaling 
A — > k\. 

The UV region corresponds to r — > in conformal coordinates, and the asymp- 
totic solution is given by: 

1 



b(r) = 
M(r) 



9 log r A 9 



4 b log (-log rA) 



log rA 



+ 6 



log(— logrA) 



+ . . . 



(2.17) 
(2.18) 



logrA ' log 2 rA 

The scale A appearing in the expansion is the only physical integration constant, 
and it is the holographic manifestation of the strong coupling scale in QCD perturba- 
tion theory. The UV boundary conditions for the metric and dilaton are completely 
specified by the choice of this scale. In practice A is determined by a combination of 
the initial conditions of A and A, given at a point r close to the boundary, as, 



A 



cxp 



^4(Ao) 



1 



Mo 



+ ■ 



(2.19) 



The ellipsis refer to contributions that vanish as one takes the cut-off away, Ao — > 0. 
The coefficients bo and b are defined in ( |2.15| ) and ( |2.16| ). The derivation of ( |2.19j ) 
follows from Appendix [H.l| and the eqs. ( p,17| - ^T^ ) above. 



2.2 IR asymptotics 

The IR properties such as confinement of the electric color charges (signaled by an 
area law for the Wilson loop) and the features of the glueball spectrum are determined 
by the behavior of jy(A) (or X(X)) for large A. In particular, the Wilson loop follows 
an area law if and only if Vy(A) grows as A 2 / 3 or faster. The same condition ensures 
a mass gap in the spectrum. In other words one has the criterion: 
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Confinement W(X) > 0(A 2/3 ) as A — oo. 

The form of the IR geometry depends on the details of the asymptotics. The 
Einstein-frame scale factor is guaranteed to decrease monotonically from the UV to 
the IR, and eventually the spacetime terminates in a singularity at some r = r . We 
classify the singularity into good and bad according to the following criterion [J2^J : 

A good singularity is screened, i.e. it is repulsive to physical modes. 

On the other hand, bad singularities are such that finite energy modes can probe 
arbitrarily deep into the region close to the singularity. Typically this means that one 
needs to specify extra boundary conditions at the singularity, i.e. the information 
provided with the classical action is not enough to compute physical quantities. 
For good singularities, all (physical) boundary conditions must be imposed in the 
UV region. Therefore we believe that only "good" singularities have a meaningful 
holographic interpretation. 

The most interesting geometries are those with the singularity at r = oo, and 
with the asymptotics: 

b( r ) ~ e -U) , A(r) ~ e 3/ Hr) (_J ; r — oo (2.20) 
Here, the length scale L is set by the same integration constant that fixes A in eq. 



In such solutions the curvature of the string-frame metric vanishes in the extreme 
IR. These solutions occur when W(X) and X(X) behave for large A as: 

W/(A)~A 2 / 3 (logA)^\ X(A)~-i-§— + A -oo, (2.21) 

2 8 a log A 

which in turn requires the potential to grow as 

V(X) ~ A 4/3 (logA)^, A — oo, (2.22) 

These solutions are confining iff a > l 6 . 

The parameter a determines the asymptotic spectrum of normalizable fluctua- 
tions around the solution, which corresponds to the spectrum of composite states 
(glueballs) of the gauge theory, with masses that scale as: 

m n ~ n (a " 1)/Q . (2.23) 

For a linear glueball spectrum, m 2 ~ n, one should choose a = 2. 



^Solutions such that b(r) decays as a power-law as r — > oo are not confining. 
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The borderline confining case, a = 1 has interesting properties: the asymptotic 
geometry in the string frame reduces to flat space with a linear dilaton. The spectrum 
has a mass gap and it is discrete up to a certain energy level, above which it becomes 
continuous. We will see that this case also has special thermodynamic properties. 

Solutions with a singularity at a finite value r$ of the conformal coordinate are 
also confining, and correspond to W(X) growing as X Q with Q > 2/3. Close to the 
singularity r = r the scale factor vanishes as 



b(r) ~ (r - r) s , Q = + (2.24) 

The glueball spectrum has quadratic growth, m 2 n ~ n 2 , as in the hard wall models. 

The case 5 < 1 (Q > 2\/2/3) should be discarded, since in this case the singu- 
larity is a bad one according to the our classification, i.e. it is not screened from the 
physical fluctuations |28| . 

An example of this type (with 5 = 1/3) is the "dilaton flow" solution of 5D 



Einstein-dilaton gravity with a negative cosmological constant, discussed in p2| , [43 
which was argued to be dual to an 5*0(6) invariant mass deformation Af = 4 SYM. 
Although this description can be adequate in the UV, calculation of any physical 
quantity requires extra knowledge about the details of the singularity, which is not 
available in the classical gravity approximation. 

2.3 The superpotential vs. the potential 

The action fl2.1|) is defined in terms of V(X), not W(X) or A(A). Therefore it is 
important to know what other large-A asymptotics for W and A (A) can occur for a 
given V(X). This problem is analyzed in Appendix [E] (for W) and in Appendix H.3 



(for X), where the form of the general solution of eqs. ( |2.7| ) and ( |2. 10| ) is discussed 
in detail. Essentially, for any given monotonic ^(A), a solution with a good infrared 
singularity, if it exists, is unique. 

In the UV region, A — > 0, all solutions to eq. ( |2.7| ) have the same expansion, given 
by eq. ( |2.13|) with the same coefficients u>j. In other words, in the UV all solutions 



to Einstein's equations flow to the same log-corrected AdS (eqs (|2. 17| - p. 18f) ) . 

In the IR, the situation is more complicated. We consider a potential V(X) 
defined over the whole range < A < oo, and such that for large enough A it is well 
approximated by the form: 

V(X)^V 00 X 2Q (\ogX) p (2.25) 

for some real P and Q 7 . We will assume V(X) is a positive, monotonic function, to 

7 Although we parametrize the IR asymptotic in this particular form, all our discussion also 
applies to any potential that has intermediate IR growth between two values of Q or P that share 
the same behavior, for example V^(A) ~ A^(log A) p (loglog X) R . . .. However if Q or P take values 
that mark the boundary between two different qualitative behaviors of the solution (e.g. Q = 4/3, 
see below), extra work is needed to understand the intermediate asymptotics. 
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avoid the presence of conformal fixed points at finite A. Thus, we take Vx > and 
Q > 0. 

All the interesting physics is found for Q < 4/3. As shown in Appendix ||, if 
Q < 4/3, there exist three classes of solutions to the superpotential equation: 

1. Special: A single solution such that W(X) ~ ^JV(X) for A — > oo. 

2. Generic: A continuous family with leading asymptotics 

W{\) ~ CA 4/3 A -f oo (2.26) 
where C is an arbitrary constant. 

3. Bouncing: A continuous family which never reaches the asymptotic large-A 
region: the variable A attains a maximum value A*, then decreases again to 
zero towards a region where 

W ~ C\- 4/3 , A -> 0, (2.27) 



On the other hand, if Q > 4/3 only the bouncing solution exists 8 , and the dilaton 
never reaches infinity. 

The special solution has asymptotic behavior: 



W (\)*W eo \ t >Qog\) p P A-,00, ^oo = y 4(1 f_ F g g2) . (2.28) 

It presents a good IR singularity for Q < 2^/2/3, and it is confining for Q > 2/3, or 
Q = 2/3 and P > 0. The generic and bouncing solutions, on the other hand, always 
have bad singularities. 

As we shall discuss in Section 5, the special solution is also the only one that can 
be obtained in the zero-mass limit of black-hole solutions of the same bulk theory. 
This gives another characterization of the special solution, and singles it out as the 
only physically sensible choice |29|] . 

Thus, given a potential V(A) with asymptotics fl2.22|) , there is a single solution 



with the large-A behavior (|2.21| ) corresponding to a "good" singularity. All other so- 
lutions have bad singularities in the IR, and cannot be lifted to black-holes. Requiring 
the absence of bad singularities is what ultimately fixes the integration constant of 
the W equation in the zero-temperature system, or equivalently, the integration con- 
stant of ([nop. 

As shown in Appendix [li).2|, changing this integration constant adds a pertur- 



bation in the metric and dilaton that goes as r close to the boundary. Thus, this 



'This includes the case when the potential grows faster than (2.25) for any Q, e.g. if V(X) 
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corresponds to changing the expectation value of the corresponding dimension 4 op- 
erator, i.e. TrF 2 . In other words, the integration constant in the superpotential 
controls the value of (TrF 2 ) in the gauge theory and there is a unique value such 
that no bad singularities appear. 9 . 

To summarize: the physically interesting situation when a good solution exists 
and corresponds to confined color, is the case 2/3 < Q < 2y/2/Z. 



3. Finite-temperature solutions and thermodynamics 

We now consider the dilaton-gravity system described in the previous section with a 
good potential according to the aforementioned criteria and study it at finite tem- 
perature T. As usual, this can be implemented by going to Euclidean signature and 
compactifying the Euclidean time (that for simplicity of notation will be still called 
t) on a circle with period (3 = 1/T. This breaks the Poincare invariance of the vac- 
uum to spatial rotations, and allows for a larger class of solutions. According to the 
AdS /CFT prescription, the partition sum is constructed by considering all solutions 
with fixed UV boundary conditions 10 . From now on we will introduce a subscript 
"o" for the quantities related to the zero-temperature solution. 
The thermal solutions are of two types: 

1. Thermal gas solution: this is the same as ( |2.2j ), 

ds 2 = b 2 (r) (dr 2 - dt 2 + dx m dx m ) , $ = $ (r), (3.1) 

except for the identification t ~ t + i/3. It corresponds to a gas of thermal exci- 
tations above the same vacuum described by the original solution, from which 
it inherits all the non-perturbative features (confinement, spectrum, values of 
condensates, etc. ) 



Black hole solutions: they are of the form 

dr 2 



ds 2 = b{r) 



— f(r)dt + dx m dx r 



$ = $( r ) (3.2) 



and are characterized by the presence of a horizon where /(r/j) = 0. This 
implies that such a solution, if it exists, corresponds to a non-confined phase, 



9 This situation has an analogue in the case of constant potential, V(X) = 12/£ 2 : also in this 
case there is a single "good" solution, AdS^ spacetime with constant dilaton. Even in this theory 
it is well known p2[ |4j| that there is a continuous family of solutions with a running dilaton, all 
of which have bad singularities in the interior. This is presented in detail both at zero and finite 
temperature in appendix [§]. 

10 Later in this section we will be more specific about what we mean by "fixed UV boundary 
conditions." 
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since the confining string tension is proportional to M.m r {^J g xx {r)g t t{r)) = 



[ 44| . In the Euclidean version, deconffnement is signaled by a non-zero value 
of the Polyakov loop, as discussed in ||. Since we want to study the theory on 
S 1 x R 3 we consider black-holes with flat horizon topology. 



Notice that in general the functions b{r) and <3>(r) appearing in eq. (|3.2| ) are 
different from their zero-temperature counterparts, and have also a nontrivial 
temperature dependence. 

In the rest of this section we will discuss the features of the black-hole solutions 
to the general system (|2.1| ), and the thermodynamics in the canonical ensemble. 



3.1 5D Einstein-dilaton Black holes 

We require that the solution has the same UV asymptotics as the one at zero temper- 
ature: an AdS boundary at r = where b(r) ~ ijr and e* vanishes logarithmically; 
we have to impose /(0) = 1, so that the black-hole solution (|3.2|) coincides with the 



zero-temperature and thermal gas solutions, (|3.1|) in the UV limit r — > 0. 

A black-hole solution with a regular horizon is characterized by the existence of 
a surface r = r h , where the dilaton and scale factor are regular, and 

f(r h ) = 0, f(r h ) < 0. (3.3) 

The Euclidean version of the solution is defined only for < r < rh- The horizon 
r = rh is a regular surface if Euclidean time is identified as r — > r + 47r/|/(r/ l )|. This 
determines the temperature of the solution as: 

T = J/M. (3.4) 
An 

The independent field equations are: 



b 2 b 4 • , / b 

6^-3- = -^ , i + 3- = 0, (3.5) 

b 2 b bf b 2 , . 

6 w +3 - b +3 l7 = 7 v - (3J) 



Integrating once the second equation of (|3.5|) , we obtain: 

/ = -£> (3-7) 

for some integration constant C. This shows that / cannot change sign. For a black- 
hole, f(r) has to decrease from / = 1 at the boundary to / = at the horizon, 
therefore C > 0. 
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The general solution for / is 

!{r)=i - c Lwf' m 

where we have chosen the second integration constant so that /(0) = l 11 . Setting 
C = and f(r) = 1 we recover the zero-temperature Einstein's equations. 
The quantity C is related to the horizon location as: 

C = T^rV ( 3 - 9 ) 

JO bp) 7 

Note that b(r) is regular in the whole region of integration. We can compute the 
temperature by eq. Q3.4p: 



n 1 4?r t s / \ f rh du A7ib 3 (r h ) 

(3 = - = - 77 — = 4tt6 3 (r h ) / _ = — (3.10) 

The horizon area is given by 

A(r h ) = b 3 (r h )V 3 , (3.11) 

where V3 is the volume of 3-space, and it is related to the entropy as usual by 
S = A/4G 5 . 

In the particular case V($) = 12 /£ 2 , $ = 0, we have the AdS'-Schwarzschild 
solution in Poincare coordinates, 

r \ „ 1 



6(0 = -. /(r) = l- - , T=— , A=[-)V S . (3.12) 
r \r h / nr h \r h J 

Notice that in this case the scale factor is temperature-independent. This is not true 
in general: for V depending non-trivially on different f(r) will result in different 
6(r). 

Near the AdS boundary (UV), this difference can be made more precise: 

• As shown in Appendix [F.3.1| , near the AdS boundary b(r) and A(r) have AdS 
asymptotics, an expansion in inverse logs of the same form as b Q (r) and A (r), 
eqs. ( |2.17| - |2.l8D , specified by an integration constant A. In particular, A(r) — > 



in the UV for all the solutions. Fixing the UV boundary conditions therefore, 
means specifying the scale A appearing in this expansion. 

Here, and from now on, by stating that two solutions obey the "same UV bound- 
ary conditions" , we require that "the scale A appearing in the expansion in the 
perturbative UV log must be the same" 12 . 



11 Recall that b(r) ~ r _1 as r — > 0, so the second term in eq. vanishes at the boundary. 

12 As it should be clear, it is not enough to specify that the metric be asymptotically AdS and 
that X(r) asymptotes some fixed value as r — ► 0, since for all solutions A(0) = 0. 
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• Assuming for b(r) and A(r) the same value of the integration constant A, as 
for b {r) and A (r), then: 



b(r) = b (r) 
A(r) = A (r) 

/(r) 



1 + + 



45 _ r 4 log Ar 
1 + — G J — 



r -> 0, 



Cr 4 
X -4^ + 



0. 



(3.13) 
(3.14) 
(3.15) 



where C is defined in (|3.7|), and £ is a temperature-dependent constant with 
the dimensions of energy. Eq. ( |3.15| ) is obtained from the expression (31?) and 
the fact that b(r) — > £r _1 as r — > 0; eqs. ( p,13| - ^T4|) will be derived explicitly 
in Section 1772]. 



According to the standard rules of the correspondence, the quantity Q(T) is 
interpreted (up to a multiplicative constant, to be determined later) as the difference 
between the vev's of the corresponding dimension-four operator in the black-hole 
and in the vacuum solution. Since we have assumed that $ couples to TrF 2 as 
J e~®TrF 2 , the appropriate operator is the gluon condensate A~ 1 TrF 2 . The precise 
relation between Q and (TrF 2 ) will be obtained in Section f|. 

Integration constants The quantities C and Q are related to two of the five in- 
dependent integration constant of the system of Einstein's equations. For C, this is 
clear from its definition, eq. ( |3.7| ): it determines the temperature of the black-hole. 
Q can be regarded as the integration constant for the thermal generalization of the 
superpotential equation, given in Appendix |F|. In the zero-temperature case it was 
fixed to single out the special solution, with the "good" IR behavior; in the black-hole 
it is fixed by the requirement of regularity of the horizon 13 . Two more integration 
constants are fixed by setting /(0) = 1, and by choosing the scale A in the UV 
perturbative expansion, i.e. by requiring that the solution has the same boundary 
behavior as at T = 0. The last integration constant, as in the T = 0, is unphysical 
and is due to reparametrization invariance. It can be eliminated by rewriting the 
solution using A as a coordinate. As we show in Appendix |F.2| , a better way of 
counting integration constants is by giving "initial values" directly at the horizon. 
This results in the following statement: 



13 In our analysis we cannot determine uniquely what is the value of the gluon condensate in 
a given background, but only the differences between TrF 2 in two backgrounds with the same 
asymptotics. In order to compute the v.e.v. of the gluon condensate unambiguously in a given 
background one would need to perform the full procedure of holographic rcnormalization, which for 
the asymptotics we are considering is not yet fully developed but will be available soon, |49|. 
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For any positive and monotonia potential V(X) that grows no faster 14 than A 4//3 as 
A — > oo, and for any value \h, there exists one and only one black-hole such that: 

1. A — > \h at the horizon 

2. it has the same UV asymptotics as the zero-temperature solution (\2. 1 r l\-^T8J 



The existence, for each A^, of a black- hole solution with regular horizon that 
extends all the way to the UV AdS boundary, is shown in Appendix [H.4| using 



the method of scalar variables. Uniqueness, on the other hand, follows from the 



discussion in Appendix |F.2 



Thus, the value of the dilaton at the horizon, A^, is the most natural candidate to 
classify the black-holes, and all the thermodynamic quantities like e.g. temperature 
and entropy which are single valued functions of A^ (the same does not necessarily 
hold if one writes them as a function of the horizon position, i.e. T(r/j) is not 
necessarily single valued. We have found numerically examples of this behavior). 

3.2 Thermodynamics 

In this section we compute the free energy differences of various solutions at a given 
temperature. This will allow us to compute all other thermodynamic quantities. The 
details of the relevant calculations can be found in Appendices |C.1| and |C.2 . 



The free energy at fixed temperature (3 1 of a given solution is given by: 

/3F = S £ , (3.16) 

where S e is the regularized Euclidean action evaluated on the solution. The action 
needs to be regularized, due to the divergences near the AdS boundary. To achieve 
this, we take 3-d space to be a torus with finite volume V3, and cut-off the radial 
direction in the UV up to a minimum radius r = e > 0, so all the integrals are limited 
below by e. Free-energy differences will be finite (and proportional to V3) as e — > 
since the large-volume divergences do not depend on the detailed solutions but only 
on the asymptotics. 



As shown in Appendix |C.1| , the regularized Einstein action (|2.1|) evaluated on a 
black-hole solution is: 



S e = S e E + S e GH + S e count = 2(3 a 
where we have defined: 



3b\e)f(e)b(e) + -f(e)b*(e) 



+ S, 



count 1 



(3.17) 



a = M?N?V 3 . 



p c 



(3.18) 



14 This restriction is necessary to ensure that (vacuum and black-hole) solutions that extend to 
arbitrarily large values of A exist. 
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The counterterm action S e count is required to make the above expression finite in the 



limit e — > (|5| . As we can see, eq. ( |3.17[) depends solely on the metric evaluated at 
the UV cutoff: the contribution to the bulk Einstein action coming from the horizon 
region vanishes. 

Instead of explicitly calculating S e count) we can define the free energy by subtract- 
ing a given reference background, following the prescription of jKJ. We should take 
as reference the thermal gas background with the same temperature and the same A 
as the black-hole, obtained by setting f(r) = 1, and replacing b(r) with b a (r). How- 
ever, this is correct for the unregularized action that extends all the way to r = 0. 



When we deal with regularized geometries, we must make sure that [46[ 



i) the intrinsic geometry of the 4-dimensional boundary be the same for the two 
solutions, 

ii) the boundary values of scalar field, A and A G are the same. 

To satisfy i), we must demand that the proper lengths of the time circles of the 
solutions (|3.1|) and (|3.2j) and the proper volume of 3-space be the same at r = e. 
Denoting by (3 and V3 the period of the time coordinate and the volume of 3-space 
in the thermal g clS CcLSC, this requirements imply: 



Pb (r) 



cut-off 



Vsb 3 (r) 



V* b 3 



cut-off 



cut-off 



(3.19) 



cut-off 

The condition ii) means that we must require A(e) = A D (e). This actually implies 
that the two backgrounds are characterized by different values A and A of the scale 
that defines the UV boundary conditions. This makes the calculation of differences 
such as 6(e) — b a (e) quite complicated, since one cannot use directly the UV expansion 
( |3.13|) , which relies on the two scales being equal. However, since we are dealing 
with only one scalar field, we can equivalently keep A = A and set the cut-off of the 
background at a different location r = e, such that 15 



A(e) = A (e). 



(3.20) 



Since the cut-off coordinates do not coincide in the two solutions, in the conditions 
( |3.19| ) one now has to evaluate each side of the equality at the appropriate value of 
the coordinate r. Now we can use eq. ( |3.14| ) to determine the needed shift in the 
boundary positions: 



e — e 



45 £ 



-f^ 5(1 ° geA)2 - 



(3.21) 



8 P A D (e) 

The regularized action for the thermal gas background, in the case that the IR 
singularity is of the good type, is given by: 



S e = 2(3a 3b 2 (e)b (e 



+ S, 



count 1 



a 



M 3 N*V 3 . 



(3.22) 



15 



This strategy will not work in the case of multiple bulk scalar fields. 
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As for the black- hole, eq. ( 3.22|) only receives contributions from the metric evaluated 



at the UV cut-off. Indeed, evaluating the Einstein term on shell in general gives an 
extra negative term localized in the IR, of the form Sjr = 2j3abl(r )b o (r ) < 0, where 
r is the position of the singularity in the vacuum background. This contribution 
vanishes exactly, if tq is a good singularity. This is good news, since this means that 
the details of the physics of the singular region are irrelevant for the calculation of 
the free energy (as for other physical quantities [ffifl), and only the UV data matter. 
Similarly a Gibbons-Hawking term at the singularity (which in principle could be 
there) also does not give any new contribution, since on shell it is also proportional 
to bl(ro)b (r ) and it vanishes for good singularities. 16 
Therefore the free energy (difference) T is: 

/3jF=Hrn(<S e -<S e ~)- (3-23) 

In the difference above, it is guaranteed that the contribution coming from the coun- 
terterms exactly cancel even at finite e, since these terms are build out of invariants 
of the induced boundary geometry and the boundary values of the scalar field 17 . 
Therefore this subtraction prescription makes it unnecessary to know the explicit 
form of S coun t. 

Using the results ( |3.17| ) and ( |3.22|) in Q3.23|) , together with the relations ( |3.19| ), 
we have: 



^ = alim \6b 2 (e)^/f{e) 



6(0 V75 



b 2 ( 



+ f(e)b 3 (e) 



(3.24) 



Using the UV expansions (|3.15| - |3.13|) and the relation (|3.21| ) n eq (|3.24|) , and 
taking the limit e — > 0, we obtain the final result for the free energy, 



C 



1. 



T = a[ 15 SCO - T = 15eG(T) - -TS, 



where the entropy S is given by the area of the horizon: 

S 



A = = o c -. 



(3.25) 



(3.26) 



In the second equality we have used eq. (|3.18|) and G 5 = 1/ (WnMpN 2 ), and in the 
third one eq. ( |3.10| ). 



16 For a recent example of a similar study, where the free energy does receive contributions from 
the deep IR, see [||. 

17 This cancelation may not hold in the case of the counterterms that diverge as loge, i.e. those 
that give rise to the conformal anomaly. However, the free energy difference due to these countert- 
erms, if any, is of the order e 4 loge. Therefore these terms do not result in any finite contribution 
to AT 
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The black-hole energy E can be obtained either by the thermodynamic formula 
E = T + TS = —dp(/3J-), or by computing the ADM mass, and consistency requires 
that the two computations give the same result. This is indeed the case: as shown 
in Appendix |U.3|, the black-hole mass is given by 



E = a (l5 g + ^C^J =J r + TS. (3.27) 

The presence of the gluon condensate term in eq. (|3.25| ) is the source of the 
breaking of conformal invariance, since in a conformal theory the relation T = —TS/A 
is exactly satisfied. As we will see in Section [|, the holographic computation of the 
conformal anomaly in flat space matches the field theory result. 



Notice that second term in eq. ( 3.25 ) is negative for any T. It is the presence 



of a non-trivial gluon condensate that may allow a change in sign of the free energy, 
corresponding to a first order phase transition. While the calculation of Q(T), in 
general, is only possible numerically, in Section |5| we will give a general argument to 
determine whether or not a given Einstein-dilaton theory exhibits a phase transition 
at some critical temperature. 

Finally, we note that the Gibbons-Hawking term contributes importantly to 
this expression. A simple calculation shows that the Einstein term in the action 
contributes as Te = —5aQ — T 5/4 whereas the Gibbons-Hawking contribution is 
Fgh — +20aQ. Note that the Einstein term itself is negative-definite, therefore the 
GH term is crucial for the existence of a phase transition. This is unlike the usual 
Hawking-Page phase transition in global AdS, where the change in sign of the free 
energy is due solely to the Einstein term. 



4. The holographic conformal anomaly (in flat 4-space) 

The expectation value of the gluon condensate plays an important role in the ther- 
modynamics of the system we are investigating. As we will see, it can be related to 
the thermal version of the anomalous trace of the stress tensor. In this section we will 
show that a holographic calculation of the trace anomaly in flat space to lowest order 
in A matches the four- dimensional result as advocated in [33]. This is a non-trivial 



check of the validity of the gauge/gravity duality in our setup, in particular of our 
identification of the holographic /3-function. 

In four-dimensional Yang-Mills theory, breaking of scale invariance is expressed 
by the operator equation: 



T; = ^TrF 2 , (4.1) 
where A t is the 4D 't Hooft coupling. 
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Defining the pressure p and energy density p, 

the trace of the thermodynamic stress tensor can be obtained immediately from eq. 

W>b = P - 3p = Q0M 3 p N 2 c G(T), (4.3) 

where we have used the definition of a, eq. (|3.18|) . The left hand side of ( |4.3| ) is 
the trace of the renormalized thermal stress tensor, (T^)n = (T£) — (T^) D , and it is 
proportional to Q ~ (TrF 2 ), in qualitative agreement with ( |4.1| ). 



To check the detailed agreement between eqs. ( |4~lD and (|4.3|) , we need to derive 
the precise relation between Q and (TrF 2 ). In what follows we work to lowest order 



in A. To compute (TrF ) holographically, we use the prescription established in [[47 
we recall that for any fluctuation of a bulk scalar with a canonically normalized 
kinetic term, 

S fluc = \t-' i J ^d 5 x(d X f, (4.4) 
and coupling to a boundary operator 0(x) as 

Scaup = J d 4 xxO, (4.5) 

one can read off the vev of the operator 0(x) from the UV asymptotics of x{ x i r ) : if 
the UV expansion is of the form 

X ~ r A -Xo(x) + r A +Xi(x), (4.6) 

with A + being the canonical dimension of the dual operator, and A_ = d — A + for 
a c/-dimensional gauge theory, then the vev of the operator is given by the formula: 

(0(x)) = (2A + -d) X i(x). (4.7) 

Now, we apply this prescription to the dilaton fluctuation. In our setup the 
dilaton $ is related to the 't Hooft coupling X t by e* = A = ftA t 18 . The coupling of 
the dilaton to the YM field strength is given by 

S coup = -Jd*x ^-Tr F 2 = ~\\ d 4 x e"* Tr F 2 . (4.8) 

Thus, the dilaton fluctuation couples as, 

5S coup = j I d A x 5$ e- $ Tr F 2 . (4.9) 



18 We will keep the multiplicative constant k unspecified, as it will drop out of the calculation, 
i.e. matching of the anomalies does not depend on the value of k. 
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From the bulk action (|2.1| ) we learn that and the canonically normalized 
fluctuation (|4.4j) are related by 19 



X = ( ^(M p £) 3 N 2 ) 5$. ( 1-10) 



Using the relation e* = n\ t we find the coupling 



d 4 x x O = / d 4 x x 



K 



{MplfNl ) 2 Tr F 2 



4A V3 V p 



(4.11) 



The particular dilaton fluctuation we are interested is the difference 5$ = $ — $ D : 
this allows to compute the difference between the thermal and vacuum values of 
(Tr F 2 (0)). Notice that this difference is well defined, contrary to e.g. (Tr F 2 (0)) o , 
which suffers from UV ambiguities. Moreover, it is a purely normalizable fluctuation 
close to the boundary (i.e. it only contains terms like the second one in fl4.6p , since we 
assumed that the black-hole and vacuum backgrounds obey the same UV boundary 
conditions. 

The fluctuation 5$ is obtained from eq. fl3.14j ). To leading order in powers of A 
(or equivalently in inverse powers of log(rA)): 

te~^f?^]og(rA), (4.12) 

Recalling the leading dilaton asymptotics (6oA) _1 = — log(r A), and using Q4.10Q 
and ([4.12 ) together with ( |4.7| ) and ( 4.11 ) yields the relation between Q and Tr F 2 : 



240 

(Tr F 2 ) T - (Tr F 2 ) = —M 3 N 2 Q. (4.13) 

KOq 

Eq. ( |4.13|) is the holographic computation of the vev of the r.h.s. of the trace 
identity ( |4.1| ), to leading order in the expansion in A. Using the expansion f3(\) = 
— 6qA 2 + . . . together with the relation A = kAj, we obtain from ( |4.13| ): 



Pv^t) /iv /,'2\ r;fi i \-2 

4A? 



(TrF 2 ) fi = Q0M 3 N 2 Q. (4.14) 



which exactly matches the l.h.s, eq. ( |4.3| ). This is a nontrivial consistency check of 
our setup. 

Notice that the matching of the conformal anomaly is independent of both k 
and bo, i.e. it is insensitive to the relative normalization of the dilaton field A with 
respect to the true 4D Yang-Mills coupling, A t . 



19 In doing this computation one should actually use a gauge-invariant fluctuation. This can be 
defined as <5$g.z. = (5$ — <&/A6i/j, where ip is the part of the metric fluctuation that couples to 
T^, and is proportional to 8 A. However, close to the UV boundary, Sip ~ r 4 , 5$ ~ r 4 logr and 
$/A ~ (logr) -1 , therefore S&g.I. — » as r — > 0. 
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5. Thermal phase transitions 



In this section we will derive the building blocks necessary to obtain the main theo- 
retical result of this paper: confining backgrounds exhibit a first order phase transi- 
tion, whose features precisely mimic those of the large- N c Yang-Mill deconfinement 
transition. We also show the converse: non-confining theories do not have a phase 
transition at finite T. 

The primary information we will need to extract is the dependence of the temper- 
ature on the horizon position r^, or better on the value of the dilaton at the horizon, 
Xh- Although this is difficult to obtain for generic Xh, it is nevertheless possible to 
determine the asymptotic form of T(Xh) for A^ very small and very large (in a sense 
that will be defined below), corresponding to very large and very small horizon area, 
respectively. 

5.1 Horizon in the UV region 

As we have discussed in Section 3, the black- hole metric approaches asymptoti- 
cally the ArfS'-Schwarzschild solution for small r, with b(r) ~ £/r and f(r) ~ 
1 — (C/4)r 4 /£ 3 , and with a logarighmically running dilaton, A(r) ~ (— 6 logrA) -1 . 

For a large enough value of the constant C, f vanishes at a small enough r 
such that this approximation is valid all the way to the horizon. In this case, f(r) 
vanishes at r\ ~ 4£ 3 C _1 . This approximation gets better for smaller values of Vh (i.e. 
horizon closer to the AdS boundary). Then, the temperature and entropy are given 
approximately by the AdS formulas (see eqs. ( |3.12| ) and ( |3.26| )) : 

T~ — , S = 4ira(— 1 > (5-1) 
nr h \r h J 



with a defined by eq. (|3.18|) . These expressions can be converted into functions of 
Xh, by using the UV asymptotics of A(r), ( p,18|) evaluated at the horizon 20 : 



r h ~ -A~ 6 e ^ (5.2) 

In particular, T(A/ l ) is a decreasing function near the boundary. 

The relation ( |5.1| ) is corrected by the logarithmic running in the UV. Neglecting 
non-perturbative contributions, the scale factor is given by eq. Q2.17] ), from which 
we can obtain the logarithmic corrections to the thermal function f(r) through eq. 
( |3.8|) and (|3.9| ). Through (|3.10| ) we obtain the corrected relation between T and rh 
(see Appendix |D| for the details): 

4 1 



T 



nr h 



+ ... 



(5.3) 



9(logrvV) 2 

20 This is justified in the limit of small (small Xh), since as we have discussed in Section 3, 
the deviations from the zero- temperature solution are 0(r 4 ) for small r, so eq. ( 2.1S| ) is a good 
approximation all the way to the horizon. 
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The results above allow us to compute the temperature dependence of the term 
G(T) in the free energy coming from the gluon condensate. From eq. ( |3.25|) we can 
write 21 : 

l „ r _ 1, 



15ag = F+-TS = - J SdT+-TS, (5.4) 

where S = 47rcx6 3 (r/ l ). For small r^, we can use the UV asymptotic form for the scale 
factor, eq. ( p.!7|) , and the relation ( |5.3|) , to obtain the asymptotic form of Q. The 
calculation is carried out explicitly in Appendix |D], and the result is: 

7r 4 e t 4 

g ~* T^oo. (5.5) 

45 (log^-) 2 

This equation has important consequences. It implies that the gluon condensate 
contribution to the free energy is subleading at large T, by a factor (logT) -2 , with 
respect to the term TS/4 ~ T 4 . Therefore, for very large BH's we can write: 

~ -(ttV 3 )T 4 , T^oo. (5.6) 

On the other hand, at very high temperatures, pure SU(N C ) Yang-Mills theory be- 
haves as a free gas with ~ iV 2 degrees of freedom, and its free energy is approximated 
by the Stefan-Boltzmann formula 22 

^--^iV 2 T 4 , T^oc. (5.7) 

Comparing the last two equations allows us to fix a = V3N 2 / (45tt 2 ) . From the 
definition of a, eq. (|3.18|) , this is equivalent to fixing the 5D Planck scale in AdS 
units, in a model independent way: 

(M/) 3 = -L. (5.8) 
4571^ 

5.2 Horizon in the IR region 

Now, we are going to answer the question: what is the behavior of the function T(rh), 
for large r^l We will find that the answer depends exclusively on whether or not the 
corresponding zero-temperature solution (more specifically, the special solution in 
the classification of Section 2) is confining. In the non-confining case, the black-hole 
temperature decreases to zero for large r^; on the contrary, for confining theories, 
T{rh) asymptotically grows with r^. This distinction has dramatic consequences for 



L It is enough to write the relation below as an indefinite integral, since adding a constant does 



not affect the result for the high temperature asymptotics (5.5). A more precise relation, with 



integration limits, will be given in Section 5.2.3 



22 This is unlike the case of AT = 4 super YM. There the theory in the UV is strongly coupled and 
a non-trivial strong coupling calculation is needed to establish the correct coefficient, Eq]. 
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the thermodynamics of the model: as we will show in the next section, it implies that 
confinement is in one-to-one correspondence with the existence of a phase transition 
at a finite T c . 

As a byproduct of our analysis, we will find that when the horizon is deep in 
the IR ( i.e. the black- hole mass is very small), the black- hole geometry is well 
approximated by the geometry of the zero-temperature special solution (provided it 
exists). More precisely: the special solution is the only zero-T geometry that can be 
lifted to a black-hole of arbitrarily small mass. 

To find the large behavior of the black-hole temperature, we must solve eqs. 
( |3.5| - |3.6|) for r and "close to the singularity" of the zero-temperature solution. 



This question is somewhat ambiguous, since in different zero-T backgrounds the 
singularity can be at a finite or infinite value of the conformal coordinates, and it is 
unclear how exactly to identify the "asymptoti" region. 

A coordinate-independent resolution of this ambiguity consists in solving the 
equations using $ as a radial coordinate, and identifying the asymptotic regions 
according to the value of 't Hooft coupling A = e . Indeed, in the zero-temperature 
background this quantity covers the whole range from in the UV to +00 in the deep 
IR, no matter what is the position of the IR singularity in conformal coordinates. So 
it makes sense to define the black-holes whose horizon is in the deep IR (with respect 
to the zero-temperature background) as those where A attains a large value at the 
horizon. In any given black-hole solution A(r) is a monotonically increasing function, 
and Xh = A(rVi) is the maximum value it can attain. More precisely, we consider 
"large" A as the region A ^> A such that the potential V(A) is well approximated by 
its asymptotic form: 

V(X)~V 00 X 2Q (log\) p , A»A , (5.9) 

for some real P, and Q > 0. The actual value A of course depends on the specific 
form of V(A). Recall that confining theories correspond to Q > 2/3 or Q = 2/3, 
P > 0, and to make sure that there exists a special solution with a good (i.e. 
repulsive) singularity we must further restrict Q < 2a/2/3. The zero-temperature 
superpotential for such solution obeys the asymptotics: 

W (\)~W 00 \ Q (log\) p / 2 , A»A (5.10) 



As shown in Appendix [F.2| , the horizon value A^ uniquely determines the tem- 
perature, once the UV asymptotics are kept fixed. In the next two subsections we 
will determine the behavior of T(\ h ) for X h in the asymptotic region, A^ ^> A . An 
efficient way to attack this problem is to generalize the superpotential technique 
described in Section 2 to the finite-temperature case. This will allow us to give a 
solution of the system, for large Xh, in the entire large-A region. Another method 
that uses diffeomorphism invariant variables is described in Section [7[ 



- 32 - 



5.2.1 Solution in the large-A region 

We will now investigate the features of the black-hole solutions, in the case when the 
horizon is deep in the IR. Consider the situation where the horizon value of the 't 
Hooft coupling A^, is large enough to be deep in the IR asymptotic region, A^ A , 
defined in the previous section. 

In this case it is possible to find analytically a good approximation to the solution 
in the region A <C A < A^, and to match it to the zero-temperature solution in the 
part of the asymptotic region far from the horizon, Ao A This will fix all 

the integration constants, and give the temperature unambiguously as a function of 
X h for X h > A . 

An explicit (approximate) solution for large Xh is found by defining a thermal 
superpotential, i.e. a function W(&), that generalizes the zero-temperature superpo- 
tential W (Q) to the black-hole background, such that the scale factor and dilaton 
equations (in domain wall coordinates) reduce to the form: 

4 dW 
A'(u) = --W, <&'(«) = (5.11) 

The details of this formalism are developed in Appendix [F|. As shown there, the 
superpotential W / ($) and the logarithm g of the thermal function / appearing in the 
metric satisfy a pair of coupled differential equations in the variable $, which take 
the form: 



d*g + ?&)a*W+%W = ^W, g = logf (5.12) 
o$gJ 9 

4 4 n 64 

--W(d$W)(d$g) - - (d^Wf + — W 2 = e- 9 V($). (5.13) 

In these variables the system of Einstein's equations splits into two separate 
decoupled systems, and solving eq. ( p,12| - |5.13] ) determines the metric and dilaton 
through eqs. ( 5 . 11| ) solely from the knowledge of W(<&), as in the zero-temperature 



case. Here, though, the superpotential cannot be taken as an independent function, 
but it depends on temperature through the coupling to g($) in eq. ( |5.13|) 

As shown in Appendix |F.3.2j , the solution of the system (|5.12 - [5T3|) in the whole 



region where ( |5.9|) is valid has the form: 

/ ($) ~ 1 - C^j e -^(*fc-*), W($) ~ A Q (log A) p/2 , $ < $ < $ 



hi 

(5.14) 



where the constants K and R are given by: 

9Q W ) ' 2 V* ' 9Q 
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More precisely, the approximate equality signs in eq. Q5.14 ) stand for dropping terms 
of 0(1/$) in W and log[l - /]. 

Notice that the asymptotic solution for the superpotential has automatically 
the same asymptotics of the special solution W Q (A), eq. Q5.1Q ): the presence of a 
regular horizon at large A^ has lifted the degeneracy present in the zero-temperature 
superpotential equation, selecting a single solution, which happens to be the unique 
one with the good singularity! We will come back to this point later in this Section, 
and in Section 7. 

The approximation W ~ W D is valid in the whole asymptotic region. Thus, the 
metric and dilaton have the same asymptotics as in the T = case, all the way to 
the horizon. Moreover, /(<&) is very close to unity already in the asymptotic region, 
specifically where $o <^ $ <^ $h- Thus, for all smaller values of $, including those 
outside the asymptotic region, /($) — 1 and the superpotential equation reduces to 
the zero-temperature one, with special solution W ($). This implies that W ~ W Q 
is a good approximation in the whole range of $, not only asymptotically. The 



equations ( |5.11| ) for the metric and dilaton therefore reduce to the zero-temperature 
equations, whose solutions are classified by the value of the integration constant A. 
This implies that we can fix the remaining integration constant of the system (the 
scale A) to the same value as in the T = background when we integrate (|5.11| ). 
With this choice, all the integration constants in the asymptotic solution are fixed. 
Thus, for large A^ the complete thermal solution for the scale factor and dilaton is 
everywhere well approximated by the zero-temperature form: 



A{r) ~ A {r), 



$(r) ~ $ (r). 



(5.16) 



From eq. ( 5.14| ), expressed in r-coordinates (by setting f(r) = /($(r))), we observe 
that, at any fixed r, taking the limit $>h — > +oo brings the metric to match the 
vacuum form with f(r) = 1. Therefore the limit A^ — > oo corresponds to the (point- 
wise) limit of the black-hole going to the zero-temperature metric. 

In the particular case of power-law asymptotic and singularity at r = oo, corre- 
sponding to Q = 2/3, P = (a — l)/a (see j28[), the solution for large r is : 

r 



A(r) 
$(r) 

m 



3 
2 

1 - 



l-a 



f r 

l)log( z 



exp 



(5.17) 
(5.18) 

(5.19) 



where L is the same length scale appearing in the zero-temperature solutions ( 2.2Q| ). 
The relation between temperature and horizon position for these black-holes is: 



T 



\f(r h )\ 

47T 



a-1 



1 r 

it, 

4vr L a 



1 + 



(5.20) 
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Notice that, for a > 1, i.e. when the theory is confining, the temperature increases 
as a function of rv In the next section we will see that this behavior is characteristic 
of general confining asymptotics. 



A no-go theorem 

As discussed in detail in Appendix [E|, the asymptotic behavior (|5.9|) does not fix 
uniquely the asymptotics of the zero-temperature solution, but rather allows various 
possibilities summarized in Appendix |E.4| . An important conclusion that can be 



derived from the discussion above is that black-holes which probe the asymptotic 
large-A region necessarily match the "special" zero-temperature solution with asymp- 
totic ( |2.28|) , as Ah — > oo. All other solutions cannot be lifted to black-holes with 
arbitrarily small mass (i.e. whose metric is arbitrarily close to the zero-temperature 
metric) . For an alternative derivation of the same results, see Section [T3[ We thus 
have the following: 

No-go theorem: The only vacuum solutions of Einstein- dilaton gravity, with a po- 
tential satisfying ( \5. 9j ), that can be continuously lifted to a regular black-hole with 
arbitrarily small mass, are the ones stemming from a superpotential with the "spe- 
cial" asymptotics ( \2.28j) . 



In particular, this implies that the two continuous families of zero-T solutions 
with "generic" and "bouncing" superpotentials (all of which have bad singularities) 
cannot be promoted to black-holes with an arbitrarily small mass (i.e. such that the 
metric looks like that of those zero-T solutions almost everywhere). 

The zero-T superpotentials with bad singularities were already discarded as 
pathological in ||28|| , since they exhibit singularities which are not screened from 



physical fluctuations, and do not have a well defined eigenvalue problem. If we sub- 



scribe to the criteria laid out in 2S , the no-go theorem above gives another reason 



to discard these kinds of solutions: in [29] a singular background was considered 



acceptable if the singularity can be "hidden" behind a regular black-hole horizon by 
an infinitesimally small deformation of the metric 23 . 



23 The physical reason behind this criterion was that, if the singular background has a holographic 
interpretation as the zero-temperature vacuum of a 4D theory, then it should be possible to obtain it 
continuously as the T — ► limit of a thermal state, described holographically by a black-hole. This 
reasoning assumes that a small black-hole, which is infinitesimally close to the vacuum background, 
has also an infinitesimally small temperature. In the following section we will see that this is not 
necessarily true: there are cases (e.g. in the presence of a deconfining phase transition at finite T) 
when the M — > limit of the black-hole corresponds to a high temperature limit. Thus, in the most 
interesting cases of confining backgrounds, the criterion of p9| loses part of its physical motivation. 
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5.2.2 Big black-holes and small black-holes 

Now, let us determine the temperature of the asymptotic solution obtained, for 
general P and Q in (|5.9|). It is easiest to do so as a function of We have (see eq. 



4vrT 



df(r) 



dr 



dm 



<i$ du 
du dr 



dm) 



dW 



(5-21) 



The quantities df/d& and dW/d& are obtained from eq. ( |5.14 ). On the other hand, 
A(<$>) can be obtained by combining the two eqs. ( |5 . 1 1|) , 



dA 



9d\ogW 



which upon integration which gives: 



4 2P 

A($) ~ $ H log$. 

1 ; 9Q 9Q S 

Inserting these expressions in eq. ( |5.21| ) we obtain: 



T(® f 



K Q (gj^/£) ^ P/ 2 +(2P/9Q) 

4tt h 



*h > $0 



(5.22) 



(5.23) 



(5.24) 



From equation ( 5.24 ) we deduce that the thermodynamic behavior of black- hole 
solutions is very different according to the IR confining properties of the zero-T 
solution: 

1. if Q > 2/3, or Q = 2/3 and P > (i.e. the T = theory confines): 

the temperature increases with so very small black-holes (large A/J corre- 
sponds to high temperature solutions. 

2. if Q < 2/3, or Q = 2/3 and P < (the T = theory does not confine): 

The temperature decreases with and the limit $^^00 corresponds to 
T -»• 0. 



Since, in all cases, T decreases with in the UV regime, it follows that: 



In confining theories, T($^) must have a minimum value T min below which there is 
no black-hole solution. In this region the only remaining solution is the thermal gas. 
In particular, T m i n cannot be zero: this would imply the existence of a point where 
f = 0, but from eq. ^3. we observe that if f vanishes at some point where b is 
regular, it must vanish everywhere. 

What we have shown implies that confining theories admit (at least) two branches 
of black-holes with the same temperature, starting at T = T min . For T < T min there 
are no black-hole solutions, whereas for high enough temperature there are exactly 
two black-holes: 
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Figure 1: Typical plots of the black-hole temperature (a) and free energy (b) as a function 
of the horizon position r^, in a confining background. The temperature features a minimum 
at r min ; that separates the large black- hole from the small black-hole branches. 

1. the big black-hole has its horizon closer to the AdS boundary, satisfying the 
relation (|5TT|) , 

r* ~ -L . (5.25) 

2. the small black-hole has its horizon deep in the interior. With the asymptotics 
( |5.19| ), the horizon position is related to the temperature by (inverting ( |5.20| ): 

r small ^ L l TL \ ; a>1 ( 5>26 ) 

The typical situation (with two branches of solutions) is depicted in figure |l|. 
Generically, the big black-holes are thermodynamically stable whereas the small 
black- holes are unstable. Using c v = TdS/dT and (|3.26|) , we obtain, 

dT 3ST dA 

GA/j Cy CLAh 
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If dA(Xh)/d\h < 0, the condition for thermodynamic stability i.e. cy > coincides 
with dT/dXh < 0, which is true for big black holes; on the other hand for the 
small black-holes dT/dXh > 0, hence if A(Xh) is monotonic, they have negative 
specific heat and are unstable. Although typically dA(Xh)/dXh < is satisfied, since 
generally increasing Aft means going deeper in the IR, it is hard to tell whether this 
is true for all values of Aft: recall that there is a "hidden" dependence on Aft in 
A(X), due to the fact that the form of the scale factor is temperature-dependent. 
The condition dA(Xh)/dXh < is certainly obeyed asymptotically both for small 
Xh and for large Aft, where all black-hole metrics reduce to the zero temperature 
solution: for small Aft, A(Xh) ~ (feoAft) -1 , and for large Aft, A(Xh) ~ —2/3 log Aft. 
However one might have some range of Aft where dA(Xh)/dXh > 0, thus the small 
black-holes are stable, or the big ones unstable. We have found some numerical 
evidence of this behavior, where the violation of monotonicity of A(Xh) is always in 
a very limited range of Aft, and only in the small black- hole branch. Since the big 
black-holes are the thermodynamically dominant solutions above T c , we don't expect 
them to become unstable for T > T c , since the dilaton potentials we use can always 
be written in terms of a superpotential: this guarantees the validity of the Positive 



Energy Theorem ||39|| , which in turn guarantees stability of the vacuum. 

We now illustrate the general behavior just described with a few simple concrete 
examples, for which we solved numerically Einstein's equations, starting from an 
explicit form of the dilaton potential. In the case of the explicit potential that was 
considered in p8| , the numerical results for the thermodynamics were discussed in 



|30H , and they indeed agree with the general results discussed above. Here we present 
the result for a simple class of potentials, which have the IR asymptotics of the form 
( |5.9p for general P and Q: 

^ = 1 + A + V X X^ [log (1 + V 2 X 2 )f (5.28) 

where we have set the AdS scale to unity. Also, we have set the coefficient of the 
linear term to one. This can be done by rescaling A and redefining the coefficients 

The numerical results for the function T(Aft) for Q = 2/3, fixed Vi, and various 
values of P are shown in figure |2|. They confirm our general analysis: confining 
potentials (P > 0) feature a minimum black-hole temperature, whereas in non- 
confining ones (P < 0) the temperature can be arbitrarily low. 

As discussed earlier, the temperature is guaranteed to be a single- valued function 
of Aft, but the same is not necessarily true if we consider T as a function of the horizon 
position rft. A rather striking example of this fact is shown in figure [| it displays 
the plot of the curve T(rft) in a case (Q = 2/3, P = 2) where the singularity of the 
vacuum solution is at a finite value of the conformal coordinate. For comparison, 
the curve T(Aft) in the same model, shown in figure ^| (b), is single- valued. One 
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Figure 2: The temperature as a function of the horizon value of A in the model specified 
by the potential ( |5.28 ), for Q = 2/3 and various values of P. The other coefficients are 
fixed to V\ = 10, V*j = 100. The confining models (P > 0) feature a minimum temperature 
at finite A^; in the non-confining model (P = —1) the T(Xf l ) monotonically decreases to 
zero; In the borderline case (P = 0) T^X^) decreases monotonically to a finite value as 
A — > oo. 



unexpected feature is that in this case there are black-holes whose horizon is well 
beyond the position of the zero-temperature singularity at r = r$. Nevertheless, as 
the temperature increases along the small black-hole branch, we have — > r$ as 
expected. 
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Figure 3: Temperature as a function of (a) in the model ( 5.28| ) for Q = 2/3 and P = 2. 



The temperature diverges for A^ — > oo, for which — > ro. It is a single- valued function of 
A/j, but not of r/j. 

5.2.3 An integral representation for the free energy 

As described above, in many cases there are more than one black-hole solutions, 
some of them having negative specific heat. This is reflected in the fact that A/j 
(or r/j) as a function of T is generically multi-valued. The formula for the free 
energy that we derive here will encompass all of the different branches under one 
integral representation. This form will be used to prove the proposition in the next 
subsection. It is also a very convenient form for numerical evaluation of the free 
energy, in the process of determining the thermodynamics of the system numerically. 
Here we present the discussion for the case of two branches of solutions, one small 
and one big black-hole for simplicity of the presentation. However the final result is 
easily generalized to multi-black-hole cases. 

Let us denote the free energies of the small and the big black-hole by Fs and Tb 
respectively. 

Integrating the first law, S = —dTjdT for the big BH, one obtains, 

Tb{T) = T min - [ S B dT, T > T min (5.29) 



J T 

where T min = J"(T min ). 

In order to determine the integration constant in (|5.29|) one can make use of the 
same formula but on the small BH branch. This is clearly suggested from fig. § 
where we depict a generic form for the function FiT) in case of two-branches. First, 
note that T m i n is the same on both branches (see figure 

•F B^miri) = 3~ S^min) = -F min- (5.30) 
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Thus one has: 



F S {T) = Fmin - / S s dT, T > T min (5.31) 

The small black-hole free energy vanishes in the limit of zero black-hole size (where 
the metric coincides with the zero-temperature background), which for confining 
backgrounds is the T — > oo limit. This allows to write T m in as: 

3~min = — Ss dT, (5.32) 




Figure 4: Black hole free energy 



Combining (|5.29|) with (|5.32 ) one obtains an integral representation for the free 



energy that only depends on the area of the horizon, and is valid on both branches. 
It can be put in a simpler form in the A^ variable. Using ( |3.26|) one obtains, 



f Xh ~ dT ~ 

Tb(K) = -4vrcr / b 3 (\ h )-~-d\ h , \ h < \ miw (5.33) 
J co d\h 

where \ m i n is the horizon position of the minimum temperature black-hole. Note 
that the two branches are combined in the integral as both b and T are single- valued 



as functions of A^, but are not as functions of T. In appendix [H.10| we present 
further useful formulas regarding this integral representation. We remark that RHS 
of (|5.33|) is finite everywhere except at A^ = where Tb —oo. Finiteness near 
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the singularity A^ = oo follows from the fact that, for all of the confining cases, b 3 
vanishes exponentially faster than dT(Xh)/dXh- 

We note the remarkable fact that the free energy is completely determined by 
the knowledge of area of (the small and the big) black-hole horizons. This means 
that the entire thermodynamic properties of the dual field theory is encoded in the 
horizon function of T. We stress that the area of the big BH horizon 

only is not sufficient; this misses the integration constant ( 5.30| ). It is therefore not 



sufficient to determine T c for instance. An alternative way to say this is as follows: as 
we showed in eq. ( |3.25| ) there are two contributions to the free energy: the entropy 
and the condensate. Here we learn that, we need both branches to disentangle the 
two contributions. 

5.3 Confinement and phase transitions 

In this section we show the one-to-one connection between color confinement (in the 
vacuum background) and the presence of a deconfining transition: 

Proposition: 

i. There exists a confinement- deconfinement phase transition at finite T, if and 
only if the zero-T theory confines. 

ii. This transition is of the first order for all of the confining geometries, with a 
single exception: 

iii. In the limit confining geometry A (r) —>■ —r/L (as r — > oo), the phase transi- 
tion is of the second order and happens at T = 3/(47rL). 

iv. All of the non-confining geometries at zero-T are always in the black-hole phase 
at finite T. They exhibit a second order phase transition at T = + . 

We outline our demonstration in the coordinate system where A is chosen as the 
radial variable. Being diffeomorphism invariant, our arguments apply to all of the 
confining zero-T geometries that are described in Section ^p 4 . 

We first consider the geometries that confine the color charge at zero T. In 
Section ( |5.2j ) we have shown that there exists an extremum of T(Xh) if and only if 
the zero T theory confines. Thus, for confining potentials there exists such points 
X m i n that satisfy, 



dT 



dX 



= where T(X min ) = T min > 0. (5.34) 



One can state the part iii of the theorem in a diffeomorphism-invariant way by stating that the 
borderline geometry is defined as lim;^oo(^ + \) log A = pq |. 
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Figure 5: Temperature as a function of A^ for the infinite r geometries of the type A — > r a . 
Black holes exist only above T m i n whose precise value depend on the particular zero-T 
geometry. 

In figure [5| a schematic plot of T(Xh) for the typical examples of confining, non- 
confining and the borderline cases is shown. 

Now, it follows from the positivity of the entropy and the first law of thermo- 
dynamics S = —dJ^/dT, that the extrema of T(Xh) coincide with the extrema of 
To see this we write — S = dJ- ' jdT = (dJ-/dXh)/(dT/dXh) and observe that 
in order for S > at all A^, the extrema of T(Xh) should coincide with the extrema 
of J-{Xh). Hence, for the confining potentials there exist at least one extremum of 
J-'(Xh). See figure || for a schematic plot of T(Xh). 

For simplicity, here we shall assume that there exists a single A m j n in the confining 
geometries, and no extremum in the non-confining ones. We shall comment on 
the multi-extrema cases in the next subsection and carry out a detailed analysis in 
Appendix 0. 

Let us first reproduce here the integral representation derived in the previous 
subsection: 



The integrand is positive definite for A^ > X min and negative definite for A^ < A 
Therefore, evaluating it at X min , one has, 




(5.35) 




(5.36) 
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Figure 6: Free energy density as a function of Xh for the infinite r geometries of the type 
A -> r a . 

On the other hand, when one evaluates (|5.35|) on the boundary, one finds 

lim T(\ h ) = -oo. (5.37) 



This follows from the UV asymptotics described in Section [O. Therefore there must 
exist a point A c where it vanishes: 

^(A c ) = 0, (5.38) 

see figure fig, |6|. This proves that there exists a phase transition if the zero-T theory 
confines. The transition temperature T c is greater than T min because A c < \ m in an d 
dT/dXh < at A c . This result confirms our intuitive picture that, as the temperature 
is increased, first the small and big BHs form at a temperature T min , where the 
minimum energy configuration is still the thermal gas, and as T is kept increasing, 
the big black-hole takes over the thermal gas phase at a higher temperature T c . The 
true free energy of the system, i.e. the function F(T) evaluated on the minimum 
energy configuration, is shown schematically in fig. [7]. 

On the other hand, for the non-confining theories (see the case a < 1 in figs. 
H and U), T(\h) is always monotonically decreasing. This is because, in the UV 
T oc exp[l/&oA/j]. Hence it is monotonically decreasing in the UV and there exists 
no extremum of this function. Therefore its derivative cannot change sign. From 
( |5.35|) , it follows that J-"(\h) < for all A^ and that there is no phase transition for 



the non-confining potentials at any finite-T. 
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Figure 7: A schematic plot of the free energy in the typical case of a confining geometry. 
"B" and "TG" denote the big black-hole and the thermal gas respectively. At T c there is a 
first order transition. This plot corresponds to evaluating F(T) of fig. ||| on the minimum 
energy configuration. 

We have so far proven part i. of the proposition. Parts ii. and iii. are proven as 
follows. The order of the phase transition is determined by the latent heat: 

L h = E(X C ) = S(A C )T(A C ). (5.39) 



It follows from ( |3.26| ) that, S(X C ) is non-zero unless A c coincides with the singularity. 



Therefore Lh > and the phase transition is first order for the standard confining 
geometries. For example, this is the case for a > 1 in infinite geometries, (see figs. |5] 
and|). 

On the other hand, in the borderline case a = 1, T c = T min and, at A c the entropy 
vanishes, because A c coincides with the singularity and b(X c ) — > there. Thus, in this 
case the transition is of second order. It would be interesting to find other examples 
of second order phase transitions in the Einstein-dilaton system. 

The last part of the proposition follows simply from the fact that, for all the non- 
confining geometries, J-{\h) < for all A^ and it vanishes only at the singularity. 
At this point both T and S vanish as the area of the horizon vanishes. Thus there 
is a trivial second order phase transition at T = + and the system is always in the 
black-hole phase for any fmite-T. 
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5.3.1 Geometries with multiple extrema 

We demonstrated our proposition under a single assumption, that the function T(A/ l ) 
has at most one single local extremum. We did not find any counter-example to this 
assumption in our numerical studies, and it is a logical possibility that, with the given 
assumptions for V(A) (that it is a positive and monotonically increasing function of 
A that limits to a constant at A = and diverges exponentially as A — > oo), multiple 
extrema cases never occur. However, we can not rule out these possibilities by 
analytic arguments, therefore they should be considered in order to complete our 
demonstration. Moreover, as we discuss below, they bear interesting possibilities for 
new types of phase transitions. 

In general, there may exist theories which admit more than one small and one 
big black-hole. In these cases, the functions J-{Xu) and T(Xh) are complicated and 
admit many extrema. As a result F(T), evaluated on the entire set of solutions (not 
only the lowest energy solution) may be a complicated multi-valued function with 
many cusps and crossings, see fig. [T7| in Appendix G for an example. 

The proof of part i. of the proposition extends without changes to the general 
case with multiple extrema. To prove points ii. iii. and iv. in full generality, however, 
we must make an additional assumption on the behavior of the entropy as a function 
of Xh- we must assume that a generic black-hole in a given branch has larger entropy 
than a generic black-hole in the next (with larger Xh) branch. This is a weak version 
of monotonicity of S(Xh), and it is a sufficient condition for the full proposition to 
hold, although it might not be necessary 25 . 

The details of the general case are presented in Appendix |G]. The upshot of 
the analysis in |GJ is that, regardless how complicated the system is, F(T) evalu- 
ated on the minimum energy configuration always have a similar form to fig. [?| as 
in the single extremum case. Since, in the infinite volume limit, only the lowest 
energy configuration is relevant for the thermodynamics, our demonstration in the 
previous subsection directly carries over to cases of multiple extrema. The precise 
statement is that regardless how complicated the function T(Xh) is, there exists a 
confinement- deconfinement phase transition if and only if the corresponding zero-T 
theory confines. 

However, our analysis in Appendix |Gj shows another interesting possibility: in the 
multiple extrema cases, there may exist first order phase transitions between different 
deconfined vacua. The temperature of these transitions are always higher than the 
deconfinement temperature T c . In general, there may be an arbitrary number k, of 
such transitions with 7*. > Tk-i > ■ ■ - T\> T c . In the dual geometric picture, these 
transitions occur between different big black-hole geometries. For a number k of such 
transitions, the function T(A^) should possess k local minima and k — 1 local maxima. 

25 Strict monotonicity is too strong a condition, since we have found numerical counterexamples 
where this is not satisfied. 
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Figure 8: The function F(T) for a multiple extrema case, evaluated on the minimum 
energy configuration. Here k = 1, i.e. there are two big-black- hole geometries denoted by 
B\ and B<i with a first-order transition between them at T^. TG denotes the thermal gas 
geometry, that takes over below T c . 



This means that, there should be k different pairs of small and big black- holes. This 
is a necessary condition but it is not sufficient. The sufficient condition follows from 
the particular shape of F(T) that leads to such transitions between different big 
black-hole branches. An example is discussed in Appendix We present the free 
energy for the minimum energy configuration, for the case k — 1 in fig. [S]. The fact 
the these transitions are always first order follows from discontinuity in F'(T) at TJ, 
as shown in fig. |] and Appendix |G|. 

It was argued that there could be a series of phase transitions in large- N c gauge 
theory that would correspond to a partial, step-by-step breaking of the center 
of the gauge group, [Q. At large N c there is room for an arbitrary number of 
such steps. The order parameter corresponding to the 1-th such transition would be 
(Tr V), namely the 1-th power of the Polyakov loop. It is plausible that such phase 
transitions may be in the same universality class as the ones described above. 

Let us remark however, that, neither in the lattice studies of nor in our 
numerical investigations, one has encountered such transitions. They may exist as 
an exotic possibility in our set-up. 
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5.4 Similarities with the Yang Mills deconfinement transition 

We have found that backgrounds which exhibit confinement, also exhibit a decon- 
finement phase transition at some finite temperature T c , above which the black- hole 
phase dominates. The qualitative features of the phase transition and of the ther- 
modynamics of the deconfined phase are remarkably similar to those found in four 
dimensional pure Yang Mills theory at large iV c . 

Below we list some of the mo del- independent features of the gravity phase tran- 
sition that match the gauge theory side. We will analyze the quantitative agreement 
in concrete models in a separate work [36]. 



It is a confinement-deconfinement transition. In particular in the high temper- 
ature phase the confining string tension vanishes. 

The Polyakov loop is the order parameter for the confinement-deconfinement 
transition of SU(N C ) YM theory. The vev of the Polyakov loop (V) is zero in 
the confined phase and it acquires a non-zero expectation value above T c . Here 
we see the analogous behavior in the dual geometric picture [|]: Holographically, 
the Polyakov loop is described by a classical string embedding that wraps the 
Euclidean time direction. In the thermal gas solution the time circle is non- 
shrinkable, hence the action of the string is infinite, giving a zero vev for the 
Polyakov loop. On the other hand, in the black-hole solution, the time circle 
shrinks to zero size at the horizon and one obtains a non-trivial vev. 

Both in large- N c Yang-Mills and in the gravity theory the phase transition is 
first order, with a latent heat that scales as On the gravity side the latent 
heat per unit volume is given by (see eq. fl3.27p ) 

L h = N 2 c Ml (l5£(T c ) + 3vrT c 6 3 (r c )) = Y c 2 ^r 3 £(T c ) = N 2 c ^r 3 T c b 3 (T c ). 

(5.40) 

As in large- N c Yang- Mills, the quantity F/N 2 serves as an order parameter: it 
is of order one in the deconfined phase, and zero in the confined phase 

The high temperature behavior is the same of a free gluon gas, up to logarithmic 
corrections. This behavior can be seen in the temperature dependence of the 
pressure p = — JF/V3, energy density p, and entropy density s : from eq. ( |5.G| ) 
and standard thermodynamic relations, and with the choice (|5.8| ) for the Planck 
scale, we have: 



3p(T) 

rp4 



3 s(T) 
4 T3 



7T 



2 



]5 K T^oo. (5.41) 
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In Figure [| we present an explicit example of the behavior of these quantities 




Figure 9: The energy density p, entropy density s and pressure p in the black- hole phase 
for T C <T < 5T C . 

in the deconfined phase, up to a temperature of 5T C , derived from our model 
with a potential of the form flET2BD with Q = 2/3 and P = 1/2. 

The deviation from conformality is expressed by the trace of the thermal stress 
tensor: 

' U/T* T^T C 

(5.42) 

: ^(logT/T c )- 2 T^oc, 



rp4 



135 

as can be derived by combining eqs. (|4.3|) and (p.5|). A concrete example of 



the temperature dependence of this quantity in the black-hole phase is shown 
in Figure (|10). 



The speed of sound is given by cj. = (dp/dp)s = S/C v , where C v = dE/dT 
is the specific heat. As expected in the deconfined phase of pure Yang-Mills 
in 4D, this quantity approaches from below the conformal value, c 2 s —>■ 1/3 
as T — ► oo. This can be seen using the high-temperature expansion of the 
thermodynamic quantities derived in Appendix |. A concrete example in our 



setup is shown in Figure 11. 



6. The axion background at finite temperature. 

The effect of a non-trivial vacuum angle in Yang- Mills is captured by including in the 
bulk a five- dimensional axion. The axion a is dual to the instanton density Tr[FAF]. 
In particular its UV boundary value is the UV value of the QCD #-angle. Moreover, 
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Figure 10: Temperature dependence of the trace of the stress tensor in the black-hole 
phase for T c < T < 5T C . 
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Figure 11: The speed of sound in the black-hole phase for T c < T < 5T C . 



its profile a(r) in the vacuum solution may be interpreted as the "running" #-angle 
in analogy with the dilaton, that we interpret as the running coupling constant. This 
was explained and justified in . 

The axion action is suppressed by 0(1/ N^) with respect to the action for the 
other fields, (O) : 



Mi 3 



d 5 xy^ Z(\) (d^af 



(6.1) 
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where Z(X) has the following asymptotic expansions, [28 



lim Z{\) = Z a [1 + 0(A)] , lim Z{\) ~ A 4 (6.2) 

A^O A^oo 

The scale Z a determines the topological susceptibility, while the strong coupling 
asymptotics are dictated from glueball universality. 

At zero temperature the axion solution that is compatible with known properties 
of large-iV c YM is 

f r ° d T 

a{r) = (6 UV + 2nk ) ^ e ^ (A) , k e Z (6.3) 

JO e 3A Z(\) 

where 9uv is the UV value of the 6 parameter defined as an angle in the range [0, 27r), 
ko is an integer that labels oblique confining vacua and is determined by minimizing 
the ^-dependent vacuum energy 

Ml (6 uv + 2nkY 
E(6 UV , k) = — — £ _ 3A (6.4) 

Jo dr z{\) 

Finally r is the position of the singularity of the zero-temperature geometry in the 
radial coordinates. 

From this solution we can extract the topological susceptibility and topological 
density condensate as 

v = Z{0)M * (rr[FAF l ) = -g!E! (^ + M °) ( 65 ) 

^ rrp drZ(0) ' \ A ' L T /N r if £3 rrp , Z(0) 

Jo e* A Z(X) Jo ar e^ A Z(\) 

At finite temperature, below the deconfining phase transition the situation is 
similar to T = 0, since it is the same vacuum solution describing the physics. Above 
T c however, once should switch to the black-hole solution. In the black-hole solution 
the axion background satisfies, 

a+ ( 3A + t + (d x log Z)))j d = (6.6) 
an equation that can be integrated once to 

( " : (6-T) 



/ e 3A Z{\) 

where C a is a constant. Integrating once more we obtain 



dr' 

aif) = 9 UV + C a I - e3A z(A) (6. 



o 



Unlike the zero temperature solution, the non-trivial solution here has a singularity 
at the horizon. Indeed, f(r) ~ fo(r — r^) while both Z(X) and e A are regular 
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there. Therefore the function J Q r j e3 ^' z ^ diverges logarithmically as r — > r^. Since 
the background solution must be regular everywhere, we must necessarily impose 
C a = 0. 

Therefore in the deconfined phase, the axion background is constant 

a(r) = 9 UV (6.9) 

and the topological density condensate (proportional to C a ) vanishes 

(F A -F) dcconfincd = (6.10) 

In fact this can be generalized to higher derivative terms containing the axion. This 
shows that in the deconfined phase, at large N c all moments of the topological density 
vanish. This is in agreement with general expectations at large N c , since in the 
deconfined phase such moments obtain contributions only from instantons, that give 
vanishing e~ Nc contributions as N c — > oo. This expectation is in accordance with 



lattice calculations , ET1 . 



7. Reduction and solution of the system in scalar variables 

Einstein's equations are hard to solve for a generic dilaton potential. This is the case 
even in numerical evaluation. Here we shall present a method to reduce the degree 
of the system of equations from 5 to 2 by introducing variables that are explicitly 
invariant under general coordinate transformations. 

There are a number of nice features of this reduction. First, the thermodynamics 
will only depend on the reduced system. Therefore all thermodynamic observables 
are determined by solving a coupled system of first order equations. This also clarifies 



the number of physical parameters in the theory, see Section |7.6| . It also allows us 
to find analytic solutions (Appendix fJJ). 

Aa another bonus, the UV expansion of the finite temperature metric and dilaton 
used throughout the paper, eqs. (|3.13| - |3.14D , are most easily derived with this 



formalism. This is done in Subsection 7.2 



Finally, this form of Einstein's equation is used in Appendix [H.4| to show the 



existence of black-hole solutions with arbitrary and AdS UV asymptotics. 

This section can be read independently of the rest of the paper. All of the 
ingredients necessary to derive the thermodynamics from the Einstein's equations in 
a diffeomorphism-invariant manner are presented in this section in a self-contained 
way. 

7.1 Scalar variables 



The basic idea was introduced in [2£] for the zero-T case and reviewed in Section 



. Here we present a generalization of [28] to the black-hole ansatz. We propose to 
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solve the Einstein's equations by introducing the following scalar variables: 

v ; 3 A'' y ! AA' K ' 
Note X and Y are invariant under radial coordinate transformations. These variables 
obey the following first order equations: 

dX - 4 fl X* + Y)(l + 3 dhgV ^ (7 2) 

~m -~3 (1 ~ x +Y) { 1 + 8x^r)> (7 - 2) 

r/y 4 y 

— = -_(l-X 2 + y)-. (7.3) 

d$ 3 y J X K ' 



As shown in Section [7.4|, the thermodynamics of the dual field theory are completely 



determined by knowledge of X and Y as a function of $. Roughly speaking, Y is 
dual to the entropy and X to the energy of the gluon fluid. 

It is crucial that the system above always admits a special solution, Y — 0. 
This corresponds to the thermal gas solution (or the zero-T solution for non-compact 
time), that is present for any V . This is because, for the solution Y = 0, the equation 
( |7.2| ) reduces to the corresponding zero-T equation ( |2.10| ) whose only solution with 
fixed IR asymptotics is the thermal gas. 

The metric functions are given in terms of X and Y by integrating eqs. ( |7.1| ). 
Let us introduce a cut-off $ , that plays the role of the regularized UV boundary. 
We call the value of the scale factor A at this point as Aq. On the other hand, the 
black- hole asymptotics requires that g vanishes (/ — ► 1) on the boundary. With 
these initial conditions, integration of ( |7.1|) gives 

A = A + £ ^d<£, (7.4) 

*4y ~ 

(7-5) 



We note that one does not need the cut-off $o in solving (|7.2j ), (|7.3| ). As we prove in 
the sequel, the physical observables only depend on the functions X and Y . Therefore 
they will be independent of the cut-off $o- 

In appendix [H.l| we prove that the reduced X-Y system solves the full equations 
of motion ( |A.17[ - [A]20D in the original w-variable. There, we also provide formulas for 
the derivatives of the metric functions A 1 , g' and $' in eqs. ( |H.1|) , (|H.2|) and (|H.3|) , 
hence completing the full five-degree system of equations. 



Let us finally note that, one can invert the equations (|7.2p and ( |7.3|) for the 
dilaton potential: 

V(Q) = Vo(l + y-X 2 )e~£-£«*( x "£). (7.6) 

This equation will be useful later on. We note that for Y = it reduces to the 
corresponding zero-T equation in 



- 53 - 



7.2 UV Asymptotics 

In what follows we prefer working in the A = exp($) variable instead of $ 26 . The 
black-hole deformation is obtained by turning on Y near the boundary. In other 
words Y should vanish as one approaches the boundary. In addition, the condition 
that the BH solution approaches the thermal gas solution requires that X —* Xq on 
the boundary. We recall that the UV asymptotics of the function X is presented in 
eq. ( |2.12D . Let us now write, 

X(X) = X (\) + 6X(\), A^O, (7.7) 

where 5X <C X for small A. Studying the small A asymptotics of the explicit solution 
for Y, given in eqs. ( [H.9| , [H.10D , one learns that Y vanishes non-perturbatively in A; to 
be precise Y ~ e -4/fcoA^-4b This can a } so j-^ seen assuming that Y is exponentially 
small and than solving ( |7.3| ) in the vicinity of A = 0. Then, it follows from ( |7.2[ ) that 
5X also vanishes with the same exponential factor. 

One derives the asymptotic behavior of the functions Y and SX by solving { \7.'Jj ) 
and ([7.2Q near A = 0. We spare the details of this calculation to appendix |H.5| . The 
result is 

y(A) = Y o e-^(& Ar 4fe , (7.8) 



SX(\) 



Y /2 -C 



^(6 A) (7.9) 



X 

Here Y and Co are integration constants. They retain finite values as the cut- 
off is removed by sending Ao — * 0. These values can be computed by matching the 
solutions above to the full solution of ( |7.2|) and (|7.3|) . Generally these integration 
constants are non-trivial functions of temperature and this dependence is determined 



by the regularity condition at the horizon. This is explained in Section \l.3\ below. 

The physical meaning of these integration constants will become clear below: Co 
determines the energy, Yq determines the entropy and the combination, 

Go = Co - (7.10) 

determines the vev of the gluon condensate in the gluon plasma. 

Finally, we would like to know the UV expansions of the metric functions A, f 
and $ in the radial variable u or r. These can be determined using the asymptotic 
expressions for X and Y given by ( |7.8| ) and ( |7.9[ ) in the formulae for the derivatives 
A', f and $' given by fn| , |H.2|JH.ij|) . 

To make use of these equations, we first define the fluctuations in A and $ in 
the domain wall frame as, 

8A(u) =A(u) - A {u), 5${u) = $(u)-$ (u). (7.11) 



26 For notational simplicity, we shall allow for an abuse of notation by referring to, in fact, different 
functions when we write e.g. V(A) and V(&), related by V(x) — > V(e x ). 



- 54 - 



One then finds 
1 



5A 
f 



G e b o*(b \) 



-ib 



1 



G (Af) 4 e 4n/ W£)- iti/9 + 



-16/9 



l-Y o e^(& A)- 4b + 



aa u /i. 



9 _^t_ 

5$ = -G e fc oA(fe A) 



-46-1 



+ 



l-Y (A£) 4 e 
^G (A£) 4 e Au/e ( 



-nil) 



-16/9 



+ ■ ■ 



(7.12) 
(7.13) 



-u/i) 



-16/9+1 



+ -... (7.14) 



In the last equations we used the results in Appendix [ti.l| to convert the expression 
in the u variable, with the scale A defined in eq. ( |2.19| ). 

Finally, we want to write an expression for the fluctuation in the conformal frame 
scale factor and dilaton. To this end, we can use the relation (valid close to the UV 
boundary r — > 0): 



- = e 



u/i 



-u/l) 



-4/9 



(7.15) 



see Appendix |H.1| . However, it is not enough to just re-express eq. ( [7.12Q and ([7.14|) 
in terms of r through ( 7.15Q : as shown in Appendix (|A.4|) one gets an extra shift in 
5 A when changing the frame. The final result is: 

5 Air) = -Go(Ar) 4 , 5$ = -G (Ar) 4 logrA, f(r) = 1 - Y (Ar) 4 (7.16) 
5 4 

Comparison of these expansions to (|3.15|) and (|3.13| ) relates the coefficients Y 
and C defined quantities Q and C that enter in the free energy, eq ( |3.25|) as 



Y r 



4(M) 4 ' 

7.3 Asymptotics near the horizon 



5 gi 

^0 — -- 



2 (A£Y 



(7.17) 



The dependence on T on the constants Co and Yo is determined by the geometry of 
the black-hole near the horizon. Near the horizon, the black-hole solution should be 
regular. In particular one requires that the various metric functions and the dilaton 
behave as, 



f(r) 
A(r) 
$(rl 



h(r h -r) + 0(r h -rf 

A h + At(r h - r) + 0(r h - rf 

$ A + $i(r fe -r) + 0{r h ~rf 



(7.18) 
(7.19) 
(7.20) 



In terms of Y and X, the requirement of regular horizon translates into the 
following conditions: 



F($) 
X($) 



Y h 



+ Y 1 + (D($ h -f), 



~Y h + Xx($ h - $) + 0{§ h - f)\ 



(7.21) 
(7.22) 
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where Y\ and X% are yet undetermined constants and 



H=-^. (7.23) 



When one solves the coupled system of eqs. (|7.2| ) and (|7.3|) near any point $j, the 
solution will generally be parameterized by two integration constants. At first sight, 
it seems that these two parameters can be taken as and However a more 
careful look at the system of equations reveals that demanding a regular horizon 
reduces the number of parameters to a single one, that fixes Yh in terms of <&h- This 
can easily be seen by substituting ( |7.21|) and (|7.22|) in ( |7.2j ). One obtains: 



Y - 9 V '^ h) (7 24) 

Yh ~ 32 mo • (7 - 24) 

Similarly X\, Yi, etc. are determined from the sub-leading terms in the expansion 
near the horizon. See Appendix [H.8|. One finds the higher order coefficients in the 



expansion of ( |7.21| ) and ( |7.22| ) order by order. The important point is that there is 
no room for an arbitrary integration constant, the requirement of a regular horizon 
completely determines all of the coefficients in terms of This shows that we 
have a single parameter family of solutions of the X-Y system, parameterized by the 
location of the horizon 

The near horizon solution presented here is continuously connected to the near 
UV solution that is presented in Section [7]2|. This fact can easily be derived by using 
the analytic structure of the equations of motion (|7.3|) and ([7.2|) . See appendix H.4 



7.4 Thermodynamic functions and relations 

The thermodynamics is completely determined in terms of the integration constants 
Co and Y . In the following, we present the derivations one by one. 

Temperature 

The temperature of the dual gauge theory is given by the derivative of / at the 
horizon, 

T = JZM. (7.25) 
Here, we shall express it in terms of the solution of the X — Y system. The com- 



putation is straightforward and the details are presented in Appendix [H.9| . One 
finds, 

T = - (A If ^M. (7.26) 

where 6(A) is determined from X by ( |7.4|) . This equation gives T in the physical 
units of A. 
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Entropy 



The entropy of the field theory is determined by the area of the horizon as in 

S = A7cab 3 (X h ), (7.27) 



where b(X) is again determined from X by (|7T 

We note that, the equations ( |7.26|) and ( |7.27|) combine to yield the entropic 



contribution to the free energy density as follows, 

ST = AaY i 3 A\ (7.28) 
This equation clarifies the physical meaning of the integration constant Yo. 
Free Energy and Energy 

One can obtain an exact expression for the free energy in terms of the scalar 
variables. This is done by converting the eq. ( p,17| ) in A using the equations 



( [H.l| . [H.2| , [ri.3|) . This calculation is explained in detail in the appendix FT6. The 
result is expressed very simply in terms of the constants of motion defined in the 
previous section: 

T = -pV 3 = aA 4 i 3 (6C - 4Y ) , (7.29) 

where the second equation relates pressure to the free energy of the system. Using 
eq. ( [7.1 0| ) and ( |7.17| ), the above expression coincides with ( |3.25| ) 



The energy follows directly from (|7. 29|) and ( |7.28| ) as 



E = pV 3 = 6aC A 4 ^ 3 , (7.30) 

where we defined the energy density as p. This clarifies the physical meaning of the 
integration constant Co- 

Specific heat 

The specific heat is given by, 



where we used ( [7.35| ), see below. 



Speed of sound 
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The speed of sound in the medium is defined by c 2 s = dp/ dp, where the pressure 
p is given in eq. ( [7.291 ). By using thermodynamic relations, one can show that: 

c 2 s = (7.32) 



Using ( |7.31|) we can derive a relation for c s directly in terms of Yo: 

1 cHogYo 



c 2 s d log T 



1. (7.33) 



Using the high-T asymptotics of Y given in eq. ( 7.36|) , we see that c 2 s — > 1/3 for 



T/T c ^> 1 as required from the conformality in this limit. 

We refer the reader to the results presented in Appendix [I] for the high temper- 
ature behavior of the thermodynamical functions discussed in this section. 

A relation between Co and Yo 

One can also relate Co and Y , by using the definition of the entropy, 

(T.34, 



It follows from ((T3§, ( fTT29|) and (|3~26| ) that 



C (T) = ^Yo(T) - \ [ T ?Mdt. (7.35) 



Therefore, knowledge of Yo as a function of T determines the coefficient Co analyt- 
ically. 

Equation (|7.35|) also helps us determine the thermodynamics at high-T. As T 
increases, approaches zero, hence the geometry of the black-hole becomes the 
geometry of an AdS black-hole. For the AdS black-hole one has, b — £/r, f — 
1 — (r/r^) 4 and T = Xj-nr^. Therefore from the definition of Y in |7.1| ), of Yo in 
( |H.21[ ) and the conversion between A and r coordinates near the boundary ( |H.7| ) we 



obtain, 



ttT\ 4 T 



Y (T)H-H jT>1- ( 7 - 36 ) 



Using this in ( |7.35| ) one obtains, 

C " ~ T J° /2 - | » 1. (7.37) 

We present the explicit high-T behavior of this function in Appendix Q. We also 
show that this implies, through eq. ( [7.10Q , that the gluon condensate divided by T 4 
vanishes at high-T. 
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A useful equation for s/T 3 



One can also obtain a relation between the dimensionless thermodynamic ob- 
servable s/T 3 and the dilaton potential as follows: 

s An 2 (12M p £) 3 e-^o h f^W 4n 2 (l6M p £) 3 [W{\ hj , 

(7.38) 



T 3 £ 6 V(X h ) 3 27£ 3 \V(X h 

In the second line we used the definition of the thermal superpotential in terms of 
X, eq. ( p. 11| ). We refer to App. [H. 11| for the derivation. 



This relation is valid for an arbitrary background. In our case, where the small 
extreme high T limit is free YM theory, the Planck mass is given by Q5.8D . As a 
consistency check we find in the case of a conformal theory, where V = 12 / £ 2 and W 
are constants and X = 0, that ([7.38 ) reduces to the conformal thermodynamics: 



^ = 4vr 4 (M/) 3 , (7.39) 



which again agrees with the free YM gas upon use of 

Using ( [7.38Q together with the relation between T and &h ( |H.68| ), one can di- 
rectly relate the thermodynamic observable s/T 3 as a function of T/A to the dilaton 
potential V($). We note that this relation side-steps the construction of the con- 
stants of motion Co and Y described above, a procedure that is numerically quite 
non-trivial. As s/T 3 as a function of T/A is one of the main observables that can be 
used in comparison with the lattice data in principle, (|7.38| ) is quite utile 27 . 



7.5 Matching the zero-T solution 

As we discussed above, the zero-T solutions of the X-Y system correspond to the 
special case Y — 0. We analyzed the entire set of solutions in this case, in Appendix 
|H.3| . The conclusion of this analysis is a rephrasing, in terms of the scalar variables, 
of the general classification in terms of the superpotential that we have discussed 
in Section 5.2.1| and in Appendix [E[ For any dilaton potential asymptotic freedom 



in the UV i.e. V(A) — > Vq + V\ A + • ■ ■ , as A — > and which exhibits exponential 
asymptotics in the IR, V(X) — > A 2< ^(logA) P as A — > oo, there are three different 
classes of solutions to X, with different IR behavior (as A — > oo): 

i. Solutions with "special" type of IR asymptotics: X — > — 3Q/4. We denote it 
A*($), 

ii. Solutions with "generic" type of IR asymptotics: X — > — 1, 



27 However, one usually uses s/T 3 as a function of T/T c to compare with the lattice data. In 



this case the utility of (7.38) is reduced as one needs additional knowledge of T c /A, which in turn 



requires the free energy, i.e. knowledge of Co and Y0. 
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iii. Solutions with "bouncing" type of IR asymptotics: X — ► 0. 

The second case is not desired because the fluctuations of the bulk fields fall into the 
singularity, hence it is not of repulsive type, whereas the last case corresponds to A 
being a non-monotonic function of the RG scale, hence it can not yield a sensible 
RG flow. Therefore we based our holographic construction on the special class, case 
i. 

Now, we consider the black- hole solutions to the same potential, with Y ^ 0. A 
priori, there is no guarantee that a regular black-hole solution does not correspond to 
cases ii or iii as the deformation is taken away i.e. Y (by sending the BH horizon 
to the singularity). However, as presented in Section p>.2.1| , we have the following: 
No-go theorem: The only vacuum solutions of Einstein- dilaton gravity, with expo- 
nential asymptotics given above, that can be continuously lifted to a regular black-hole 
correspond to the special class of solutions, case i. 

In the language of scalar variables, the proof is simple. In the previous subsection, 
we gave the condition for regularity of the horizon, eq. Q7.24Q . We can write this in 
terms of X($) using eq. ( |7.22j ) 

= (7.40) 

For all regular BHs one should be able to push the horizon down to the singularity 
of the zero-T solution ro, by continuously sending — - > vq. First, we rule out case 
iii: Since V(<$>), V'($) > f o r all one finds from flT^Ol) that X($ h ) < 0. On the 



other hand, in the region close to the singularity in case iii, one finds X > 0, see 



Appendix [H.3| . Thus, it is not possible that regular black-holes can be continuously 
connected to case iii, as the horizon is taken close to the singularity, ^ — > oo. In 
this limit, ( [7.40 ) clearly fixes X — > —3Q/A. Therefore, the function X exhibits the 



desired asymptotics of case i. 

On the other hand, one can solve for Y in this asymptotic region, using the 
analytic solution of ([H.10|). One finds, 



m - c 6 _ d _ e ^ ™ *h » 1 (7-41) 

Here c > and d > are given in terms of Q and the location of the horizon is given 
by the integration constant C\ as 

d = -e c *\ 
c 

We can show that, in the limit $^^00 the entire Y function becomes a spike 
centered at $ = 00: From ( |7.3| ) we observe that the RHS is positive definite, hence 
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Y is a monotonically increasing function of $. From ( 7.41 ), it is also clear that 
Y($h) — > as $h — > oo, for any $ ;§> 1 but $ < oo, in the asymptotic region. 
Combining these two facts, we learn that Y($>) should vanish for all $ ^ oo, including 
the UV region. On the other hand, it diverges exactly at <E>/j = oo. 28 

Thus we proved that, in the limit &h — > oo of any regular BH, X limits to X* 
and Y(<&) — > for any $ < oo. This corresponds to the zero-T solution Y — 0, with 
the integration constant of X equation tuned to X(oo) = —3Q/A. In other words, 
X($) — ► in the entire range of $ as $^ — > oo. 

7.6 Parameters of the solutions 

Here we examine the integration constants in the Einstein equations for a generic 
black-hole solution. We solve the system by requiring the asymptotic behavior of 
X and Y near the horizon, as discussed in Section |7.3| . This solution flows in the 
UV to X — > X and Y — > as described in Subsection |7.2| . Thus one has a single 
integration constant $^ from eqs. ( [7.2[ ) and (|7.3|). It determines the temperature by 
eq. ([□§). 



After X and Y" are determined, the metric and dilaton are obtained using eqs. 
(EnD -flEPD. The condition g -»• near the boundary fixes the integration constant 



in (|H.3|) . The two remaining constants are Aq and A Q . As described in [28 1, these 



combine to determine the mass scale A of the physical system, ( [2.19| ). This combina- 



tion can be viewed as one integration constant of the two equations (|H.1| ) and (|H.2|) . 
(The other one is irrelevant due to a shift symmetry in r, see ( ||28|| ).) As we require 
that the finite-T solution approaches to the zero-T on near the boundary, the value 
of A£ is determined by the corresponding value at zero T. 

Thus, we conclude that one has only two parameters in the solutions, A and Xh 
which corresponds to the Aqcd and the temperature. Furthermore, one of them i.e. 
A is completely fixed by the zero-T solution. It is practically set by the mass of the 



lowest glueball in the spectrum, see J28 
8. Outlook 

In this paper we have analyzed in detail the equilibrium thermodynamics of the 5D 
Einstein-dilaton system that was proposed in |28| as a phenomenological holographic 



dual of 4D large-A c pure Yang Mills. There is a variety of possible directions to 
extend our work. 

An example of an explicit background with similar asymptotics, in critical (IIB 
compactified to 5D) or a non-critical string theory would be desirable. In addition 
to justifying our phenomenological set-up based on the principles of AdS/CFT cor- 
respondence, this would allow more detailed studies on the a' corrections and how 

28 In fact, one can prove that Y(&) becomes proportional to 6(&h — <£") in the limit $^ — * oo by 
using a limit representation of the delta function. We will not need this here however. 
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they can affect our results especially in the UV. For the thermodynamics of QCD, 
the higher derivative corrections are desirable also for a phenomenological reason: It 
is well-known [[50] that rj/s is constant in any gravitational theory based on a two- 
derivative action. However, this quantity is expected to be a non-trivial function of 
T in QCD that becomes asymptotically large for large temperatures. The higher 
derivative corrections may provide the desired T dependence. 

Having set the general construction in this paper, a natural step forward is to 
compute dynamical observables (bulk viscosity, drag force, jet quenching parameter) 
that are important for the physics of the RHIC collider and the upcoming LHC 
collider . We will address this problem in the near future, p4 |. Another related issue 
is the computation of the various Debye screening masses, where a better comparison 
with lattice data can emerge. 

Another important direction involves the meson sector, that should be introduced 
through probe -D4 + -D4 branes in the background. Introduction of baryon chemical 
potential is the next very interesting step to analyze. We expect that this will 
involve the study of charged black- holes under the overall U(l) gauge field of the 
flavor branes. 
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APPENDIX 



A. Various forms of Einstein's Equations 

We use a metric signature (— , +, +, +, +). We start from the action: 



S R = -M- 



J d 5 xy/g 



R- -(<9$) 2 + \/($) 
3 



+ 2M 3 I d^x\fh K (A.l) 

JdM 



with 



K, 



/IV 



-V„n„ = -n"d n h K = h ab K ab 



(A.2) 



where h a b is the induced metric on the boundary and is the (outward directed) 
unit normal to the boundary. E.g. if r denotes the AdS conformal coordinate, 



!h~r \dr ) 



5\ 



(A.3) 



The sign of the bulk term is chosen in such a way that 1) in the Euclidean regime, 
the scalar field kinetic term is positive definite and 2) the curvature of Euclidean 
AdS is negative. With this choice, the sign of the Gibbons- Hawking term is fixed, 
as usual, by the requirement that the variation of the action does not contain metric 
derivatives. 

Einstein's equations are: 



dV 



«9 M $^$--((9$)V 



- = 0, 



with the Einstein tensor defined as, 

1 

A.l Conformal Frame 

Consider the following ansatz for the metric and dilaton: 

dr 2 



ds 2 = b{rf 

with Einstein tensor, 
36(4/6 + 6/) 



f(r) 



f(r)dt 2 + dx l dx l 



$ = $( r ). 



(A.4) 
(A.5) 

(A.6) 



(A.7) 



E rr = 



26V 



3/(6/ + 2/6) 66/ + 6/6 + 6/ ^ 

&tt — ^7 j j &ij — 7~i Oiji 



26 



26 



(A.8) 
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Laplacian, 



and Dilaton stress tensor, 



T 



2 • o b 2 2 f 2 ■ „ 6 2 / 

T rr = -$ 2 + — y , T« = — $ 2 -V 

rr 3 2/ ' 3 2 

The equations of motion are: 



b 2 b bf b 2 b 2 b 

6 V 2+3 b + 3 bJ = J V > ' 6 ^" 3 6 



2/., b 2 
— $ 2 + -V. 



1$ 2 

3 ' 



(A.9) 



(A.10) 



(A.n; 



/ b -iff b \ ■ 3 l2 dV 



(A.12) 



The second equation in (|A.12|) is not independent of the other three, and it can 
be dropped. 

The Ricci scalar is: 



2 p / 
'3 E = ~V 2 



f b bf b 2 

7 + 8 6 + 8 6? + V 



R 

A. 2 Domain- wall frame 

We define, 

b = e A , f = et>, 
and use the domain-wall parametrization of the metric, 

dr = e~ A du. 

In this coordinate frame the metric has the following form: 

ds 2 = e 2A (-fdt 2 + dx 2 ) + ^. 



(A.13) 



(A. 14) 



(A.15) 



(A.16) 



The equations of motion (|3.5| ) and (|3.6|) in the variable u take the following form: 



12A' 2 + 3A'g' - -$' 2 - e~ 9 V = 0, 
3 

A + U' 2 = 0, 

g' + ^ + 4A' = 0, 
9 



8 d$ 



0. 



(A.17) 
(A.18) 
(A.19) 
(A.20) 
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A. 3 Dilaton frame 

This frame uses $ (or A = exp $) as the radial coordinate, and it is in some sense 
"maximally gauge fixed," since only the physical integration constants of Einstein's 
equations appear in the metric. To change variables from u to A it is useful to define 
a superpotential W(Q), such that eq. ( |A.18| ) is written as: 



4 dW 
A' = -gW($), = (A.21) 

The coordinate change u = is obtained by inverting the second equation in 

( |A.21|) . The existence and properties of the superpotential in the zero-temperature 



and black hole case will be extensively discussed in Appendices || and respectively. 
The solution of Einstein's equation in this frame is: 

e —9{®) rich 2 

ds " = -m uZ + e 2M * } {-e^dt" + dx 2 ) , (A.22) 

where with a slight abuse of notation we have written g(u(§)) = g(&) and A(u(&)) = 
v4($). One interesting property of the setup we are discussing is that two different 
dilaton potentials V(§) and V(Q), related by V"($) = V(<&+K) for some constant K 
give essentially the same set of solutions and the same physics. First, given a solution 
(A(u),g(u),$(u)) of eqs. (|A.17| - |A.20D with potential it is straightforward to 

show that the functions (A(u), g(u), = (A(u), g(u), $(it) — K) solve the system 

with potential V(&). The two solutions are physically equivalent (except for a change 
in initial conditions, which as we know [p| only affects the overall scale). This is 
easily seen writing the second solution in the $-frame: 

ds 2 = e ~ 9md($>2 + (_ e m dt 2 + dx 2) (A . 23) 

(r);W) 2 V / 



where W is the appropriate superpotential. However, since A{u) = A{u) and g{u) = 
g(u), it follows that 

i($)=A($), W(Q)=W(§), g($) = g($). (A.24) 

Thus, after a change of coordinates $ — > $ = $ + K the metric (|A.23|) becomes 
identical to the solution of the original system with potential V(Q), eq. ([A.22 ). The 



initial conditions for the two systems, such that the solutions coincide, are related 

byi($o-^) = A($o). 

Thus, there is an "accidental degeneracy" in the classification of Einstein-dilaton 
gravity by the dilaton potential, since the two potentials V(X) and V(k\) lead to 
the same physical results. For this reason, the value of the proportionality constant 
between the dilaton A and the physical Yang-Mills coupling X t is irrelevant. 
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A. 4 Relating fluctuations in different frames 

In this appendix we work out the relation between the scale factor and dilaton 
fluctuations close to the boundary in different frames. For simplicity we set i = 1. 
Suppose we start with the zero-temperature and black-hole metrics, both in the 
domain wall frame: 

ds 2 = < ^- + e 2AU (fdt 2 + dx 2 3 ) , ds 2 = du 2 + e 2A " (dt 2 + dx 2 ) (A.25) 

For clarity, we have added a label (it) to the warp factor. If the two solutions obey 
the same boundary conditions at u = — oo, then as shown in Section \!.2\ the two 
scale factors are related, to lowest order, by: 

5A U = A u {u) - A u (u) ~ G u e 4n (A.26) 

for some constant Q u . The difference between the dilatons, <5$" = — 
can be related to 5 A by perturbing equation ( |A.18| ), which gives: 

(5A U )" + l($ u o y(5<f> u )' (A.27) 
9 

which can be integrated to give: 

6® u ~ -g u (-u)e 4u , (A.28) 

in agreement with eq. ( |7.14j) . 

Now we want to obtain the same quantities, namely A — A Q and $ — $ Q , in 
conformal coordinates. Naively, one may think that it should be enough to make the 
replacement u — > logr in eqs. ( [A.26Q and ( |A.28|) , since any correction to the relation 
between u and r would only affect higher orders in 5 A and 5§. This is however 
incorrect, careful analysis reveals. 

What we want to obtain is the difference SA r = A r (r) — A r (r), where the con- 
formal warp factors are such that the metrics have the form: 

rfs W^ + ^ + 4 * = + # + (A,9) 

Now, it is clear that to bring the two metrics in this form one needs two different 
coordinate transformations. We can first define e~ Au ^du = dr: if we perform this 
coordinate transformation on both metrics we get: 

ds 2 = g 2A«(«(r)) ^ + fdt 2 + dx 2^ ? ^2 = e 2A^u(r)) + 2SA^u(r)) ^2^^)) ^2 + ^2) 

(A.30) 
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So the black-hole metric is in the conformal form, but the zero-temperature one is 
not. From the above expression we read off that the function A r (r) in ( A.29Q is given 
by A u (u(r)), but a similar relation does not hold for A r . 

Let us define A r (r) = A^(u(r)). To get the correct scale factor we have to 
perform a further coordinate transformation to bring ds 2 in conformal form. We 
thus define: 

(l+5A u (u(r)))dr = dr => r = r + \g u r 5 . (A.31) 

This transformation brings the metric ds\ in conformal frame, parametrized by the 
coordinate f, with scale factor given A r (r) = A r (r(r)). Using the explicit form of 
r(f), we can write this as: 

A r (f) = A r (r - G u f 5 /5) ~ A r (f) - %- f 5 d f A r (?) = A u {u{r)) + ^f 4 (A.32) 

5 5 

In the last step we have used the fact that, to lowest order, A r (r) ~ — logr. The 
first and last steps of the above equality mean that there is an extra shift in A Q , 
compared to the naive change of variables. Renaming f — > r in ( |A.32| ), we finally 
have: 

5A r = A r (r) - A r (r) = A r (r) - A*(u(r)) - — r 4 = -Q u r 4 (A.33) 

5 5 

This result could have been guessed from the fact that 5 A is not invariant under 
linearized diffeomorphisms r — > r + £, but rather it transforms as S A — > 5 A + A Q ^. 
For ^ = Q u r 5 /5, this gives exactly eq. ( |A.33| ) 



Following the same procedure, one can write 5$ in conformal coordinates, and 
one would find that this time the leading order is simply given by the change of 
variables u — > logr. The reason is that, under r — > r + £, the dilaton fluctuation 
transforms as 5$ — > 5$ + $ £ = 5$ + G(r 4 /5)(logr) _1 . Thus the shift induced by 
the extra diffeomorphism is subleading with respect to the first term, which behaves 
as r 4 logr. Thus, we have to leading order: 

<5<F(r) ~ 8<$> u (u(r)) = ^"r 4 logrA. (A.34) 

Finally, to connect this discussion with the definition of Q given in Section 3, let 
us define Q = 4Q u /5. Then we arrive at: 

5A r (r) = Gr A , 5$ r = — Or 4 log rA. (A.35) 

8 

For consistency, one can check that the above fluctuations solve the linearized 
Einstein's equations in conformal coordinates, obtained by perturbing eq. ( A.ll| ). 
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B. The AdS§ case revisited 



In this appendix we will reconsider the Einstein system plus a scalar in the conformal 
case, with a view of exploring all potential boundary conditions at infinity and their 
effect in the bulk physics both at zero and finite temperature. This situation is 
simpler than the one we are studying but some of the lessons are similar. Although 
what we present here is mostly understood in the literature (see for example f42, f43fl), 
they are not widely known and we would like to put them in the right perspective. 



B.l Zero temperature 

We reconsider the zero temperature field equations (in the conformal coordinate 
system) of the Einstein-dilaton system with a potential 



b 2 b 

\ 2+ h 



&V($) , 6 



b 2 



1$ 2 
3 



- b- 3 l2 dV 
b 8 <i$ 







(B.l) 



(B.2) 
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We will set the potential to be a constant V 
of solutions for arbitrary initial conditions. 

The first equation can be integrated once to yield 



and we will find the UV asymptotics 



b= - 



C 2 6 4 



(B.3) 



where we have chosen the minus sign branch so that b decreases with increasing 
r, and C 2 is an integration constant that we set to be positive. When C = the 
solution is AdS 5 

b = — ^— (B.4) 
r — ro 

The constant C, and two extra boundary conditions for the two first order equations 
( B.3|) and the second one in ( B.l|) viewed as a first order equation for the dilaton are 
the full set of boundary conditions necessary near the boundary. 

For general C, the first order equation (|B.3|) can be integrated as 



du 



r — r 



VI + C 2 £ 2 u 8 



(B.5) 



giving 



1 1 9 l 2 C 2 



8' 2' 8' 



6 s 
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i C 2 
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+ 



r 



(B.6) 



where F is the standard hypergeometric function and b* 
the dilaton becomes 



b(r ). The equation for 



$ = + -ArcTanh, / 1 + 



(icy 



where the sign ambiguity is hidden in the sign of C. Note that for the AdS solution 
the dilaton is constant. The three integration constants C, 6*,$o are in one-to-one 
correspondence with the three boundary conditions at the boundary. As explained 
in ||28|| , one of them, that we can take to be fe* is a gauge artifact, related to the 
position of the boundary in the radial coordinate. We will therefore choose ro = 
to be the position of the boundary, 5* = oo. Then the solution becomes 

"1 1 9 £ 2 C 2 ' 



;'2'8' 



b 8 



-ArcSinh^-f 
4 b 4 



Near the boundary 



6 s 



oo and we can expand F around zero to find 

C 2 r 8 



1 - 



18£ 6 



O (r 16 ) 



valid when r 



0. In the same region 
$ = + 



(B.8) 



(B.9) 



(B.10) 



Therefore for non-zero C there is a non-zero vev of the operator dual to the scalar 
Consider now the region b — > 0. We use the transformation rule 



F 

to obtain 



1 1 9 

8' 2' 8' 



r[|]r[|] b 16- 

sr[|] (cey 



3C£ 



1 3 11 

2' 



8 ' £ 2 C 2 



r 



r [j] r [j] 
sr[i] 



$ ~ log(f 



b ~ (3C(f - r) + • 
The scalar there diverges as 

We may use the relations 

6 _ 26 2 

to calculate the curvature invariant for the metric. We obtain 

12C 2 



£{C£) 



2C 2 
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b 8 
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(B.H) 

(B.12) 
(B.13) 

(B.14) 
(B.15) 



Near the boundary, b — > oo and we obtain constant negative curvature as expected. 
In the interior, as b — > 0, we observe that the space has a curvature singularity, if 
C 7^ at a distance 5r ~ (C£)~^ from the boundary. Imposing regularity in the bulk 
imposes C = 0. Therefore the dynamics of the theory does not allow for a non-trivial 
vev associate to the operator dual to $. 
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B.2 The black-hole solution 



We will now solve again the equations with constant potential seeking a black-hole 
type solution. 

A • 



b 2 b bf 12b 2 



lr b 
6 ^" 3 6 = 3* 



l + 3j = , $ • | t -V I 1 



We may integrate once the two last equations as 



/ = - 



8C 



3D 

6V 



(B.16) 



(B.17) 



(B.18) 



and we will take the two constants of integration to be positive D > 0, C > 0. Using 
this may derive the following equations for b 
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(B.19) 
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We now introduce an auxiliary variable 
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(B.20) 



b i + 2Cb 



% + 2Cb) 



(B.21) 



where e = 1 iff ^ + 2Cb < and e = — 1 if |j + 2C6 > 0, (6 is always negative). In 
terms of this new variable we obtain 



I -2- 
b ~ b 2 



A 



C+l 



2£ 2 (eD - C^/C 2 ^!) 



(B.22) 



which may translated as a first order equation for £ 



(B.23) 



where the prime stands for derivative with respect to b. This can be integrated as 

= Cb 8 (B.24) 



t- 1 \VC 2 + D 2 + eC + D^J 
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Using the relation 



(B.25) 



the solution can be written in the following alternative form 
1 C+l {^TW-eC-D^ 



2D(- 1 



In particular, for e = 1 



—s— = Cb 8 (B.26) 
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1 " J ' VC 2 + D 2 -C - D x i — M [ VC 2 + D 2 + C + D* i — ^) =C , 6 8 



2D(-1 \ V C + l/ \ V C + l 
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(B.27) 



c c 

1 + ; > , a_ = 1 , > (B.28) 

VC 2 + ~ VC 2 + D 2 ~ v ; 



while for e = — 1 



a+ 

= C6 S 



(B.29) 

We will now investigate several special cases of this general solution. 



B.2.1 C = 

The solution becomes 



c-i 

which using ( |B.21| ) becomes 



/} Cb s (B.30) 



b 2 = ^—±- + ^—b* (B.31) 

AC 2 £ 4 b 4 2CD£ A 

Compatibility with the other equations determines C = This solution has no 

horizon (f=l). It is the same solution found in the previous subsection, where D 
plays the role of the $ condensate. 

B.2.2 D = 

The equation (P-26| ) becomes the trivial one £ = 1. From ( p.20|) we obtain I = 2^ 
which is solved by the AdS scale factor b ~ -. Finally the rest of the equations give 

p 9rv 4 



-71- 



This is the standard AdS black-hole solution with a flat horizon. 
The two previous cases indicate that the constant C controls the temperature of 
the solution while the constant D controls the <£> condensate. 



B.3 Analysis of the general solution 

The function / is given by 



/ 



2D£ 2 (D - eC^/C^T) 



(B.33) 



The boundary b — > oo, is always at £ = 1. To test whether we have a regular horizon 
we need the the trace of the Einstein tensor which is given by 



E = — 



b 2 



f b bf b 2 



2 D(3( - 13) + lOeCyC 2 - 1 
P D-eC^C^l 



(B.34) 



To continue we distinguish two cases: 



B.3.1 e 



There is no horizon in this case as / never vanishes. There is also no singularity as 
E never blows up. The scale factor becomes a constant asymptotically 



B.3.2 e = 1 

The horizon (/ = 0) is at 



(B.35) 



It is singular unless D = 0. At the horizon 6^0. Therefore, even at finite temper- 
ature the $ condensate must vanish in order to have a regular solution. 



C. The Black hole action and ADM mass 



C.l The on-shell action 

We want to compute the regularized action evaluated on a solution of Einstein's 
equations, ( |A.4| ). We start from eq. Q2.1| ) : 



S 5 = S E +S GH , S E = —M / <Pxy/g 



J R--(a$) 2 + y($) 



, S GH = 2M 3 / d 4 xVhK. 

JdM 

(C.l) 

We work in conformal frame, with the metric given by eq. ( |A.7| ). 



Taking the trace of equation ( |A.4| ) we obtain for the Ricci scalar 



R= t m2 ~l V - 



(C.2) 
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This leads to the on-shell Einstein action: 



2 f 2 f@ r h 
S E = -M 3 / d 5 x^ V{<$>) = -M 3 V 3 dt dr b 5 V{$) 

3 J 3 Jo J 



2M 3 V 3 I dt P dr ^-(fb 2 b) 



o 



o 



dr 



From ( |A.3| ) the components of the unit normal to the boundary are n r 
Hi = 0. The trace of the extrinsic curvature is: 



K 



V7 

2b 



b f 



(C.3) 
(C.4) 



We find: 



S E = 2M 3 V 3 I dt 



" dr ^(fb 2 b) = 2(3M 3 V 3 {f{r h )b 2 {r h )b{r h ) - f{e)b\e)b{e)) 

(C.5) 
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(C.6) 



Putting together eqs. ( |U.5| ) and ( |C.6| ) we obtain for the full result: 



S e = S E + S e GH = 2(3M 3 V 3 



3b 2 (e)f(e)b(e) + \f{e)b 3 {e) 



(C.7) 



where we used that /(r^) = 0. 

The calculation for the zero-temperature background is exactly the same, except 
that the integral of the Einstein-Hilbert action extends on the region (0,r ), where 
r is the singularity. Thus in evaluating ( |U.5| ) one gets: 



S E = 2pM 3 V 3 (b 2 o (r )b o (r ) - b 2 {e)b {e)) 



(C.t 



The IR contribution vanishes whenever b 2 b —>■ as r — > r . This is always true for 
good singularities. 

C.2 Evaluation of the free energy 

We start from eq. fl3.24Q , which we rewrite below: 



T = a 66 2 (6)V75 



b(e)yW)- 



b 2 (e): 



bl(e 



+ f(e)b 3 (t 



(C.9) 
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where the limit e — ► is understood. In terms of 5b = b — b , and 5e = e — e, the 
previous equation reads: 



a 



66 2 (e)v^ 



(bo + 5b) (e) yW) - r f° + /f, ( '\ (bo + 5eb )(e) 
(b + 5eb o y(e) 



+f(e)(b +5bf(e) 
(CIO) 



For small e, b (e) ~ e/r, and by eqs. Q3.13| ), ( p,15|) and (|3.21|) we have: 



.e 6 . „3e 2 x . Ce 4 . 45 £ f _ 2 



Sb(e)^Q^ 5b(e)^G—, /( e )~l- T _ 5e = --^e & (logeA) 2 . (C.ll) 

We see that the only divergent contribution inside the square brackets, i.e. fo , cancels. 
What is left is of order e 2 times eventually some logarithmic corrections. Therefore, 
to this order we can replace the overall prefactor b 2 {e)^J '/(e) by £ 2 /e 2 . Thus, to 
lowest non-vanishing order: 



a e l 



c 



(C.12) 



The last term in the parenthesis requires more care: due to the extra logarithm in 5e 
we cannot just replace b Q (r) by ijr. On the other hand we can use the zeroth order 
Einstein's equation ( |A.11| ) to write it as: 



4 • 



bo(e)5e 2 g(e) - £(e) ) = b (e) 5e -& (e 



b 2 , 
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(C.13) 



where in the last line we used $ (e) = — (elogeA) 1 (cfr. eq. (|2.18|) ). Notice that 
the logarithm in 5e has canceled. Finally, we get: 



a 4 



(C.14) 



C.3 The black-hole mass 



The mass of a solution, with respect to a reference background, can be defined 
following the procedure: first consider a time slicing of the 5-dimensional metric, of 
the ADM form 



ds 2 



N 2 dt 2 + 7ij(dx l - N'dt^dx 3 - N 3 dt) i,j = r, 1, 2, 3 



(C.15) 



where 7y is the induced metric on each 4D slice E t . Then the mass of the solution, 
with respect to a reference background with the same asymptotic behavior at spatial 
infinity, is given by j|6 



E 



8ttG 



5 Jt, 



(C.16) 
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where T^ is a 3- dimensional surface at spatial infinity embedded in the 4D constant- 
time slice T t , "find is the three-dimensional induced metric, ( 3 'K is its extrinsic curva- 
ture, and 7* nd and ^'K the analogous quantities for the reference background. The 
latter should be chosen so that the geometry of the 3-dimensional surface at infinity 
and the value of the scalar field on that surface match. 

In our case, the 5D solutions are static and in conformal coordinates they are of 
the form: 

ds 2 = b(r) 2 (-f(r)dt 2 + ^- + dx m dx m ) , ds 2 = b (r) 2 (-dt 2 + dr 2 + dx m dx m ) . 
V fij) J 

(C.17) 

The boundary at infinity is at r = e, with e — > 0. Thus, we have: 

r- ( dr 2 \ 

N = 6(r)v7, j ij dx i dx j = b{r) 2 —— + dx m dx m , (C.18) 



fir) 



and on the surface r = e we have 



^ = He) 2 6 mn , „' = (0 (C.19) 
The 3D extrinsic curvature is given by: 

= v<n 4 = (vry%) = -3^^- (G.20) 



The reference background has / = 1, boundary at e = e + 5e (as in |3.21|) so that 
A(e) = A Q (e), and rescaled volume V3 = V3& 3 (e) /b (e). Also, the time slicing has to 
be the same, i.e. N(e) = N (e) in the ADM decomposition. Thus from eq. ( |C.16| ) 
we obtain: 

Using eqs. ( |U.11[ ) and performing similar steps to those that led to eq. ( |C.14| ) we 
obtain: 

E = M*N 2 V 3 (150 + |c J (C.22) 



where we have used 167rG 5 = M„ N r 



D. The gluon condensate asymptotics 



We want to show that the gluon condensate Q obeys to the asymptotics ( |5.5| ) at high 
temperatures. To do this we need to compute explicitly the relation between T and 
Th in the UV limit r n — > 0. 
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Temperature Knowing that the expansion for the metric reads as in ( 2.17|) (this 
is the zero temperature solution, however the r 4 correction of the finite temperature 
solution is subleading w.r.t. logarithms) we find that the expansion for the thermal 
factor f(r) in the UV is given by solving the second equation in ( |3.5| ) 



m = i - 



i + 



3 log Ar log Arh 
The derivative of / w.r.t. r then evaluates to 
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The temperature is obtained by evaluating the above expression at the horizon 
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Entropy We now want to calculate the entropy density s = 47rM 3 6 3 (r^). Inverting 
this relation we get r^ as a function of the temperature 
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Plugging this expression into the expansion for the scale factor (|5.5|) we obtain b as 
a function of the temperature 



b(T) = ti£T 



I 1 m / log(- log 4;) 

1 + T- ■ O ' 



91og^ V W£ 



(D.5) 



The subleading term O (log 2 (A/7rT)) will not enter into the leading calculation of 
the gluon condensate asymptotics. The entropy density evaluates to 



s(T) = 4ti 4 £ 3 T 3 



I 1 m I log(- log 4) 



4tt 4 £ 3 T 3 £(T) (D.6) 



where the last equality is the definition of the function £(T). 
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Gluon condensate asymptotics Putting together the information relating the 
gluon condensate to the free energy on the one hand through eq. (|3.25| ) and, on the 
other hand, the free energy to the entropy through JF = —dS/dT we arrive to an 
equation for the gluon condensate 



ug{T) J^p- - jf dt s(t) 

=7rV 3 T 4 £(T) - 4ttV 3 f dt t 3 £(t) (D.7) 

Jt c 
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f dt tr(t). 



The last line uses integration by parts. The expansion of the derivative of the function 
£(T) reads 



e / (T) = --^ ir (l + ...) (D.8) 



4 1 1 

The ellipsis indicates subleading terms in the log expansion. 

So that finally the gluon condensate expansion at high temperatures T > A at 
leading order can be written as 



36 log ^ 
E. The Superpotential at zero-T 

Here we analyze the general solution of the zero-temperature superpotential equation, 
eq. (below, A = e*). 

4 64 
--X\W\X)f + *-W\X) = V{X). (E.l) 

we assume V{X) > First let us observe some general properties: 



1. The solution can only exist as long as \W(X)\ > y/(27/64)V(\); 

2. The equation has a symmetry W — > —W, so we can limit the analysis to 
W > 0. 



3. For any Ao ^ 0, there are two solutions of (|E.1| ), W + (X), W-(X) passing through 
the point A , such that W+(X ) = W_(A ), and W' + Xq) = -W!_(X ). In other 
words there are two branches of solutions: one where W and W have the same 
sign (i.e. WAX) > 0) , another where they have opposite sign (W + (X) < 0) 
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4. 



At any A* ^ where |W(A*)| = ^/27/64y(A*), W = 0. 



5. A solution can go past such a point A* only if V'(A*) = 0. Indeed, suppose that 
y (A*) > 0. if the solution exists for A < A*, at the point A* we have: 

W(K) = v / (27/64)^/(A,); W\K) = ; V'{\*) > ^(A*+e) < v / (27/64)V(A* + e) 

(E.2) 

therefore the solution does not exists for A > A*. 

6. By the same argument, if V^(A*) < 0, the solution does not exist for A < A*. 

7. It follows from points 3,4 and 5 that, if V(A) is positive and monotonic, the 
two branches W+(\) and W-(X) (see point 4) are completely disconnected, 
since neither W nor W can change sign. However two solutions belonging to 
different branches can be glued together at a point W' = 0. 

8. All solutions that reach A = have either W(0) = ^(27/64)^/(0), or W'{0) = 
oo. 

In what follows we assume V(A) > and without loss of generality we take 
W(X) > 0. 

E.l Solution close to a critical point 

Let us see how the solution approaches the critical points W'(\*) = 0. For defmite- 
ness, consider V(A) monotonically increasing, and W(X) > (as in our model). The 
solution exists only for A < A*. As we said there are two disconnected branches with 
opposite sings of W . Let us are analyze them separately. 

W- Branch 

In this case, W'(X) < for all A < A*. Then eq. QE.1J ) can be written as: 

W'(X) = -^Wi-f^V (E.3) 

Let us look at this equation close to a point A* ^ where W 2 (\*) = 27/64V A (A*) = 
W?. Write W = W* + w(A), with w(A*) = 0, and expand QE.3Q to linear order in 

w(\y. 

w'(X) = -^^/2WMX) - V'{\*){\ - A,) (E.4) 

OA* 

this is still hard to solve explicitly, but we can carry out the analysis by making 
further assumptions about the possible behavior of w(X) close to A*. There are only 
three possibilities: 
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1. \w(X)\ > 0(X - A*) as A —> A*. 

then we can neglect the second term under the square root, and the equation 
becomes 

w'(X) = - — y/2WM*) (E.5) 

which is solved by w(X) ~ (A — Z*) 2 . This is inconsistent with the assumption 
|w(A)| > 0(X — A*), so this case is ruled out. 

2. |w(A)| ~ —Wx(X — A*) as A — ► A*, with w\ > . 
Eq. (|E.4|) becomes: 

w'(A) = --^-^(2^1 + V"(A,))(A, - A), (E.6) 

which is solved, for A < A*, by w ~ const(X* — A) 3//2 , in contradiction with our 
assumption. So this case is ruled out too. 

3. HA)| < 0(X - A*) as A -> A,. 

In this case we can neglect the first term in the square root, and obtain: 

w\x) = --Lyv'ix^-x) (e.7) 

which integrates to: 

^(A)~^v / ^ 7 (Aj(A*-A) 3/2 A<A, (E.8) 

This time the solution is consistent with the hypothesis. 

The three possibilities listed above exhaust all possible behaviors of W(X) close 
to a critical point A*, and the only one which does not lead to a contradiction is the 
last one. Thus, we can conclude that the behavior close to a critical point is given 
by eq. ( E.8|) . in other words, W-(X) is positive and decreasing for A < A*, and it 
reaches a finite value at the critical point, where it behaves like 

W_(A)~ W* + ^(A* - A) 3/2 , Wi>0 (E.9) 

The behavior of W-(X) close to a critical point is exemplified in figure |l| (a) 

W + Branch 

The analysis is the same as for the W- branch: the solution is defined only for A < A*, 
except that W(X) is increasing and close to A* we have: 

W + (X) ~ W* - Wi(A» - A) 3/2 , W 1 >0 (E.10) 

this type of solution is shown in figure ^ (b) 
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W(X) W(X) 




40 60 
(b) 



Figure 12: Superpotential on (a) the W_ branch and (b) the W+ branch, close to a critical 
point. The black area is the "forbidden" region below the critical curve a/271764, where 
W' would become imaginary. The solution stops where it meets the critical curve. 

Metric and dilaton close to a critical point 

Although PU'(A*) = 0, the metric is not AdS close to A*. In fact, the equation for 
A{\) reads: 

X dA - 3 X — 3 \ W ' (FU) 

Close to the critical point A*, X ~ ^(2>W]_/2W*)(\* — A) 1 / 2 , so the scale factor close 
to A* is finite, and behaves as 

A±(A)~ A*±Ai(A* -A) 1/2 , A 1 >0 (E.12) 

Notice that the upper sign corresponds to the W + branch (X < 0), in which the 
scale factor decreases towards the endpoint A*, which is therefore in the IR.The 
other branch W_(A) (X > 0) has A = A* as the UV. 

Finally, we can integrate the dilaton equation as a function of the coordinate r, 



Using the form of A(X) and W(X) close to A*, eqs. (|E.12| ), we arrive at: 

A(r) ~ A* - Ai(r* - r) 2 , A±(r) ~ A* ± A x \n ~ r\, Ai>0,A!>0 (E.14) 

Here the upper sign holds for r < r*, the lower for r > r*. From the last equation, 
we see that in the r coordinate the point r* where A(r ! „) = A* is perfectly regular, 
and we can obtain a full solution describing both branches by simply removing the 
absolute value, 

A(r) ~ A, - Ai(r* - r) 2 , A(r) ~ A* + - r), Ai > 0, Ai > (E.15) 

At the critical point the dilaton reaches its maximum value A*, then reverts the 
direction of running, and is not a good coordinate globally. Instead A(r) is monotonic 
along the full solution. The UV corresponds to r < r*, the IR to r > r*. 
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E.2 Solutions close to A = 

Here we still have two branches, but we have two completely different behaviors in 
each branch. Assume V > 0, V > close to A = 0, and a power expansion of the 
form: 

V = V + V 1 X + V 2 X 2 + ... (E.16) 

W- Branch 

In this case the general solution of eq. ([E. 1|) close to A = is: 

/ C A 4 / 3 A 4 / 3+1 \ 
W-{X)~W l^ 5 + -^- + W 1 - d - + ...J 7 C>0 (E.17) 

where Wq and W± are completely fixed by the expansion coefficients of V(A) around 
zero. Since W_ is a decreasing function for small A, from our general considerations 
we know that if V^(A) is monotonic then W-(X) is monotonically decreasing globally, 
therefore the solution will hit a critical point at some A* > and terminate. 

Solving the metric and dilaton equations in the region A ~ with the superpo- 
tential (|E . 1 7|) gives a singularity at a finite value r = tq of the conformal coordinate, 
where both the scale factor a(r) ^ and the dilaton A — > 0. 

e^r) „ (r -r)3, A(r)~(r -r)i (E.18) 

We see that e A ^ decreases to zero as A(r) — > 0, and l/A > 0. With our holographic 
dictionary, (log A <-> E and /3(X) = l/A, we conclude that Therefore in this case the 
small A region is in the IR, and the theory is IR-free. 

W + Branch 

In this branch, any solution of ( |E.1| ) necessarily satisfies 



W + (0) = ^V o (E.19) 

(one can show that any ansatz with W(0) = +oo cannot solve the equation). 
Moreover, if W is written as a power series expansion around A = 0, 

W(X) = W + WiX + W 2 X 2 + ... (E.20) 

then all the coefficients W{ are uniquely determined by the expansion coefficients Vi 
ofV(X). However, it is incorrect to conclude that that the solution in this branch is 
unique. To see this, take any function iy(A) that solves (|E.1|) to all orders in powers 
of A. Then, consider a function W(X) that, close to A = 0, behaves as 

W{X) = W(X) + w{X), -> as X -> 0. (E.21) 
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Inserting this in the eq. 
A = 0: 



TJ), and expanding to linear order in w (A) gives, close to 



4 



64 



-2AW'(0)w'(A) + — 2W(0)w(A) = 
i.e. a homogeneous, linear equation, whose general solution is 



Wc (\) = CA 16 / 9 - 



•4b 



exp 



16 W(0) 1 
9 W>(0)\ 



(E.22) 



(E.23) 



Therefore, also in the W > branch we have a one-parameter family of solutions, 
that close to A = all have the same power expansion, and look like: 



W(X) = W + w c (X) + 



(E.24) 



where W is a fixed power series (say, with no exponential part), w c (X) is given in eq. 
( [E.23|) , and the dots represent even more subleading terms (~ w 2 ). 

Due to the expasion (|E.20| ), the solution close to A = is, to leading order, an 
AdS§ spacetime with logarithmic running, 



b(r) 



1 + 



4 



9 log r A 



+ 



A(r) 



1 



b log r A 



+ 



(E.25) 



Notice that the exponent in (|E.23j ) is fixed by the first two expansion coefficients of 
V^(A), and one can easily show that, in terms of the /^-function coefficient bo : 



16 W(0) 4 
9 W'(0) ~ b 



A 16/9- 



lb. 



(E.26) 



Using the perturbative asymptotics &oA ~ (logr) -1 , this corresponds to a power- 
law correction to the logarithm expansion in (|E.25|) , that scales like r 4 close to the 
AdS boundary r = 0. Since the power series expansion of W(X) around A = is 
independent of the integration constant C in ( E.23j ), we conclude that metric that 
correspond to different solutions on the W + branch differ only by non-perturbative 
0(r A ) terms, which correspond to different values for the gluon condensate. 

E.3 Solutions close to A = oo 

Finally we analyze the solution of ( E. 1|) in the asymptotic region of large A. We 
assume for the potential a power-law behavior 



^(A)~y oo A 2Q (logA) J 



A 



oo 



(E.27) 



for some constant and Q > 0. We are interested in I4o > 0, since this case 
corresponds to a potential which is bounded from below. There are two kinds of 
solutions: 
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1. a continuous one-parameter family of the form: 



W c {\) = W 



CA 4 / 3 + 



C 



-i 



-A 2Q - 4 / 3 (logA) p + 



(4-3Q)' 

where C is an arbitrary constant of integration; 
2. a single solution that asymptotes as 

W s (X) = W OQ X Q (\ogX) p / 2 , W x 



64 
(E.28) 



27T4 



4(16 -9Q 2 ) 



(E.29) 



Notice if > 0, both types of solutions exist only if Q < 4/3: for Q > 4/3 the 
l.h.s. of the differential equation is asymptotically negative. In this case there is no 



solution that reaches arbitrarily large values of A, but rather all solutions to (|E 
are of the type described in section ( [E.l| ): they reach a maximum value A* where a 
W + and a W_ solutions join. 



With the restriction Q < 4/3 the first term in eq. ( |E.28| ) is the dominant one, 
and the singular solution grows slower than any of the solutions in the continuous 
family. 

For all superpotentials in the continuous family the metric and dilaton exhibit 
the same kind of IR singularity at finite r (where a(r) — > 0, A(r) — > <x>) , regardless 
of the value of Q and P: 



a(r) ~ (r - r) 1/3 , A(r) 



(ro — r) 1 / 2 



(E.30) 



This is similar to the singularity in eq. ( E.18| ), up to A — > 1/ A. These singularity 
always fall in the pathological class, as discussed in (fl2"8f): the singularity is not 
screened from physical fluctuations, and one can have an infalling flux of particles. 
Moreover, as shown in Appendix [FJ these singularity do not appear as continuous 
limits of black-hole solutions with regular horizons. 

On the ohter hand the singular solution ( |E.29| ) is the most interesting from a 
physical point of view: the singularity is repulsive for Q < 2y2/3 and it can be 
cloaked by a horizon. Solutions of this kinds are the ones that give rise to the most 
interesting holographic constructions from the QCD perspective. 

E.4 General Classification of the solutions 

The results of this Appendix can be summarized as follows: For any positive and 
monotonic potential V(A) with the asymptotics : 



V(X) = V + VtX + V 2 X 2 + ... V > 0, 
V(X) = V 00 X 2Q (\ogX) p , V oo >0, X 



A 

oc 
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the zero-temperature superpotential equation has three types of solutions, that we 
name the Generic, the Special, and the Bouncing types: : 

1. A continuous one-parameter family that has a fixed power- law expansion near 
A = 0, and reaches the asymptotic large-A region where it grows as 

W ~ C 6 A 4/3 A -> oo (E.31) 

where is an arbitrary positive real number These solutions lead to back- 
grounds with "bad" (i.e. non-screened) singularities at finite ro, where b(r) — ► 
and A — ► oo as 



a(r) ~ (r 



-rW 3 



A(r) 



(ro-r)- 1 / 2 



(E.32) 



We call this solution generic. An example is shown in figure [13 




Figure 13: Superpotential of the "generic" kind. The black area is the forbidden region 
below the curve 



2. A unique solution, which also reaches the large-A region, but slower: 



W(X)^W^(\o g X)^, Woo = J 4(16 27 _^ g2) (E.33) 

This leads to a repulsive singularity, provided Q < 2^/2/3 pSj . We call this 
the special solution. An example is shown in figure [14] 

3. A second continuous one-parameter family where W(X) does not reach the 
asymptotic region. These solutions have two branches that both reach A = 
(one in the UV, the other in the IR) and merge at a point A* where VT(A*) = 
v /27^(A,)/64. The IR branch is again a "bad" singularity at a finite value r , 
where W ~ A _4//3 , and 

b(r) ~ (r - r) 1/3 , A(r) ~ (r - r) 1/2 . (E.34) 



We call this solution bouncing. An example is shown in figure 15 
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Figure 14: Superpotential of the "special" solution. The black area is the forbidden region 
below the curve 




Figure 15: Superpotential of the "bouncing" kind. The black area is the forbidden region 
below the curve -y/27V/64 



Notice that, as two solutions with positive derivative cannot cross, the special 



solution (figure [14]) marks the boundary between the generic solutions, that reach 
the asymptotic large-A region as A 4//3 (figure [13|) and the bouncing ones, that don't 



reach it, figure |15| Notice that, if Q > 4/3, only bouncing solutions exist. 

In all types of solutions the UV corresponds to the region A — > on the W + 
branch. There the behavior of W+ is universal: a power series in A with fixed coeffi- 
cients, plus a subleading non-analytic piece which depends on an arbitrary integration 
constant C,„: 



W 



J2 + C w \ w/9 e 



16W 

W [1 + 0(A)] 



(E.35) 



i=l 



All the power series coefficients Wi are completely determined by the coefficients in 
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the small A expansion of V(A), the first few being: 



W = — — , Wi = — W — , W 2 = jj- 2 E.36 

8 16\ VO 1024^ 7 



F. The Superpotential at finite T 
F.l The thermal Superpotential 

A useful way of counting the integration constants and to parametrize the differ- 
ent black-hole solution is to extend the notion of superpotential to the black-hole 
backgrounds. We will call the resulting thermal superpotential W(<&), and its zero- 
temperature precursor W ($). 
First, we define 

b = e A , f = e«, (F.l) 
and use the domain-wall parametrization of the metric, 

dr = e~ A du. (F.2) 



In this coordinate frame the metric has the following form: 



ds 2 = e 2A (-e 9 dt 2 + dx m dx m ) + e~ 9 du 2 . (F.3) 

The equations of motion ( |3.5|) and ( [3.6|) in the variable u take the following form (a 
prime denotes derivative w.r.t. u): 

A" + ^$' 2 = 0, (F.4) 

g > + C + AA > = o, (F.5) 
9 

12A' 2 + 3A'g' - -$' 2 - e~ 9 V = 0, (F.6) 
3 



Equation (|F.4|) is the same as in the zero-temperature system, cfr. eqs. (2.6) 



This has an interesting consequence: we recall that ( |F.4| ) is the equation that guar- 
antees a well-defined arrow of the RG flow in the gauge theory. Here we see that 
one can define the RG flow in the same way in the gauge theory at finite-T. In 
other words, the conclusion that A(u) is monotonically decreasing still holds at finite 
temperature. 

Since eq. ([F.4j) takes the same form as when g(r) = 0, we can use this equation 
to define a superpotential just as for T = 0: the second order equation (|F.4j) can be 
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replaced by the two first order equations, and the system becomes: 

4 dW ($) 

A = -~w($), # = JlU., (f.7) 

9 " +9 ' 2 '4wm (f.8) 



9' 9 

2 



4 /diy($)\ 64 o 4 
-3 (^^J + ^ ($) 2 - -W (*) 9' = e"^($). (F.9) 

Equations ([F.7|) provide the definition of the thermal superpotential W(Q). 



F.2 Counting integration constants: uniqueness properties of BH solu- 
tions 

The system ( F.7| - [F.9| ) has to be solved for the functions A, §,g,W once V($) is 
given as input. As in the T = case, once is given the scale factor and dilaton 

are uniquely fixed up to a single physical integration constant (a choice of scale) 
However, now is not simply a solution of a first order equation, as in ( |2.7| ). 

Since we cannot decouple the equation for W, we cannot choose it as an input as we 
did in the T = case. 

As it turns out, the most natural place where to fix the integration constants is 
the horizon, rather than the boundary. In other words, the general solution of the 
fifth order system (|F.7[ - [F.9| ) is most easily parametrized by the horizon value of the 
functions involved. 

Consider a black- hole with a regular horizon located at u = Uh- As we approach 
the horizon, u — > Uh, the dilaton, the scale factor and the superpotential have regular 
expansions, 

A -> A h - (u h - u)A'{u h ) + . . . , A -> \ h - {u h - u)\'{u h ) + ..., 

W->W h -(u h - u)\'(u h )d x W(\ h ), (F.10) 

where Wh = W(Xh)- On the other hand close to the horizon we must have have: 

g = \og(u h -u)+g h + ((u h - uf) , u -> u h (F.ll) 

where gh = log[— f'(uh)]- Substituting these values in equation ( [b'.9[ ) we see that 
regularity at the horizon requires the condition: 

W h = ~e-^(A A ). (F.12) 

On the other hand, the quantities Ah, \h and g^ are free. The differential equa- 
tion determines the following terms in the expansion around in terms of these 
quantities, and of V(A) and its derivatives at A^. 

We can use the horizon quantities to fix all integration integration constants of 
the system (|K.7| - |F.9|) . The system is first order in A(u),<&(u) and W(u), and second 



-87- 



order in g(u), therefore it contains five integration constants. One of them is un- 
physical, and it is due to reparametrization invariance. It can be eliminated by using 
A in place of u as a coordinate, or it can be fixed by setting Uh to an arbitrary value. 
The four quantities \h, A h , g h ,Wh provide the remaining four integration constants. 
For an arbitrary choice the solution will be singular at the horizon, whereas regularity 
is assured by imposing the constraint ( |F.12j ). 

Notice that the value of the potential at the horizon has to be positive in order 
to get a well-behaved black-hole solution. This is similar to Gubser's criterion iden- 



tifying good singularities |29|| 



What we have shown above means that the theory has a three-parameter family 
of regular black-hole solutions, characterized by the three real numbers Xh,Ah,gt- 
However, since the initial conditions were set at the horizon, these solutions will 
not all have the same UV asymptotics. To understand what happens in the UV, 
u — > -co in these coordinates), notice that close to Uh we have W{u) > 0, g'(u) < 0. 
Moreover by eq. ( |F.8| ), g"(u) < as long as g'(u) < and W > 0. Thus, 
decreases as we go further away from the horizon. As a consequence, the extra term 
—Wg' in (|F.9|) , that did not appear in the zero-temperature equation (|2.7| ), become 
less and less important as we move away from the horizon towards the asymptotic 
region. At the same time due to the relative signs of g' and g", g(u) approaches a 
constant value g , and the r.h.s. of eq. QF.9|) approaches the r.h.s. of (|2.7|) up to a 
multiplicative constant. Therefore the solution will get closer and closer to one of 
the zero-temperature solutions, up to a rescaling of W(&): W(X) — > e g °/ 2 W (X) as 
A — > 0. The existence of a solution that connects the horizon to the UV boundary 



will be proved more rigorously in Appendix |H.4 . 

This is not the end of the story, since we want the the black-hole solution to have 
the same UV asymptotics as the zero-temperature solutions. Notice however that 



eqs. (|F.7| - |F.9|) are invariant under the two following independent transformations: 

X(u), A(u),g(u), W{u) -> X(ue- Sl ), A(ue- Sl ), g(ue~ Sl ) - 25 u e[W '(ite^F .13) 
X(u), A(u),g{u), W(u) - A(V), A(u) + 5 2 , g(u), W(u) (F.14) 

where 8\ t 2 are arbitrary real numbers. These transformations map a solution into 
another solution, and also preserve the regularity condition (|F.12|) . Therefore one 
can use these transformations to move in the space of solutions, and reach the one 
with the desired asymptotics. Specifically, one can construct a regular UV black-hole 
solution with UV asymptotics matching a given T = background in the following 
way: 

1. choose an arbitrary horizon position Uh and fix arbitrary values for the initial 
data at the horizon, namely Ah,gh, A^ 

2. fix the fourth initial data, Wh, according to the regularity condition. 
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3. Evolve the solution from the horizon to the UV. In general, as u — > — oo, g(u) 
will go to a constant guv 7^ 0. 

4. Use a symmetry transformation with parameter 5± = gxjyjl. In the new solu- 
tion g goes to the correct UV limit (namely, zero). 

5. Use a 5 2 transformation to reset the overall scale of the solution to the desired 
value (e.g. to match a given T = solution). (This does not affect g(u)) 

At the end of this procedure, the only free parameter remains the initial choice of Xh- 
Thus, for each A^, the solution with given UV boundary conditions is unique. For 
a given choice of UV asymptotics, the black-hole metric and temperature depends 
only on X h , which can then be used as an unambiguous quantity to parametrize the 
different solutions. 

F.3 Asymptotic form of the solution 

Here we determine the behavior of the solution of the finite- T generalization of the 
superpotential equations, ( |F.9| ), in the small A and large-A region, respectively. 

The nontrivial part of eqs.( [F.7|4F.9l ) is the one that determines W and g; then, 



A and $ follow as in the zero-temperature case. We can decouple the W — g system 
from the rest, as follows: First using ( [F.7|) , we rewrite eqs. (|F.8|) and ([F.9|) keeping 
<3> as the independent variable. In this way, the four equations ( [F.7HF.9] ) split in two 
independent sets of equations. For the scale factor we have (as for zero-temperature) 



4 W 

^ = -95^ (F ' 15) 



for PU($) and g($) we find: 



d *9 + JT^) d * W + d l W = ^ W (F.16) 
d^gj 9 

4 4 o 64 

--W{d*W){d*g) - - {d^Wf + -W 2 = e-«V{\) (F.17) 



Eq. (|F.16|) can be integrated to a closed expression for g{<&) in terms of W(&). 
Write 

16 w _ Q2 rxr 

9 $ [g + log(-^)] = 9 Q ^ = F w (*) (F.18) 

Thus, 

g + ]og(-d*g) = J d$F w ($), (F.19) 

and exponentiating: 

(-d*g)e 9 = ef d * Fw ^. (F.20) 
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Integrating one more time we obtain an explicit expression for g(<&) in terms of the 
(still unknown) function F\y{&)'- 



g($) = log J d$eS d *' Fw{ -*'\ (F.21) 



F.3.1 Solution of the W($)-g{$>) system in the A — >• limit 

As usual, we assume a power series expansion of V(X) as A — > 0. We take a power 
series ansatz for W(A) as well, as for its zero temperature counterpart, 

W(X) = W + W X X + W 2 X 2 + ... (F.22) 

where the expansion coefficients W\ are a priori temperature dependent. 

Using this ansatz in eq. (F.21 ), written in terms of A = e*, we can obtain the 



form of g(X) for small A. We have from eq ( F.18Q 



9Wx X 9 I Wi 



F W {X) = -^ T + - ( 7-32^^ ) +0(A), A^O. (F.23) 



Using this expression in eq. ( F.21|) we obtain: 



g(X) = log 9l + g 2 X^e ^ * [1 + 0(A)] 7 = - 32^^ (F.24) 

where g\ and ^2 are integration constants. To recover the UV boundary condition 
g(X) — > as A — » 0, we must fix = 1. (? 2 ultimately determines the temperature 
of the solution, and is only required to be negative. 

Next, we insert the asyptotics (|F.24Q in the equation ( |b\17| ), in order to determine 
the coefficients. For small A it reads: 



A -16W 1 A a A 

■-g2W{d x W)X' +2 e 5^ -[i + (X)]--X 2 (d x W) 2 +^W 2 = V(X) 



-16Wq 1 

1 - g 2 X y e 9W i x + 



(F.25) 

To any finite order in powers of A, this equation is the same as the zero-temperature 
superpotential equation, ( |E.1| ). It follows that the power series expansion of W(X) 
is not affected by the temperature, and all the coefficients are the same as in the 
zero temperature solution: Wi — Wi. The difference between the finite-T and zero-T 
solutions are of order A 7 e~ 4 ^'' oA - ) ~ r 4 , and imply a temperature dependent value for 
the gluon condensate, cfr. Appendix E. 

Since the series coefficients of W(X) completely determine the UV series in inverse 
logarithms of r of the metric and dilaton, it follows that such series has the same 
form for any temperature. 
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F.3.2 Solution of the W($)-g($) system in the asymptotic large $ region 

Next we want to solve the system of eqs ( |5.12| - |ihT3| ), for in the asymptotic region of 
large $. This is defined as the region beyond some $o where the potential can be 
well approximated by its leading asymptotic 29 , 



V($) ~ e 2Q$ $ p $ > $ c 



(F.26) 



We assume the horizon is situated in this region, i.e. we work in the limit 

$o- The IR asymptotics of the zero temperature superpotential W ($) giving rise to 

the well behaved solutions with the repulsive singularity are of the form: 



IF ~$ p / 2 e Q<I> ,$>$ . 
Let us try the ansatz for the large <3? behavior of W: 



(F.27) 



(F.28) 



With this ansatz we can directly calculate the r.h.s. of eq. ( F.18| ). To leading and 
first sub leading order we have: 



R 



16 



9Q 



2 V 9Q 2 



(F.29) 



and from eq. ( |F.20| ) we obtain: 

g^)=\og{c 2 -C^ R e Kq> 



1 + 



;f.3o) 



where C\ and C 2 are two integration constants. This solution is supposed to be valid 
in the whole region $o $ < 3>h, and we can relate $/i to the integration constants 
by going to the horizon, g — > — oo: 



C 2 -Ci(^e 







(F.31) 



so we can write eq. ( F.30| ) as: 
g(#) = log cJl 



R 



o 



(F.32) 



Finally, requiring g ~ for 1 <C $ <C fixes Ci — 1. 
The appropriate solution for / = e 9 is therefore: 



/(*) = ! 



(F.33) 



29 the actual value of $o is immaterial, and will depend on the specific choice for the potential. 
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which correctly interpolates between the desired behavior at $ <C $^ and $ = 
This solution is valid in the whole region $ > $o up to the horizon $ = $h 30 . 

Now let us look at equation ( |5.13| ). Using ( |F.33| ), and neglecting 0(l/<3> 2 ) terms, 
it becomes: 

®\ R e K(<i>-<i> h ) 



i - ^ 

For $ -C $/i this requires: 

Q = Q, P = P, (F.35) 

i.e. the superpotential must have the same large $ asymptotics as the zero-temperature 
special solution. Then, eq. ( {Fl 3 4 ) is equivalent to the ^-independent algebraic equa- 
tions: 

A A 64 f PK \ 

--QK = --Q 2 + -, ^-— + RQj=-Q P (F.36) 

which are identically satisfied, due to the definitions of K and R in ( |F.29p ! 

Therefore the asymptotic solution with horizon at $ = $^ is , to this order: 

W($) W ($) ~ e Q *$ p/2 (F.37) 



/(<&)- 1- ( — 1 exp[-^($,-$)], (F.38) 

valid in the whole asymptotic region $^ > $ 3> $o- 

In the particular case of power-law behavior, A (r) ~ —Cr a , corresponding to 
Q = 2/3, P = (a — l)/a (see p8|), the solution in r-coordinates is, for large r : 

A{r) ~ -CV a , (F.39) 
$(r) ~ ^r a + ^(a- l)logr, (F.40) 

/ ~ $-f P e 2,I> ~ r -|(o-l) e 3Cr««+|(«-l)Iogr _ ^-c^SCr- ^ ^ 

The horizon position is obtained by inverting $(7^) = 

We can verify directly that these asymptotics solve eqs. (p3|- p.6| ), for V ~ e3*. 
Using the asymptotic form b(r) = b e~ CrC " , we get, integrating eq (|3.7[) 

/r p3Cr a 
dtb~ 3 (t) = C 2 + ^ / dss-^e' 

^ C 2 + -^r^e 3 ^ (F.42) 
3aC 



:!0 



For $ > <f>/j the solution is the same, but the definition of g changes: / 
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II 



Figure 16: Two possible types of cusps at which the small and big BH in pair can merge. 
These are denoted by "S" and "B" , respectively. The full curve F(T) should be constructed 
by various combinations of these vertices. 

Substituting into the equation for the potential yields 



V(r) 





!'[ 


-\ 

b 2 






-\ 


b 2 1 


9a 2 C 2 c 


bon 2 



b 2 b 
2 ¥ + l 



ib 



To l I *-^2 



1 _ r l-a e 3Cr a 



3aC 



(3a 2 C 4 r 2 ^) + Cie 3Cra (-aCr^ 1 ) J 



C 2(a-l) e 2r° 



(F.43) 



This is the expected leading behavior for the potential. The contributions pro- 
portional to C±e ~ CieT* cancel between the two terms containing / and / 
respectively. V is independent of C\, which encodes the horizon position. Further- 
more, comparing eq. (|F.43|) with the expression for V(r) one obtains with / = 1, 
we see that Ci must be set to = 1 if the potential has to be the same as in the 
zero-temperature solution, in agreement with eq. ( |5.19|) . 

G. Multiple Big Black-holes 



In this section we generalize the proposition of sec. |5.3| to a much larger class of 
geometries, for which the functions F(Xh) and T(A^) may acquire multiple extrema 
31 , see figs. [18] and [19] for an example. In these cases there are more than just two 
BHs. However for the generic confining theories 32 , they still come in pairs of one 
small (T'(Aft) > 0) and one big (T'(\ h ) < 0) BH, connected at an extremum of 
T(Xh). Therefore there is an even number of BHs in total. This follows from the fact 
that there always exists an asymptotically AdS big BH for A^ — > and a small BH 
for X h — > oo and one extremum of T(X h ) creates two branches (a small and a big), 
by definition. 



31 We assume that these functions are C°° in A/, . 
32 except the borderline case of a = 1. 
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In particular we want to prove the following: 



1. Existence of a deconfinement phase transition, 

2. Finite latent heat (first order transition), 

3. Continuity of F function of T, 

4. Uniqueness of the deconfinement transition. 

We shall first present a graphical proof of these points. For illustration purposes 
we make the following additional assumption (that seems natural and satisfied by a 
large class of potentials): C v = TdS/dT is negative (positive) for small (big) BHs 33 . 
After the graphical demonstration below, we provide an analytic proof which applies 
to cases when this additional assumption is weakened. The assumption about the 
specific heat is actually not needed for point 1 in the above list, so that the main 
result (the existence of a phase transition) is valid for an arbitrary behavior of the 
specific heat. The function F(T) is determined from F(Xh) and A^(T). Since the 
latter is generally multi- valued, so is F(T). As a multivalued function it can have a 
very complicated form with cusps and crossings, see fig. Ofor an example. Although 



complicated, the form of F(T) is restricted by certain rules: 



1. On every piece of the F(T) curve, F'(T) < 0. This follows from the positivity 
of entropy. 

2. On the small black- hole branches F"(T) > 0, and on the big black-hole branches 
F"{T) < 0. This follows from our assumption above and from F" oc —C v . 

3. There should always be a big black-hole branch (which asymptotically be- 
comes the AdS black-hole at high-T) in the high-T (small Xh) region, on which 
F(T) — ► — oo as T — > oo. 

4. There should always be a small black-hole branch in the high-T (large ) 
region, on which F(T) — > 0. This follows from the discussion in section |5.2| . 

5. The small and the big BHs always come in pairs, hence there are equal numbers 
of branches on the F(T) curve, with negative and positive F"(T). This is clear 
from the fact that one small and one big black-hole branches off at an extremum 
of T(\ h ). 

6. These merging points of a pair of big and small BHs are represented by a cusp. 
There are two possible types of cusps as shown in fig. [16]. These particular 
shapes follow from the first two properties above. Since at the merging points 
the entropy is the same, T' is the same on the two branches of the cusp. 
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Fb 




Figure 17: Example of a curve F(T) that exhibits multiple extrema. S\ and S2 denote 
small BHs whereas B\ and B2 denote big BHs. The arrows represent direction of increasing 
Xh- 




Figure 18: T as a function of in the example of fig. 17 



For an example of a curve F(T) with these properties, see fig. |T7[ Given these 
properties, it is not hard to show that the minimum energy configuration for T > T c 



as in the case of the AdS^ BH with spherical horizon. 
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F 




31 



Figure 19: F as a function of in the example of fig. O. 



is always a big BH. This is because, the entire curve F(T) should be formed out 



of the vertices given in fig. [T|] connected with small and big BH legs. Clearly for 
any small BH, there exist a big BH that stems from the same vertex which has lower 
energy. Therefore, in the entire curve for F < 0, the lowest energy configuration 
should be a big black- hole. 

Of course, below T = T c the lowest energy configuration becomes the thermal 
gas and the curve F min (T) for the minimum energy configuration always looks like 
in fig. [7[ For the example given in fig [17|, the corresponding free energy diagram is 
constructed in|^. The reason that F(T) should cross zero (point 1 in the list above), 
hence exhibit a phase transition, follows from the 2nd and the 3rd properties above. 
Similarly, the points 2, 3 and 4 simply follow from the fact that one cannot draw a 
function F(T) that violates these points with the properties listed above. Thus, we 
demonstrated what we wanted: the proposition of section |5.5| apply to the multiple 
extrema cases as well. 

Before going into the analytic proof of these statements, let us analyze the ex- 



ample of fig. |T7| in more detail. In this particular case, there are two small and big 
BH pairs that are denoted by (S\,Bi) and (S2,B 2 ). The black-holes in each pair 
merge at two local minima of T(X h ) (that are denoted by T^ in and T^ in )) and the 



two pairs are connected at a local maxima of T(Xh) (denoted by T max ). See fig. [18 



which plots T(Xh), that corresponds to fig. [17| for a clear demonstration of these 



34 There may be more than one T c on which F vanishes and our statement applies to all of these 
points. 
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facts. We also present the function F(Xh) that corresponds to this example in fig. 



[19]. Given fig. [18] and fig. [T^, fig. [L7| follows by solving for F as a function of T 
parametrically. The arrows in fig. |I7] point to the direction of increasing A^. 

Fig. |l7| exhibits various first order transitions between different branches. First 
of all there are transitions between thermal gas (with F(T) = 0) and B\, B 2 and 
Si. We denote these points as Tf, T c and T s respectively in fig. [17|. However, 



not all of these points correspond to actual phase transitions. In order to obtain 
the true free energy of the system, one should draw F(T) on the minimum energy 



configuration. This looks is much simpler and the one corresponding to fig. |I7] is 
given by fig. |[ In particular, we see that out of Tf, T c and T s , only T c is a real 
phase transition. It is in fact the confinement- deconfinement phase transition in 
this example. Tf corresponds to a "fake" deconfinement transition, because B\ has 
higher energy than B2 at the point it crosses F — 0. Similarly, at T s , there is a fake 
transition between a small BH and the thermal gas geometry. Although these are 
fake transitions, hence uninteresting for the dual field theory point of view, they may 
bear some interest on the bulk as they describe possible transitions between various 
different geometries in asymptotically AdS spaces. The fake transitions described 



here parallels the transitions found in R 4 and F 4 corrected AdS geometry in [[5TJ. In 



addition, there are various transitions among small and big BHs (S2 and B\ in the 
figure) and different small black-holes [S\ and S 2 in the figure). 

However there is another interesting possibility that is also present in this exam- 
ple: a first order transition between two different big black-holes atT b . This transition 
is not "fake" as others, because it corresponds to first order transition between two 
minimum energy configurations, B\ and B2 in fig [17]. See also fig. ^| It will be 
interesting to investigate whether or not a dilaton potential V(X) with this property 
exists; if so, whether or not this corresponds to a meaningful phase transition in the 
dual gauge theory. 

G.l Analytic demonstration 

Suppose the function T(A^) has an arbitrary number n of minima, with n (even) odd 
for (non-)confining IR asymptotics 35 , corresponding to certain values Aj, i = 1 . . . n. 
Then, there will be n + 1 black-hole branches, each corresponding to the ranges 
(Aj, Aj+i), with Ao = and A n+ i = +00. In each branch the free energy as a function 
of T is single valued, and is given by the first law, Fi(T) = — J Si(T')dT' + Ci where 
Si is the entropy function in the ith branch (i.e. in the interval (Aj_i,Aj)), and C* 
are integration constants. We can still write the free energy in compact form as 

F{X h ) = ^ dX'S(X')^, (G.l) 



.",5 



except in the borderline confining case where there is also an even number of minima 
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as this is the unique continuous function that satisfies dT jdT = —S on every branch, 
and vanishes as A^ — > oo. Since S is nowhere vanishing for finite A, a minimum for 
T(A) corresponds to a maximum for ^(A). Using the integral expression of the 



free energy, we can give an analytic proof of the proposition in Section ^3 relating 
confinement with the existence of phase transitions. 

Let us consider confining asymptotics. We want to show that these always exhibit 
a first order phase transition. In particular we want to prove the following three 
statements: 

1. JF changes sign at some finite A. Specifically, T ^ — oo as A ^ 0, whereas 
T{K) > 0; 

2. For every A« corresponding to a minimum of T(A) , either is positive or it 
is larger than the free energy of some big black-hole with the same temperature; 

3. For every \ corresponding to a maximum of T(A), ^(Aj) there exists a big 
black-hole with the same temperature and lower free energy. 

Step 1 by itself shows that there is indeed a phase transition. It goes in the right 
direction, since at small temperature (T < T min ) there are no black-hole solutions, 
and moreover in any black-hole branch the free energy is a decreasing function of 
temperature. Therefore the black-hole free energy crosses from positive to negative 
(with respect to the vacuum) as the temperature increases. 

Steps 2 and 3 imply that the free energy never jumps in either directions: the 
branching points when a large/small black-hole pair appears or disappears have al- 
ways higher free energy than some other state in the ensemble (which therefore must 
have been dominant starting from some lower temperature). If this weren't true, it 
would be possible to have discontinuous jumps in the free energy, i.e. phase tran- 
sitions with infinite latent heat. In this step, one needs an extra assumption about 
the ordering of the black-hole branches. Step 3 also implies that ones the system is 
in the black-hole phase for some temperature, it stays in a black-hole phase for all 
higher temperatures, so there cannot be any inverse "reconfining" transition. 

The proof of Step 1 is straightforward, and identical to the argument we used in 
section 5.3 in the case of n = 1, i.e. when T(A^) has a single minimum. Thus, the 
existence of a deconfining phase transition per se is easy to prove. To prove Steps 
2 and 3 however, in the general case n > 1 we need to make an extra assumption 
about the nature of the entropy S(X): 

Assumption: 

the black-hole branches are "ordered" in the sense that Si(X) > 5j + i(A / ) for any 
A e (Xi-i, A;) and any A' G (A;, A m ). 
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This still allows some local violation of the monotonicity of S(X) within each branch. 
Although we cannot exclude violation of this assumption from first principles, it is 
satisfied in all examples we have studied numerically, where we also found cases with 
strict monotonicity of S(X) violated (these cases correspond to regions of stable small 
black-holes, as discussed in section 5.2.2). The above assumption might be relaxed, 
since it is probably not a necessary condition for the statements 2 and 3 to hold, but 
it allows the construction of a simple enough proof. 

With the above assumption, the proof of the proposition goes as follows. 

Step 1 

The fact that T < for small enough A follows from eq. ( |5.6| ), whose validity does 
not depend on the number of extrema of T(Xh)\ on the other hand, evaluating the 
free energy on the extremum (a minimum) of T(Xh) with larger A^, we obtain: 

r+°° at 

T{X n ) = S—>0, (G.2) 

since in the last branch dT/dX > 0. This is true without additional assumptions 
about the function 5(A). 

Step 2 

First, consider the extremum Aj corresponding to the absolute minimum of T(A/ l ), 
T(Xi) = T m i n . Since the last extremum X n is also a minimum, there will be 2k = n — i 
minima and 2k — 1 maxima in the region A > Aj, corresponding to 2k maxima and 
2k — 1 minima for J-{X). We are going to show that ^(Aj) — J^(A n ) > 0. Let us 
evaluate the black- hole free energy at the point Aj: 

r+°° at 

HV = l six)- 

r*i+i at f Xl+2 dT f Xn dT 

Each integral is extended over a different black-hole branch. In every branch, we can 
use the mean value theorem: 

J x S(X')— = S J+1 (T J+1 -T J ), (G.4) 

where Sj + i = S(Xj + i), for some appropriate value A J+ i with Aj < A, + i < Xj + i. 
Then, ( |G.3|) becomes: 

F\Xi) ~ FyXn) = Si + \{T i+ i — Tj) + Si + 2(T i+2 — T i+ i) + . . . S n (T n — T n _i). (G.5) 



- 99 - 



The sum has alternating sign since, for any I > 0, Tj + 2/ are local minima and Tj + 2z+i 
are local maxima. By assumption, Si+i > Si + 2 > ■ ■ ■ S n , therefore: 

J-'(Xi) — T\\n) > Si + \(Ti + i — Ti) + Si + i(Ti + 2 — Ti + i) + . . . S n -i(T n — T n __i) 
= Si + i(Ti + 2 — Ti) + Si + 3(Ti + 4 — Ti + 2) + . . . S n -i(T n — T n _i) (G.6) 

Notice that now only temperature corresponding to local minima appear. Next, 
subtract the combination S n ~i(T n — Ti) from both sides of the above inequality: 

^"(Af) — T{\ n ) — S n -\\T n — Ti) > Si + \{Ti + 2 — Ti) + £1-1-3 (Tj_t_4 — Ti + 2) + . . . 

. . . + S n -s(T n _2 — Tn_4j + S , n _i(T n — T n _2) — S n -i(T n — Ti) 

= Si + \(Ti + 2 — Ti) + Si + s(Ti + 4 — Ti + 2) + . . . + 5'„,-3(T n _2 — T n ^) + S n -i(Ti — T n _2) 

> ^+1(^+2 — Ti) + Si + s(T i+4: — T i+2 ) + . . . + S n s(T n _2 — T n _ 4 ) + S n _3(Ti — T n _ 2 ) 
— Si + \(T i+ 2 — Ti) + Si + s(Ti +4: — T i+2 ) + . . . + S'„,_3(Tj — T n _ 4 ) 

> . . . > S l+1 (T i+2 - Ti) + S i+3 (Ti - T i+2 ) = (S i+1 - S i+3 )(T i+2 - T t ) > (G.7) 

In each subsequent step we have used Si + i > Si +3 and (TJ +2 — Ti) < since we took 
Ti to be the absolute minimum. Therefore, from the first and last side of (|G.7|) we 
get: 

^(Ai) > F{\ n ) + 5 n _x(T n - Ti) > (G.8) 

This shows that at the minimum temperature when a black-hole pair appears, the 
free energy is positive, and the thermal gas background still dominates the ensemble, 
so the global free energy of the system does not jump abruptly. 

To show that the free energy does not exhibits jump at the creation of subsequent 
black-hole pairs occurring at T > T min , we proceed as follows. Consider the case 
when Aj is a generic local minimum. If all the subsequent minima Aj + 2 • • • A n have 
higher temperature, then we can proceed exactly as above and show that F(\) > 0; 
Otherwise, if there is a local minimum for some Aj + 2z > Aj, with T i+2 i < Ti, then there 
is also a big black-hole with temperature Tb = Ti, corresponding to a point A# € 
(K+21-1, Aj+2j)- It is easy to show that J-(Xi) > JF(As), using the same procedure 
that led to flGTTD : 

^(Ai) - F{\ B ) = [ B d\' S(\')T\\') 

— Si + i(T i+ i — Ti) + Si + 2(T i+ 2 — Ti + i) + . . . Sb(Tb — Ti+21-1) 
> Si + \(Ti + 2 —Ti) + Si +3 (Ti + 4 — Ti + 2) + . . . Sb{Tb — Ti+21-1) 
>S B (T B -T l ) = (G.9) 

Therefore, the black-hole pair that appears at Tj cannot dominate the ensemble right 
at Ti, and the global free energy of the system does not jump. 
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Step 3 



Finally, we show that there cannot be an increase of the free energy back above 
zero, for temperatures higher than the critical temperature. Since in each branch 
the free energy of single black-holes is monotonically decreasing, the only way the 
global free energy can increase is by jumping up at a point where a small/big black- 
hole pair disappears, i.e. at a local maximum of T(A^) (i.e. a local minimum of 
J-{\h) )• Therefore, it is sufficient to show that for any maximum \ of T(X h ), with 
temperature T i: there exist a black hole with the same temperature T B = T { but 
lower free energy, so that the system does is not forced to jump at T — %. Since 
T(A — > 0) = +00, there certainly exists at least one big black-hole with T B = % and 
\ B < Aj. Let us consider the closest one to Aj. Aj will be separated from \ B by an 
odd number of extrema (/ minima and I — 1 maxima). Proceeding along similar lines 
as in Step 2, we now compute the difference T{\ B ) — J-{\i) and show it is negative: 



In the last line, we replaced T B with Tj in the first term. Now only temperatures of 
local maxima appear, of which Tj is the highest one. Using repeatedly the inequalities 

Si-2k-i > Si-2k+i > ■■■ Si-i we have: 



Thus, when a black- hole pair disappears at temperature Tj, the system is already on 
another big black-hole branch with lower free energy. 

So far we have treated the case of confining asymptotics. The converse can also 
be proven under the same assumptions: non-confining asymptotics do not exhibit a 
thermal gas/black-hole phase transition 36 , rather the system is in a (big) black- hole 
phase for any T > 0. The proof proceeds along the same lines as in the confining 
case. 

36 0ne cannot talk about a deconfining transition here, since the zero-temperature solution itself 
is not confining 




— >S'i-2fc-l(Xj-2fc-l — T B ) + Si-2k(Ti_2k ~ ^i-2fe-l) + • • • + ^(Tj — Tj_i) 



< Si_2k-l{Ti-2k — T B ) + Si_2k+l{Ti~2k+2 — T%-2k) • • • + <S'j_i(T'j — Tj_ 2 ) 
Si-2k-l{Ti-2k — Ti) + S'j_2fc+l(7i--2fe+2 — T%-1k) • • • + Si-\\Ti — Tj_2) 



(G.10) 
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H. Details of the computations with scalar variables 



H.l Equivalence to Einstein's equations 



Here we prove that the reduced system of equations presented in the section |7| are 
equivalent to the equations of motion in the u- variable in (|A.17| - |A~20"1) . For this 
purpose we should supplement (|7.2|) and (|7.3| ) by the equations that determine the 
metric functions. These are given by the following first order equations: 



A! 
9 



1 4 (•* v 

t 

3X 4 /■* 



AY 



J — oo 



'il 4 ri' Y 
e _ 3 J-oo A . 



(H.l) 
(H.2) 
(H.3) 



Using d/du = $'(d/d$), (fT2|) , ( |PD and the three equations above in ( |A.17HA.20| ), 
it is straightforward to show that they are all solved. 

To convert the system (7.2) , (|7.3| ) , ( |H.l| - |[i.3D to the conformal coordinate system, 
one uses, 



du 
dr 



(H.4) 



Now, we use this equation to show how r and A related near the boundary. Converting 
( [ti.2[ ) into the r- variable by ( |H.4| ), changing variable to A = exp($) and using the 
equation f [7.4[) , one has: 



dX _3X 

~^ ~ T 



(H.5) 



We use the expansion of X near the boundary in ( |2.12|) to get, 
d(b \) 1 



(6oA) 2 e Au "^(6oAo)~ b e^(6 A) b (1 + 0(A)) . 



dr 



(H.6) 



Now, we use the definition of the QCD scale A in ( |2.19| ) and integrate ( |H.6| ) to obtain: 

r A = e'^(b X)- b . (H.7) 

The corresponding relation involving the domain wall coordinates is obtained by 
integrating eq. (|H.2|) . The result in terms of X(u) reads: 



1 /4 \ u 
— +(- + 6)log(6 A) = - 7 -log A/ 



(H. 
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H.2 Solution of eq. (7.3) 

Here, we note that Q7.3Q can be solved in terms of X explicitly. Define, 

A(X 2 -1) 4 
C($) = 3X ' m = ~3X- (H ' 9) 

Then, the solution is, 

Y ($) = ef* <^ d *' ^ - J* d$ / d($ / )e / *' c(l>)dl> ) • (H.10) 

This is the general solution of (|7.3| ) for non-zero Y. As already mentioned, Y = is 
also a consistent solution, which corresponds to the thermal gas. 

H.3 A fixed point analysis of the X-Y system 

In order to understand the number of integration constants in the system, one can 
perform a fixed point analysis of the XY system given by eqs. (|7.3|) and (|7.2|) . 

It is obvious that Y = is a fixed line in the phase space. This line corresponds 
to the zero temperature solutions discussed in [p8fl . We focus on this case first. 
The solution is determined by (|7.2|) at Y = 0. Clearly, there are four fixed points 
of the system for an arbitrary potential. These are given by X = ±1 and ±oo. 
For the special class of exponential potentials V ~ exp(a<£>), there is the additional 
fixed point X = — 3a/8. Furthermore, X = is also a fixed point for the types of 
potentials we study in this paper. 

We solve the system by specifying the boundary conditions (or the asymptotic 
behavior) in the IR, for large A. Suppose X = Xf at A = A/. Let us assume for 
the moment that Xf < (as is the case for the "special" and "generic" solutions 
discussed in this paper). Then, one obtains the following: 

1. If — 1 < Xf < 0, then X —>■ as A decreases, at A = 0. 

2. If — oo < X < — 1 then X — ► — oo as A decreases at A = Aj ^ 0. 

To prove the above statements, let us first show that X = —1 is a repulsive fixed 
point in the direction of decreasing A. To show this, we substitute X = — 1 + e, 
1 > e > in Q. One has, 

with the solution, 

e = c \/($)e- 8<I>/3 . (H.12) 

As $ — ► oo V — ► e 4 */ 3 in our case. Hence one falls into X = — 1 as $ - >oo. Note 
that the fixed point would instead be attractive for decreasing A, were V e a */ 3 
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Figure 20: All solutions to the X equation of motion are shown. The thick (red) curve 
corresponds to our solution Xq, that flows to the fixed point X = —1/2. One clearly 
observes the fixed points X = ±l,±oo and 1/2 in the figure. The direction of flow (as a 
function of r) is towards the left for X > and towards the right for X < 0. 



with a > 8/3. The analysis is the same for X = — 1 — e, 1 ^> e > 0, as one obtains 
the same equation, (|H.11| ). 

Now let us focus on the vicinity of X = 0, by writing X = — e, 1 ^> e > 0. One 
obtains 

de 4 1 dlog(V) 

M = 3~Ye^^- (R13) 
Let us assume that one can reach X = at finite $ = Then, one can ignore the 



first term in the RHS of ( H.13 ) above and obtain the follwing solution: 

e 2 = -log(y($)/y($ / )) (H.14) 

as $ — > <&f. This shows that X = can never be reached in the decreasing A 
direction, in finite A-time. Instead X always runs into X = 0, in the direction of 
decreasing A. 

On the other, hand QH.14|) allows to pass the X = point in the direction of 



increasing A, in finite A-time. These solutions (as functions of r) continue to the 



positive X region and hit back to the X = 1 fixed point as r increases, see figure pO 
However, in these solutions the derivative of A(r) changes sign at the locus X = 0. 
As our purpose is to find solutions dual to field theories with negative definite (3- 
functions, these do not correspond to any reasonable theories. 

One can easily carry out a similar analysis in the vicinity of X = — oo. For 
this purpose, we define Z — 1/X, and focus on the vicinity of Z — > CT by defining 
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Figure 21: Flow chart of the solutions for Y = case. The arrows show the direction of 
decreasing A. 



Z = — e > 0. To leading order, 

de_ _ 4 _ e d\og{V) 
rf$ ~ 3 ~ 2 d$ 

with the solution, 



(H.15) 



4 /■*/ 

e=-V r [$]-a/ W (H.16) 

where $j finite, is an integration constant. It is clear that, e goes to only at a finite 
point $j. 

This completes the proof of the assertions above. A similar analysis can be 
carried out in the region of X > 0. We summarize the behavior of flows in the 



direction of decreasing A in fig. 21 



The class of solutions to the X equation of motion for Y — 0, are summarized 



in fig. 20. There are five classes: 



First of all there is our main solution, Xq that flows to —1/2 as I increases. We 
recall that the fixed point X = —1/2 exists only for the potentials that has an 



asymptotic form V — > exp(4$/3). This is shown as the thick (red) curve in ^0 



Secondly, there are the solutions — 1 < X < Xq that flow to —1. These are 
also asymptotically AdS in the UV but have differ in the IR from X . The 



solutions described in [43] is an example of this class. 



• The solutions Xq < X < +1. They also have an AdS fixed point in the UV. 
As a function of r the solutions continue in the region of positive X, hence 
positive /3-function. This is not acceptable for an RG-flow in a field theory. 
The reason that such a behavior can happen in a gravity dual is because A(r) 
is determined by solving a second order differential equation whereas in the 
field theory it solves the first-order Gell-Mann-Low equation. Therefore we 
discard these solutions as un-physical for holographic purposes. 

• Finally there are the solutions — oo < X < — 1 and oo > X > 1 that posses 
negative and positive /5-functions respectively. They do not exhibit an asymp- 
totically AdS fixed point, hence are not useful for our purposes here. 

Now, it turns out that the above behavior of solutions easily generalize the case 
of Y 0. For simplicity, we focus only the region of the phase space that we are 
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X 



Figure 22: Flow chart of the solutions for the general case for non-zero Y. The arrows 
show the direction of decreasing A. 

interested in, i.e. X < 0, Y > and Y > X 2 — 1. In this case, the fixed point 
X = — 1 is replaced by the fixed line X = — yl + Y. The flow chart in the phase 
space is shown in fig. |22|. The thick (red) curve represents the fixed line of the 
system 1 — X 2 + Y = 0. As described in section |7.3| , Y diverges at the horizon 
that corresponds to a X = Xh = —3/8V'(\h)/V(\h). Therefore each different curve 
corresponds to a different temperature. 

It may be interesting to carry out an analysis as a function of r and obtain a flow 



chart that generalizes fig. BO to the case of the full X-Y system. This would determine 



the entire family of solutions to ( |7.2j , |7.3| ), most of which would be un-physical. We 
leave this question to future work. 



H.4 Near-horizon continuously connects near-boundary 

We want to show that the BH solution that is given by the initial conditions for 
X and Y at the horizon i.e. ( 7.21 ), ( [7.221) continuously extends over the solution 



near the boundary, i.e. Y ~ 0, X rs X . As shown by eqs. ( |7.8| ) and ( |7.9|) , the 
solution of the system becomes Y ~ 0, X « Xq as $ — > — oo, independently of the 
initial conditions at the horizon. We shall prove that this asymptotic UV solution 
is continuously connected to the solution near the horizon. This is achieved by 
examining the equations is motion ( |7.2| ) and ( |7.3| ). 
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First, we see from ( |7.3|) that X = Xq, Y = cannot be reached at a finite point 
$ = $/ with — oo < < Suppose Y = e near Then, one obtains 

de _4(1-X 2 ) 
~ 3X C 

with the solution 

4 

f = Ce" 



2 i 



where X/ = X($j). Thus, e = can only be reached as $ — > — oo (recall that 
X < 0) in contradiction with our assumption. Therefore we learn that Y and X 
should be finite at an arbitrary mid-point $/. 

The only way the near UV (<& < $/) and the near horizon ($><&/) solutions 
are detached would be a divergence in the RHS of (|7.2| ) and/or ( [7.3|) at $ = $/ 
as $ decreases down from In this region, V is bounded, X is also bounded as 
— 1 < X < 1. From the exact solution for F in (|H.10|) one finds that Y is also 



bounded. Hence the only way one can get a divergence is as X — > at $*. Now, we 
write X = — e in ( |7.2|) and find that e satisfies, 

where e = 1 + Y(&f) > 0. The solution is, 

which cannot be satisfied for any $ > because we assumed V A ($) monotonically 
increasing. Thus we proved that, given the initial conditions for X and Y at the 
horizon, the solution flows down to $ = — oo and connects continuously to the UV 
region where Y — > and X — > X . 

Note however that X = can be reached at a finite $ / for $ approaching $ f from 
below. This is analogous to the "bouncing" solution studied in the previous appendix. 
At finite temperature, this possibility is automatically ruled out by imposing the 
regularity condition at the horizon. 

H.5 Near boundary behavior of SX and Y: 

In order to derive the asymptotic behavior of SX and Y, it proves useful to obtain 
differential equations satisfied by these perturbations. This is done by expanding 
( |7.2|) and (|7.3| ) to first order in the perturbations and making use of ( |7.6|) . One 

dY 41 -X n 2 , 

= -o^^ H.17) 



obtains: 



d\ 3 X 



= SX 

d\ 



3Xq Xq — 1 Xq dX 
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The solution is straightforward: 



4 j-X dX 1—Xq 

Y(X) = Y(X Q ) e~ lJx o~^T 
5X(A) - ' lAoj 



K(A ) = Y e -^(6 A ) 



-46 



-Co(A )J — + Co(Ao)X 
Co(A ) = Co e-W(5 A„)- 4t . 



(H.19) 
(H.20) 

(H.21) 



Here Y and Co are integration constants. By expanding the above equations one 



obtains ([7.17 ) 



H.6 Free energy in A: details 

We will compute the free energy in the A-coordinate frame. In this frame, the metric 
of the thermal gas and the black-hole are as follows: 



ds 2 TG = B 2 (X) dt 2 + dx 2 + 



dX 2 



A)(A) 



ds 2 BH = B 2 {\) dt 2 F{\) + dx 2 + 



2 \ ' 



d\ 2 



F(\)D(\yj ■ 

Here the various metric functions are defined as follows: 



(H.22) 
(H.23) 



i rX dx l 3 Xn( X) 

Sb(A) = Sb(Ao)c» / *oTX5, 



i r x dx i_ 

X 



1 rx d> 

B(X) = B(X )e^ h o- 
F(X) = «t/o 



D (X) = 
D(X) = - 



, . 4 rX dX v 

B (X)e~* J ° T x » 



£ X 
31(A) 



e x 



, „ 4 rX d\ v 

B(X)e~* j ° T x , 



4 fA a y 

g3 JO X Jf . 



(H.24) 

(H.25) 
(H.26) 



They are obtained directly from the the expressions for the metric functions defined 
in the text in terms of the radial variable r, viz. (|7.4J) , (|7T5| ) and ( [H.2j ). We call the 
metric functions in A with the capital letters to distinguish them from the analogous 
functions of r. The relations are explicitly given by the following formulae: 



B (X) = b (r (A)) , 73(A) = b (r(A)) , F(X) = f (r(A)) 



(H.27) 



where r and r are determined by, 



ro(A) 



dX 



o XD (X) 



r(A) 



dX 



XD(X) 



(H.28) 



The expressions above completely determine the map between the r-frame and the 
A-frame. 
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H.6.1 Einstein contribution 

We first compute the Einstein contribution to the free energy. This is generally given 
by the frame-independent expression, 

Se = \m z [ ^gV. (H.29) 

JM 

M. is the manifold with a boundary. We regulate the integral in the A-frame by 
placing a cut-off at Ao- Thus, using the metric functions defined above, one obtains 
the following expression in the lambda variable, for the thermal gas solution, 

2 f°° 

S T E G = -M z (3% / fi (A)V(A) J D (A)- 1 . (H.30) 

Here (3' and V$ are the circumference of the Euclidean time and the volume of the 
space-part. They are related to the analogous quantities in the black- hole geometry, 
by matching the two solutions on the cut-off: 



P'Bo(Xo) = (3B(X )y/F(X ), V>B (\ ) 3 = V 3 B(\ ) 3 . (H.31) 
Now, we use the expression for the potential p8| , 

V(\) = ^(l-XZ)Jo^)f, (H .32) 

in ( [H.30| ) and see that it can be rewritten in the following form: 

= 4M/?y< B „(A„)VH^ r^l^lel/if^. (H.33) 



To obtain this expression, we used ( |H.25p , ( |H.24| ) and ( [H.32| ). We see that the 



integrand is a total derivative, hence the integral can be carried out exactly. One 
has, 

rrn 2 lit , A 4 r X v d\ 4 M d\ 1 ~ X Q °° 

S T E G = --M 3 /3V 3 %(A )V3 /o x »i e 5^T^ . (H.34) 

^ A 

It is easy to see that the contribution from oo vanishes, as the expression above 
scales a A -2 for large A for our confining IR asymptotics. Thus one only has the 
contribution at A : 

S™ = ~M 3 f3'ViB (\ ) 4 e-lti° x °f. (H.35) 

This is our final expression for the Einstein term on the TG geometry. 

Let us now compute the analogous contribution on the BH geometry. From 
( fTggD we get, 

S» H = \m 3 (3V, 5(A) 5 V(A) J D(A)- 1 . (H.36) 
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Following the same steps as before, we substitute the expression for D(\), B(X) and 
V(X) from ([EL25D , flEL24[) and Q), and obtain, 



^ = -^MW 3 5(Ao) 4 e-t^(^$)f 



A 



A 



g 3 J A A 



Again, the integrand is a total derivative and one has, 



Sf H = --M*(3VzB{X Ye 



4 fA „ dA 4 fA dA l-X 2 + y 



A g3 JA a 



^H.37) 



(H.38) 



Ao 



This can be simplified further: using (|7.3|) , one realizes that the integrand in the 
exponent is a total derivative of log F(A). Thus, one has, 



qBH 



■-M 3 (3V 3 B(X ) 4 e 



4 f A v d\ 
"3 JO A X 



Y(X C 
Y(X h 



[R.39) 



But Y(Xh) = oo by regularity condition at the horizon (see section ([7. 3D ), hence we 
have the final expression for the Einstein contribution on the BH geometry: 

-M 3 (3V 3 B(X ) 4 e—^^ ( x ~*)t. 



qBH 

J E 



[RAO) 



The Einstein contribution to the free energy follows from the difference of ( |H.35| ) 
and QH.40|) . In order to match the two expressions we use the matching conditions 
( [H.31D . Finally we also use (|H.26|) to obtain: 



5S E = S B E H - S T E G = jM 3 PV 3 B(X ) 4 (e- 



4 r X 0(y Y \d'\ 4 r A Y sd\ 

3 JO { X -x)T—e~3 Jo ( - A °~2X' , T 



. (H.41) 



As Ao is very small and the integrands in the expression above are very small in that 
region, one can expand the exponentials and obtain, 



5S E = -M 4 f3V 3 B(X c 



X 



25X 



[RA2) 



The functions Y and 5X are given in section (\7.2\) . Using these expressions it is 
straightforward to carry out the integrals. One obtains, 

2 



SS E (Xo) = --M'pVsCoBiXoTe-^iboXo) 



(H.43) 



Finally we remove the cut-off by sending Ao to and using the definition of the 
integration constant A: 



lim B(X )e~^(b X )- b = A i. 

Ao^O 



(H.44) 



Thus the final expression for the Einstein contribution to the free energy reads 

2 



;C M d V 3 (A 



H.45) 



- 110 - 



H.6.2 Gibbons-Hawking contribution 

We move on to the Gibbons-Hawking term that is given by the frame-independent 
expression, the second term in ( |2 . 1| ) . 

We first define two more functions of A in addition to ( |H.25| ): 

D A0 (X) = -Is (A)e-f^f^o(A) ) Da{x) = _l 5(A)e -!/o A f *<*). (H46) 

Just as in (|H.25|) these are obtained from mapping the derivative of the scale factor 
from r to A frame: 



A4o(A) = ^(r (A)), D A (X) = ^ (r(X)) 
ar ar 



(H.47) 



where the functions r(A) and r (X) are given in ( |H.28| ). 

Computing the trace of the extrinsic curvature in the A-frame on the TG solution 
( |H.22p , one obtains, 

S T G G H = SM 3 (3'V^B (X ) 3 D A0 (X ). (H.48) 
Using the expressions above and ( |H.24| ), one has, 



S T g G h = -- D M 6 pXB»{X Ye 



4 f^O d\ 



f X (A) 



^H.49) 



Similarly, the Gibbons-Hawking term evaluated on the BH solution reads, 

S%§ = M 3 (3V 3 B(X ) 3 D A (X ) (8 + 4F(A )) , (H.50) 

which gives, 



S B c H h = ->W 3 5(A ) 4 e-!^ fm ( J + l F(Ao) 



(H.51) 



Just as in the previous subsection, we compute the Gibbons-Hawking contribu- 
tion to the free energy by taking the difference of ( |H.49| ) and ( |H.51| ). Following the 
same steps as outlined above ( |H.42j ), one arrives at the following result: 



5S GH = SBi - Sl G H = ~M*[3V 3 B(X ) 



8 



Y 
X 



25X + -y(A ) . (H.52) 



Evaluating the integrals as before and using the small A asymptotics of Y (see section 
|7.2j ) one arrives at the final expression for the Gibbons-Hawking contribution to the 
free energy: 

Te = - £ (C - Y /2)M 3 V 3 (A if. (H.53) 
Thus, the total free energy is obtained from combining ( |H.45| ) and (|H.53|) as, 



^=-(6C -4y )MV 3 (A 



(H.54) 
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H.7 Fluctuations of A and A in the A-frame: 



In order to compute the gluon condensate, we need to know SX(r) = A(r) — Ao(r), 
near the boundary. This quantity is defined in the r-frame. In the A-frame the 
quantities A and Ao are the same, by definition. However one still finds a non-zero 
value for <5$(r) after a careful change of variables from r to A frame. 

The quantity Ao(r) maps to A(r (A)) where r is defined in ( |H.28| ). Thus, 5X in 
the A-frame is evaluated by: 



S\ = A (r) - A (r) = A - A(r (A)) 



dX 
dr 



Sr(X) = D (X)Sr{X), 



^H.55) 



where D Q is defined in ( H.25|) and Sr is given by, see (|H.28 



8r(X) = r(A) - r (A) 



dX [D(X)- L - D (X) 



(H.56) 



It is straightforward to work out (|H.56|) from ( |H.25|) and (|H.24|) . Let us define, 

B(X) = B (X) (1 + 5A(X)) . (H.57) 

Then, one obtains, 



Xo(X) 3 



5X 



rX Bp (A) 

e x aW . 



(H.58) 



The last term is subleading in the limit A — > 0, so we neglect it for the moment. The 
second termcan be computed from ( |H.24 ) and one finds, 



5A(X) = ~G e-^(b X)-* b - 2 



+ 



Finally the first term follows from (|7.9|) . All in all, one finds, 

5r(X) = ^ K g e~^(b X)- 5b - 2 + --- 



(H.59) 



(H.60) 



where we also used the definition of A in ( H.44| ). Now, 5X follows from ( |H.55 ) as, 

5X = D {X)5r{X) = ^(A£)V^(6 A)- 4fe , (H.61) 

which yields 



5$(r) = ^ = ^Go(A£) 4 (r/£) 4 log(rA). 



^H.62) 



This is the correct coefficient that produces the mathching of the conformal anomaly 
in section HI. 
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One can also compute 5b(r) in the r-frame, from the expressions that we obtained 
above in the A-frame. By definition, 



Sb(r) = b(r)—b (r) = b(r) — b (r )- 



dbo 
dr 



6r = b (X)6A(X)- 



dB 
~d\ 



d {X)Sr{\), (H.63) 



where 5A(X) and Sr(\) are given by (|H.60|) and (|H.59| ) above. We see that the leading 
terms cancel each other out, and one has to take into account the sub leading terms 
in (|H.60|) and ( |H.59 ). This is a straightforward but lengthy computation and best 
carried out by a symbolical evaluation program. We shall only present the result 
here: 

6b(r) = -G (A£) 4 (r/£) 4 . (H.64) 
5 

This is the result in the r-frame. In the u-frame, the coefficient becomes \Gq. 

H.8 Higher order terms in the near horizon expansion 

Demanding regularity at the horizon determines the higher terms in the expansion 
of (|73ll ) and ( pp . One finds, 



( 



16 V V($ h ) V($ h )2 
H.9 Derivation of eq. (7.26) 



Y _ 9 2 

1 64 V V(* h ) 



1. (H.65) 



Here, we compute temperature in the scalar variables. By definition, T is given by 
4vrT = \f{r h )\ = f |. Using the definition of Y in (Q, eqs. (Q, (gj), 

(|7.5|), (|7.4| ) and du/dr = exp(A) we obtain, 

T = — e Ar ^ f x e !^ f (h.66) 



Now, using Y equation of motion ( |7.3| ) we see that the integrand in the last expo- 
nential is a total derivative and can easily be integrated. One finds, 



T 



Til 



Using the UV asymptotics of Y in ( |7.8| ) and 
RHS of flFTBTD as 

Y ° ■r 1 e 4A °-^(&oAo)- 4 1. 



(H.67) 

with b = exp(A) one identifies the 



fH.68) 



The expression in the square brackets defines A, see (|2.19 ). Thus, one finally obtains 
eq. (CT . 
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H.10 Integral representation for the free energy and the energy 

Here, we provide further formulas regarding the integral representation of the free 
energy in section p.2.3| . For the big BH one has, 



^s(Afc) = -4ttM 3 V 3 / ^{Xhj-^dXh, X h < X min . (H.69) 

Note that the two branches are combined in the integral, as b(Xh) is a single valued in 
the entire range Xh G {Ao, oo}. One can also put this equation in various useful forms. 
For example, one can write it in terms of the Yo(Xh) function. For this purpose, we 
make use of the relation between T and Y in ( |7.26| ). Use this in ( |H.69| ), carry out 
the integral and note that Y (oo) = to get: 



f^h A A 



F B (X h ) = A\Ml) d V 3 12 / Y {\ h )-~-d\ h - 4Y (X h ) . (H.70) 

V Joe dX h ) 

Comparison of this equation with ( 7.28|) reveals an alternative expression for the 

energy density: 

f Xh ~ dA 

p(X h ) = A 4 (M£) 3 12 / Y (A fc )-=-dA A (H.71) 
Joo dXh 

H.ll Derivation of eq. (7.38) 

Reorganizing ( |7.6| ) we obtain, 

12 , o . 4 M viiA , 4 rX dX /-i v-2 \?\ 4 rX dX . . 

V = — 1 -X 2 + Y e-3/0 x ~ limes L ~^- x - y )~zL xx . H.72 

Now we use ( |7.3| ) to perform the first integral in the exponent under the limit sign 
and ( [7.4p to obtain: 

V= ^i-X 2 + Y)e-lfo x f lim^L^e" 4 ^. (H.73) 

Now, we can take the limit using the definition of Y ( |7.8| ): 

limF(e)e 4A(e) = Y (A£) 4 . 
e->o 

This yields a relation between enthalpy as a function of A^ and the dilaton potential: 

V(A) = ^(1 - X 2 + Y]*£^e-*£ x * - 4A W. (H.74) 
i z Y (A) 

A remarkable simplification occurs if we take the horizon limit A — > Xh- Using ( 7.21 ) 
and ([7.24 ) one obtains, 

4 rX x dX eV(X h ) e 4A ^ / tt , 

12 Y (A£) 4 K ' 

Finally, using the Bekenstein-Hawking formula for the entropy density 

s = 47rM 3 e 3A(Ah) , 

and ( |735D , we obtain eq. ( fT5S| ). 
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I. High T asymptotics 



We can determine how the various quantities we described above approach their ideal 
gas values for large-T, by studying the next-to-leading order corrections to the AdS 
black- hole. To leading order Y in ( |7.10 ) is determined by the AdS black- hole: 



Y = 7T 



T 1 
A(>lAp 



T 



1 



7rr h 



(1.12 



The subleading corrections to the AdS scale factor, close to the boundary, is presented 
in section ( |5.1|) . From ( |7.8|) and ( |3.10|) , one arrives at the following expression: 



Y (T) = T 4 



1-1 



16 b log(log(T)) | 



3log(T) 9 l og 2 (T) 



(1.2) 



Let us first compute how (p — 3P)/T 4 approaches to at high-T. This determines 
the high-T asymptotics of the gluon condensate. It follows from ( |1.2[ ) and integration 
by parts in (|7.10|) that, 



p — 3p An z 



135 



+ 8 6 



log ( log [ 



V (S) '°s 3 (f ) 



+ 



(1.3) 



This computation can easily be extended to the large-T asymptotics of p, p and s. 
We find: 



P_ 

rp 4 



S 

J*3 



_P_ 

rp 4 



rr 



4vr 2 



45 135 
4vr 2 16tt 2 



45 

37T 2 

45 



135 
4tt 2 



log (f- 



m \og (log (J 

^ iog 2 (£) 

647r 2 6 lQg (log^ 

135 log 2 (£ 

167r 2 6 log (log (j)) 



log £ 



45 



log 2 £ 



+ ■ 



+ ■ 



(1.4) 



(1.5) 



(1.6) 



We note that the pressure, entropy density and the energy density approach their 
ideal gas limits from below as they should. 

It is also useful to derive the high-T asymptotics of the speed of sound. It is 
obtained from eqs. fll^ ) and (|7.33|) : 



4-3 = 1 



1 



32b lo S ( lQ g(i) 



log 3 ( I 



+ 



(1.7) 



- 115 - 



We note that the speed of sound approaches to its ideal gas value 1/3 also from below 
as it should. 

We remark that, all of these expressions are completely independent of the pa- 
rameters of our theory and any modification to the dilaton potential. It follows 
directly from demanding an asymptotically AdS solution dual to the UV of the field 
theory. Moreover, the coefficients of the first terms are even independent of the 
parameters of the gauge theory, i.e. the /3-function coefficients. Thus, we expect 
this form hold universally for any large-A c gauge theory that exhibits logarithmic 
running in the UV. 



J. Analytic solutions 

One can easily obtain analytic solutions to the system, by restricting to the fixed 
points of ( |7.2| ) and (|7.3|). One obvious fixed point of ( [7.31 ) is Y = 0. This takes us 
back to the zero-T analysis which was studied in [^]. In the following we always 
consider the case 7/0. We present our solutions in the domain-wall coordinate 
system, see section |A.2| . They can easily be converted into the conformal frame using 

(U). 

J.l Analytic solutions: zero potential 

Another obvious fixed point of the system is A = const., Y = const and l + Y — X 2 = 
0. From (|7.6| ), we see that this corresponds to vanishing dilaton potential, thus it is 
not very interesting for our purposes regarding holography. Nevertheless, it produces 
the following interesting analytic solution. 

It is straightforward to solve ( |H.1| ), ( |H.2| ) and ( [H.3Q . One finds, for the dilaton: 

A = e*= (d-4X 2 ^)*. (J.l) 
Then the metric can be found as, 

ds 2 = (C x - 4A 2 ^) ^ [dxidx 1 - (c x - 4A 2 ^) ^ dt 2 J + (d - 4A 2 ^) ^~ du 2 . 

(J.2) 

There is a curvature singularity at, 

I ~ AX 2 ' { ^ 6) 

where the dilaton blows up and the metric shrinks to a point (for A 2 < 1) or a line 
(for A 2 > 1). Note that A < in all of our solutions. 

In order to understand the physics of this solution, one has to distinguish these 
two cases. For A 2 > 1 the same point Uq coincides with the event horizon. Therefore 
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there is a curvature singularity at the event horizon Uh = «o where the geometry 
shrinks to a point. We note also that Y > in this case. Therefore, from (|H.3|) we 



see that g (or /) is monotonically decreasing. It decreases from 1 at the boundary 
to at the horizon. 

In the other case, X 2 < 1, there is no event horizon. There is a curvature singu- 
larity at Uq, where the geometry shrinks to a line. Also in this case g is monotonically 
increasing. 

J. 2 Analytic solutions: exponential potential 

A less obvious fixed point of ( [7.2| ) is when X = const, and the dilaton potential is 
exponential: 

V = V (l-X 2 )X~l x . (J.4) 

The proportionality constant will become clear below. 

Also in this case one can find the most general analytic solution to the system. 
This case is of more interest because of the following reasons. We find below that 
in this case Y does not need to be constant or zero and it can be a function of 
$. However, as ( |J .4| ) does not depend on Y, we find that it is a moduli of the 
exponential potential. This fact will allow us to obtain both thermal gas solutions 
and the black-hole solutions to the same potential. Moreover, we note from [2S| that 



the leading IR behavior of the dilaton potential, in most of the confining theories 
is an exponential. Since the confinement- deconfinement phase transition is expected 
to take place in the IR of the theory, ( |J.4|) can be taken a first approximation to 
understand the finite temperature dynamics of the interesting confining theories. 
One first solves (|773| ) to obtain Y as, 

- A"-C 2 ' lJ ' 5) 

where C 2 > and we defined, 

4(1 — X 2 ) /T . 

Note that for a monotonically decreasing g we need Y < 0. As we also require Y — > 
near the boundary, where A — > 0, one should take X 2 < 1, hence a < 0. We stress 
that, this case covers the rest of the physically interesting, constant X solutions, as 
the physical solution in the previous subsection covered the range X 2 > 1. 

As we demonstrate below, the solution (|J.5|) describes a black-hole for Ci > 



and a thermal gas for Ci = 0. To check that indeed Y is a moduli of (|J.4| ), one 
inserts ( pT5| ) in (|7.6|) and finds ( |J.4|) after nice cancelations. 

Both for the black-hole and the thermal gas, one finds the same A behavior for 
the dilaton: 

A = e*=(c 1 -4X 2 ^) 3§r , (J.7) 
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and the scale factor, 

e A = e A °X^. (J.8) 

We note that these are the same as in the previous subsection, as they follow directly 
from the fact that X is constant. 

One finds the location of the horizon by solving for / from ( |7.5| ): 

/ = e 9 = 1 - C 2 X~ a . (J.9) 

We find that indeed / — > 1 on the boundary, (A — > 0) as a < 0. There is an event 
horizon located at (using ( |J.7|) ), 



t4-%- (j - 10) 



The curvature singularity is located at A = oo i.e. , 

Mo _ _Cj_ 
I ~ AX 2 ' 



(J.ll) 



We note that Uh < uq and indeed there is a well-behaved black-hole solution to the 
system. The metric of the black-hole is given by, 

ds 2 = e 2A ° (d -AX 2 -^ ^dxdx 1 - (l -C 2 (Ci -AJPjf^j dt 2 } 

(II 1 — X"^ \ ^ 
l-C 2 (C 1 -4X 2 -)^j du 2 . (J.12) 

The temperature of the black-hole is determined by requiring regularity of the 
Euclidean continuation at Uh- 

1 Aix 

(3= T = \f'{u h )\e^r (J ' 13) 



One finds, 



1 y2 

3~ 



P = \_ 2 X2 • (J-14) 



The physically most interesting case corresponds to the value X = —1/2, see [^S 
Very interestingly, in this case the temperature is only given by the integration 
constant A : 

1 A-Ki 

(3=L = ^L. (j.15) 

Otherwise the temperature is determined by the combination of A Q and C%, namely 
the string tension and the location of the event horizon. 
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The thermal gas solution is found by setting C2 = in (|J.5| ), hence / = 1. The 
dilaton is given again by (|J.7| ) and the metric is, 



ds 2 = e 2A ° - 4X 2 ^) ^ {dxdx 1 + dt 2 } + du 2 . (J.16) 



Here we required the same integration constant for A as the black- hole solution (|J.12|) . 



This is because they should have the same asymptotics at the boundary. Euclidean 
time is compactified with circumference, j3. We note that there is a curvature singu- 
larity ar u Q that is given by ( |J.11|) . It is the same locus as the curvature singularity 
of the black-hole solution-that is cloaked behind the event horizon- resides. 

Computation of the energy of the solutions 

Here we prove that the analytic solutions describe above do not demonstrate a 
Hawking-Page transition. Hence they are not interesting for holographic purposes. 
The action is given by (|2.1| ). One finds that the trace of the intrinsic curvature is 
given by, 

K=^l(8A' + g') (J.17) 

in the domain- wall coordinate system. Thus, the boundary contribution to the action 
becomes, 

S bnd = -M 3 V 3 (3 {e s+4A (8A' + g')}^ , (J.18) 

where u b denotes the regulated boundary of the geometry infinitesimally close to 
—00. 

The bulk contribution to the action, evaluated on the solution can be simplified 

as, 

= 2M 3 V 3 (3 p du^- (fe 4A A') 

= 2MV 3 /3 {f(u s )e 4A ^A'(u s ) - f(u b y A ^A'(u b )} . (J.19) 

Here u s denotes Uq or Uh depending on which appears first. Thus, for the black- hole 
solution u s = Uh, whereas for the thermal gas u s = Uq. 

The first term in ( |J.19|) deserves attention. Clearly it vanishes for the black-hole, 
as f(uh) = by definition. However, it is not a priori clear that it also vanishes for 
the thermal gas. A straightforward computation using (|J.8| ), (|J.7|) and, 

A' = -\\~^ (J.20) 

shows that it indeed vanishes for our physically interesting case X 2 < 1. Therefore, 
one obtains the following total expression for the action from ( J.18| ) and ( |J.19 ) by 
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dropping the first term in ( J.19|) : 



S = -2M 3 V 3 (3e 9{Ub)+4A(Ub) (hA\u h ) + ^g'{u b )\ . (J.21) 

In order to compare the energies of the black-hole and the thermal gas geome- 
tries, we fix the UV asymptotics of the thermal gas geometry by requiring the same 
circumference for the Euclidean time at u b : 



P = z^vTR- (J.22) 

Now, it is straightforward to compute the energy of the geometries. For the 
black- hole flJ.12|) , one finds: 

S BH = -2MV 3 f^j e 4Ao (C 2 (3 + 2X 2 ) - 5A°) . (J.23) 

Here X b is the value of the dilaton on the regulated boundary u b . As a < and 
A — > near the boundary, it is a divergent piece that should be regulated. 
For the thermal gas one finds, using ( J.22| ), 



*,-<;■ -IOMV 3 (^)e 4Ao f Y~ A 6 ) (J.21) 



We note that the divergent terms in (|J.23|) and (|J.24j) cancels in the difference and 
one finds, 



Sbh - S TG = -2MV 3 ( e 4A °C 2 [ZX 2 + I ) . (J.25) 
We note from (|J . 14|) that the temperature is given by, 



l-x 2 



4_ 



e Ao = Y^X~2 C 2- (J-26) 



As (|J.25| ) is always negative, we observe that the BH solution always minimizes the 
action, hence if it exist it is the dominant solution. Therefore, there is no phase 
transition in this geometry. 



References 



[1] J. M. Maldacena, "The large N limit of superconformal field theories and 

supergravity," Adv. Theor. Math. Phys. 2, 231 (1998) [Int. J. Theor. Phys. 38, 1113 
(1999)] |[ArXiv:hep-th/9711200l . 



[2] E. Witten, "Anti-de Sitter space and holography, " Adv. Theor. Math. Phys. 2 (1998) 



253 ||ArXiv:hep-th/9802150 



- 120 - 



[3] S. S. Gubser, I. R. Klebanov and A. M. Polyakov, "Gauge theory correlators from 
-critical string theory," Phys. Lett. B 428 (1998) 105 |ArXiv:hep-th/9802109| . 



non- 



[4] E. Witten, "Anti-de Sitter space, thermal phase transition, and confinement in gauge 
theories," Adv. Theor. Math. Phys. 2 (1998) 505 ||ArXiv:hep-th/9803l"3T|; 



[5] 



J. M. Maldacena and C. Nunez, "Towards the large N limit of pure N = 1 super 
Yang Mills," Phys. Rev. Lett. 86 (2001) 588 ||ArXiv:hep-th/0008001|| ; 
I. R. Klebanov and M. J. Strassler, "Supergravity and a confining gauge theory: 
Duality cascades and xSB-resolution of naked singularities, " JHEP 0008 (2000) 052 
ArXiv:hep-th/0007191|| . 



[6] 



S. W. Hawking, C. J. Hunter and M. Taylor, "Rotation and the AdS/CFT 
correspondence," Phys. Rev. D 59 (1999) 064005 |ArXiv:hep-th/9811056l1 ; 
A. Chamblin, R. Emparan, C. V. Johnson and R. C. Myers, "Charged AdS 
black-holes and catastrophic holography," Phys. Rev. D 60 (1999) 064018 



ArXiv:hep-th/9902170 



M. Cvetic and S. S. Gubser, "Phases of R- charged black-holes, spinning branes and 
strongly coupled gauge theories," JHEP 9904 (1999) 024 ||ArXiv:hep-th/9902195|| ; 
T. Harmark and N. A. Obers, "Thermodynamics of spinning branes and their dual 
field theories, " JHEP 0001 (2000) 008 |ArXiv:hep-th/9910036| ; 
D. Mateos, R. C. Myers and R. M. Thomson, "Holographic phase transitions with 



fundamental matter, " Phys. Rev. Lett. 97, 091601 (2006) |[ArXiv:hep-th/0605046 
D. Mateos, R. C. Myers and R. M. Thomson, "Thermodynamics of the brane," 



JHEP 0705, 067 (2007) |[ArXiv:hep-th/0701132 



[7] O. Aharony, J. Sonnenschein and S. Yankielowicz, "A holographic model of 
deconfinement and chiral symmetry restoration, " Annals Phys. 322 (2007) 1420 
|ArXiv:hep-th/060416il 



[8] C. P. Herzog, "A holographic prediction of the deconfinement temperature, " Phys. 
Rev. Lett. 98, 091601 (2007) ||ArXiv:hep-th/0608151|| . 



[9] N. Evans and E. Threlfall, "The thermal phase transition in a QCD-like holographic 



model, " | ArXiv:0805.0956| |hep-th 



[10] J. Babington, J. Erdmenger, N. J. Evans, Z. Guralnik and I. Kirsch, "Chiral 

symmetry breaking and pions in non- super symmetric gauge / gravity duals, " Phys. 
Rev. D 69, 066007 (2004) ||ArXiv:hep-th/0306018|| . 



[11] T. Sakai and S. Sugimoto, "Low energy hadron physics in holographic QCD," Prog. 
Theor. Phys. 113 (2005) 843 |[ArXiv:hep-th/0412141j| . 

[12] S. Kuperstein and J. Sonnenschein, "Non-critical supergravity (d > 1) and 
holography, " JHEP 0407 (2004) 049 || ArXiv:hep-th /0403254 1| ; 

S. Kuperstein and J. Sonnenschein, "Non- critical, near extremal AdS$ background as 
a holographic laboratory of four dimensional YM theory," JHEP 0411 (2004) 026 



ArXiv:hep-th/0411009 



- 121 - 



I. R. Klebanov and J. M. Maldacena, "Superconformal gauge theories and 
non-critical super strings," Int. J. Mod. Phys. A 19 (2004) 5003 



ArXiv:hep-th/0409133 



F. Bigazzi, R. Casero, A. L. Cotrone, E. Kiritsis and A. Paredes, "Non-critical 
holography and four-dimensional CFT's with fundamentals, " JHEP 0510, 012 (2005) 
ArXiv:hep-th/0505140j ; 

R. Casero, A. Paredes and J. Sonnenschein, "Fundamental matter, meson 
spectroscopy and non-critical string / gauge duality, " JHEP 0601 (2006) 127 



[ArXiv:hep-th/0510110l| ; 



A. Paredes, "On unquenched N = 2 holographic flavor," JHEP 0612 (2006) 032 
ArXiv:hep-th/0610270lt 

C. Csaki and M. Reece, "Toward a systematic holographic QCD: A braneless 



approach, " JHEP 0705 (2007) 062 |[ArXiv:hep-ph/0608266 



[13] J. Polchinski and M. J. Strassler, "Hard scattering and gauge/string duality," Phys. 
Rev. Lett. 88 (2002) 031601 |ArXiv:hep-th/0109174| ; "Hard scattering and 



gauge/string duality, " Phys. Rev. Lett. 88 (2002) 031601 ||ArXiv:hep-th/0109174 



[14] J. Erlich, E. Katz, D. T. Son and M. A. Stephanov, "QCD and a holographic model 
ofhadrons," Phys. Rev. Lett. 95, 261602 (2005) [arXiv:hep-ph/0501128]; L. Da Rold 
and A. Pomarol, "Chiral symmetry breaking from five dimensional spaces, " Nucl. 



Phys. B 721, 79 (2005) |[ArXiv:hep-ph/0501218 



[15] C. A. Ballon Bayona, H. Boschi-Filho, N. R. F. Braga and L. A. Pando Zayas, "On 
a holographic model for confinement / deconfinement," Phys. Rev. D 77, 046002 



(2008) |[ArXiv:0705.1529] [hep-th] . 



[16] A. Karch, E. Katz, D. T. Son and M. A. Stephanov, "Linear confinement and 
AdS/QCD, " Phys. Rev. D 74 (2006) 015005 ||ArXiv:hep-ph/0602229| . 



[17] O. Andreev, "Some Thermodynamic Aspects of Pure Glue, Fuzzy Bags and 



Gauge/String Duality, " Phys. Rev. D 76 (2007) 087702 |ArXiv:0706.3120| |hep-ph 



[18] K. Kajantie, T. Tahkokallio and J. T. Yee, "Thermodynamics of AdS/QCD, " JHEP 
0701 (2007) 019 ||ArXiv:hep-ph/0609254|| . 

[19] J. Adams et al. [STAR Collaboration], "Experimental and theoretical challenges in 
the search for the quark gluon plasma: The STAR collaboration's critical assessment 
of the evidence from RHIC collisions, " Nucl. Phys. A 757 (2005) 102 



ArXiv:nucl-ex/0501009 



B. B. Back et al., "The PHOBOS perspective on discoveries at RHIC," Nucl. Phys. 
A 757 (2005) 28 ||ArXiv:nucl-ex/0410022| |; 

I. Arsene et al. [BRAHMS Collaboration], "Quark gluon plasma and color glass 
condensate at RHIC? The perspective from the BRAHMS experiment," Nucl. Phys. 



A 757 (2005) 1 [ ArXiv:nucl-ex/0410020 



- 122 - 



K. Adcox et al. [PHENIX Collaboration], "Formation of dense partonic matter in 
relativistic nucleus nucleus collisions at RHIC: Experimental evaluation by the 
PHENIX collaboration, " Nucl. Phys. A 757 (2005) 184 ||ArXiv:nucl-ex/0410003|| . 

[20] M. Luzum and P. Romatschke, "Conformal Relativistic Viscous Hydrodynamics: 

Applications to RHIC results at ^/sWJj = 200 Ge V, " Phys. Rev. C 78 (2008) 034915 
ArXiv:0804.4015||nucl-th| . 



[21] G. Policastro, D. T. Son and A. O. Starinets, "The shear viscosity of strongly coupled 
N = 4 super symmetric Yang-Mills plasma, " Phys. Rev. Lett. 87 (2001) 081601 



[ArXiv:hep-th/0104066 



P. Kovtun, D. T. Son and A. O. Starinets, "Viscosity in strongly interacting 
quantum field theories from black-hole physics," Phys. Rev. Lett. 94 (2005) 111601 
[arXiv:hep-th/0405231]. 

[22] E. Shuryak, "Physics of Strongly coupled Quark-Gluon Plasma," 



ArXiv:0807.3033] [hep-ph 



D. T. Son and A. O. Starinets, "Viscosity, Black Holes, and Quantum Field 



Theory, " Ann. Rev. Nucl. Part. Sci. 57 (2007) 95 ArXiv:0704.0240 hep-th 



M. Natsuume, "String theory and quark-gluon plasma, " ArXiv:hep-ph/0701201 



[23] F. Karsch, D. Kharzeev and K. Tuchin, "Universal properties of bulk viscosity near 
the QCD phase transition," Phys. Lett. B 663 (2008) 217 ||ArXiv:0711.0914| | hep-ph || . 

[24] H. B. Meyer, "A calculation of the bulk viscosity in SU(3) gluodynamics," Phys. Rev. 
Lett. 100 (2008) 162001 |[ArXiv:0710.3717] [hep-Iatj| "Energy-momentum tensor 
correlators and viscosity, " |[ArXiv:0809.5202 hcp-lat 



[25] C. P. Herzog, A. Karch, P. Kovtun, C. Kozcaz and L. G. Yaffe, "Energy loss of a 

heavy quark moving through N = 4 supersymmetric Yang-Mills plasma, " JHEP 0607 



(2006) 013; |[ArXiv:hep-th/0605158 



[26] H. Liu, K. Rajagopal and U. A. Wiedemann, "Calculating the jet quenching 
parameter from AdS/CFT, " Phys. Rev. Lett. 97 (2006) 182301 
[ArXiv:hep-ph/0605178]| . 



[27] S. S. Gubser, "Drag force in AdS/CFT," Phys. Rev. D 74 (2006) 126005 
|[ArXiv:hep-th/0605T82]l 

[28] U. Gursoy and E. Kiritsis, "Exploring improved holographic theories for QCD: Part 



I, " JHEP 0802 (2008) 032 |ArXiv:0707.1324 ||hep-th 



U. Gursoy, E. Kiritsis and F. Nitti, "Exploring improved holographic theories for 



QCD: Part II, " JHEP 0802 (2008) 019 [ArXiv:0707.1349] [hep-th 



[29] S. S. Gubser, "Curvature singularities: The good, the bad, and the naked," Adv. 



Theor. Math. Phys. 4 (2000) 679 ArXiv:hep-th/0002160 



- 123 - 



[30] U. Gursoy, E. Kiritsis, L. Mazzanti and F. Nitti, "Deconfinement and Gluon Plasma 
Dynamics in Improved Holographic QCD," Phys. Rev. Lett. 101, (2008) 181601 



ArXiv:0804.0899l hep-th 



[31] S. S. Gubser and A. Nellore, "Mimicking the QCD equation of state with a dual 
black-hole, " Phys. Rev. D 78 (2008) 086007 || ArXiv:0804.0434| [hep-th] . 

[32] G. Boyd, J. Engels, F. Karsch, E. Laermann, C. Legeland, M. Lutgemeier and 
B. Petersson, "Thermodynamics of SU(3) Lattice Gauge Theory," Nucl. Phys. B 
469, 419 (1996) ||ArXiv:hep-lat/9602007|| . 

[33] B. Lucini, M. Teper and U. Wenger, "Properties of the deconfining phase transition 
in SU(N) gauge theories," JHEP 0502, 033 (2005) ||ArXiv:hep-lat/0502003|. 



[34] U. Gursoy, E. Kiritsis, G Michalogeorgiakis and F. Nitti, in preparation. 

[35] H A. Chamblin and H. S. Reall, "Dynamic dilatonic domain walls, " Nucl. Phys. B 
562, 133 (1999) ||ArXiv:hep-th/9903225| . 



T. Hertog and K. Maeda, "Stability and thermodynamics of AdS black-holes with 



scalar hair," Phys. Rev. D 71, 024001 (2005) |[ArXiv:hep-th/0409314l 



[36] U. Gursoy, E. Kiritsis, L. Mazzanti and F. Nitti, "Improved Holographic Yang-Mills 
at Finite Temperature: Comparison with Data," [ArXiv: 0903. 2859] [hep-th] . 



[37] T. Hertog, "Towards a novel no-hair theorem for black-holes, " Phys. Rev. D 74, 
084008 (2006) |[ArXiv:gr-qc/0608075] 



[38] S. W. Hawking and D. N. Page, "Thermodynamics Of Black Holes In Anti-De Sitter 
Space," Commun. Math. Phys. 87, 577 (1983). 

[39] P. K. Townsend, "Positive Energy And The Scalar Potential In Higher Dimensional 
(Super) Gravity Theories," Phys. Lett. B 148, 55 (1984). 

[40] K. Skenderis and P. K. Townsend, "Gravitational stability and renormalization- group 
flow," Phys. Lett. B 468, 46 (1999) ||ArXiv:hep-th/9909070|| . 
O. DeWolfe, D. Z. Freedman, S. S. Gubser and A. Karch, "Modeling the fifth 
dimension with scalars and gravity, " Phys. Rev. D 62, 046008 (2000) 
||ArXiv:hep-th/9909i34l 

[41] E. Vicari and H. Panagopoulos, "Theta dependence of SU(N) gauge theories in the 
presence of a topological term, " || ArXiv:0803.1593 hep-th || . 

[42] A. Kehagias and K. Sfetsos, "On running couplings in gauge theories from type-IIB 
supergravity," Phys. Lett. B 454, 270 (1999) |ArXiv:hep-th/9902125| . 



[43] S. S. Gubser, "Dilaton-driven confinement, " [ArXiv:hep-th/9902155 



[44] Y. Kinar, E. Schreiber and J. Sonnenschein, "Q anti-Q potential from strings in 
curved spacetime: Classical results, " Nucl. Phys. B 566, 103 (2000) 
|ArXiv:hep-th/981lT92l 



- 124 - 



[45] M. Bianchi, D. Z. Freedman and K. Skenderis, "Holographic Renormalization, " Nucl. 
Phys. B 631, 159 (2002) ||ArXiv:hep-th/0112119|| . 

I. Papadimitriou and K. Skenderis, "Thermodynamics of asymptotically locally AdS 
spacetimes, " JHEP 0508, 004 (2005) ||ArXiv:hep-th/0505190|| . 

[46] S. W. Hawking and G. T. Horowitz, "The Gravitational Hamiltonian, action, entropy 
and surface terms," Class. Quant. Grav. 13, 1487 (1996) |[ArXiv:gr-qc/9501014] 



[47] I. R. Klebanov and E. Witten, "AdS/CFT correspondence and symmetry breaking," 
Nucl. Phys. B 556, 89 (1999) |[ArXiv:hep-th/9905104]| . 

[48] S. S. Gubser, I. R. Klebanov and A. W. Peet, "Entropy and Temperature of Black 
3-Branes," Phys. Rev. D 54 (1996) 3915 |[ArXiv:hep-th/9602135]| . 

[49] U. Gursoy and I. Papadimitriou, to appear. 

[50] A. Buchel and J. T. Liu, "Universality of the shear viscosity in supergravity," Phys. 
Rev. Lett. 93 (2004) 090602 ||ArXiv:hep-th/0311175|| . 

[51] D. Anninos and G. Pastras, "Thermodynamics of the Maxwell-Gauss-Bonnet anti-de 
Sitter Black Hole with Higher Derivative Gauge Corrections, " arXiv:0807.3478 
[hep-th] . 



- 125 - 



